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The method of the delayed impulses has been applied to the study of the instability of the 
meson. The decay curve of the mesons at rest can be traced, as in the case of a radioactive 
element, and an evaluation of the lifetime given. 


METHOD 


HE technique of the delayed impulses from 
Geiger-Miiller counters as devised by R. 
Maze, has been first used by us for the detection 
of very broad A showers.'! The delay was intro- 
duced between one of the counters and the 
nearby coincidence set, in order to compensate 
the inevitable delays occurring in the transmis- 
sion of impulses from far distant counters. 

In the study of the decay of the meson at rest 
we have to consider the following pair of events: 
First, the passage of a meson as detected by a 
counter telescope (a and 6 in Fig. 1) and its 
subsequent stopping in an appropriate absorber 
A, then the emission by the stopped meson of its 
decay electron detected by a third counter (c in 
Fig. 1). The second event is separated from the 
first by the duration of the life at rest of the 
particular meson considered. The two impulses 
coming from the coincidence a, b and from the 
counter c can be brought in coincidence only if 
the first is delayed by the proper amount of time 
before entering the circuit set. So the number of 
coincidences recorded with a particular delay 
will measure the number of mesons which have 
lived, at rest, just that length of time. If we 
measure the number of coincidences per hour for 

* Research associate at the University of Chicago. 


'P. Auger, R. Maze, P. Ehrenfest, and A. Freon, J. de 
Phys. [7] 10, 39 (1939). 


increasing values of the time delay, we shall ex- 
pect to find the exponential decrease which 
characterizes the radioactive decay curves. 


EXPERIMENTAL ARRANGEMENT 


The experiments described here were begun 
in 1940 in Paris, but because of the circum- 
stances definite results have been obtained only 
at the end of 1941. 

A beam of mesons was defined by counters a 
and } which were separated by a lead block 10 
cm thick (B in Fig. 1). The beam fell then on an 
aluminum absorber 4 cm thick covering the 
whole solid angle. Aluminum was chosen in 
preference to lead because of its relatively greater 
transparence for electrons, under thicknesses 
which are equivalent in stopping power for 
mesons. The decay electrons issuing from this 
absorber were detected by a side counter c, placed 
outside of the beam defined by a and b. All the 
counters were of the so-called slow type (argon 
+hydrogen filling) and had an effective area of 
about 100 cm*. 

Counters a and 6 were connected with a double 
coincidence set, and the impulse issuing from this 
set, as well as the impulses coming from c, passed 
by two separate time delay circuits before being 
sent to the final coincidence set. Any time differ- 
ence ranging from 0 to 10-5 sec. could be intro- 
duced between them. A difference of one division 
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Fic. 1. a, 6, counter telescope; C, first coincidence stage; 
¢, side counter; D;, Dy, time delay circuits; F, final 
coincidence stage. 


on the dials of the variable condensers of the 
delay circuits corresponded to a shift of (1.5 
+0.05)10-7 sec. 

The resolving power of the final coincidence set 
has been measured by the study of the coinci- 
dences obtained when artificial simultaneous im- 
pulses are sent into the delay circuits. No coinci- 
dences were obtained when the difference of 
reading on the dials was of more than 7 divisions; 
with smaller differences all the impulses gave 
coincidences. From the value of the division 
(measured with an oscillograph) the resolving 
time was calculated to 1.05 X10~° sec. 

If instead of artificial impulses the impulses 
coming from the counters a and bd are sent into 
the delay circuits, the distribution curve of the 
number of registered coincidences with varying 
time differences is not square but bell-shaped 
because of small variations in the latent period 
separating the real passage of the particles and 
the beginning of the counter impulses. 


MEASUREMENTS 


The number of mesons passing through the 
telescope was 500 per hour and the number of 
triple coincidences a b c without any time delay, 
(coincidences which are due to showers) was of 
the order of 10 per hour. This last number 
suffered no measurable change when the ab- 
sorber A was set in place, or taken away. 


AUGER, MAZE, AND CHAMINADE 


Three series of measurements were made, with 
differences of 14, 20, and 30 divisions in the 
readings of the delay circuits dials, the coinci- 
dences from a b being of course retarded against 
the impulses from c. The measurements were 
made in each case with and without the alumi- 
num absorber and are given in mean number per 
ten hours. The differences are taken as propor- 
tional to the numbers of mesons having presented 
lifetimes at rest within the intervals defined by 
the delays and by the resolving time of the final 


stage. 


Dela divisions 14 20 30 
seconds 2.1X10-§ 3.0X10-* 4.5x10-* 
With absorber 0.8 0.73 0.35 
Without absorber 0.17 0.29 0.2 
Difference 0.6340.2 0.4440.2 0.15+0.1 


From these numbers we can estimate the mean 
lifetime of the mesons at rest (time during which 
a meson has a probability of one-half for decay) 
as equal to 1.3 microseconds, with an accuracy of 
about 50 percent. This result is in better agree- 
ment with the results of Nielsen, Ryerson, Nord- 
heim, and Morgan? and with those of Rasetti* 
than with the measurements of Rossi, Greisen, 
Stearns, Froman, and Koontz:4 


Note added in proof: In the August issue of this Review 
M. de Souza Santos has reported some results on the life- 
time of the meson, obtained with a method which seems 
very similar to ours. However the curve he obtains for 
the number of coincidences as a function of the time delay 
introduced in the transmission of the meson telescope im- 
pulse does not show an exponential decrease as in our 
experiments, but a Gaussian shape around a maximum. 
The counting rate, as indicated by the low value of the 
statistical errors, is surprisingly high for the final resolving 
power claimed, even though the maximum value (for a 
delay of 5 microseconds) has not been effectively reached. 

We should like to point out the similarity of M. de Souza 
Santos’s meson decay curve with the bell-shaped curve 
which we have obtained when testing the resolving power 
of our coincidence outfits by the use of delay circuits. 
Registering the coincidences between the impulses of two 
counters in a vertical telescope, after their passage through 
two delay circuits, a symmetrical curve was found, with 
a high maximum (almost 100 percent of the real number 
of particles) for a certain adjustment of the delay circuits. 
We considered this adjustment as corresponding to equal 
delays on both pulses. 


2 Nielsen, Ryerson, Nordheim, and Morgan, Phys. Rev. 
57, 158 (1940). 

F. Rasetti, Phys. Rev. 60, 198 (1941). 
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A gamma-ray spectrograph of high resolving power has been employed in measuring the 
quantum energies and relative intensities of the gamma-rays from Na™. The energies thus 
determined are 0.84, 1.31, 1.66, and 2.90 Mev with relative intensities 0.28, 0.41, 0.45, and 
1.00. These results indicate excitation levels in the Mg™ residual nucleus at 1.3, 2.9, and 3.7 
Mev in partial agreement with experiments on proton scattering by magnesium, but in dis- 
agreement with the level schemes suggested by previous gamma-ray measurements. 

Measurements were also made on the gamma-radiation of the thorium active deposit. 
Close agreement with previous results was considered adequate confirmation of the calibration 


of the spectrograph. 


INTRODUCTION 


HE energies and relative intensities of the 

gamma-rays emitted from Na* have been 
studied by several investigators.'~’ Their results 
do not agree very well, and the presence of 
gamma-radiation in the neighborhood of 2 Mev 
has been in dispute. This paper contains an ac- 
count of new measurements of these gamma-ray 
energies and intensities, and a level scheme 
for Mg*™ is suggested. The gamma-rays from 
Th(C’+C”) have also been observed as a check 
on the accuracy of the method. 


THE SPECTROGRAPH 


The energies of the Compton recoil electrons 
arising from gamma-rays were measured by semi- 
circular focusing in a magnetic field. A diagram 
of the spectrograph is shown in Fig. 1. It is very 
similar to that described by Curran, Dee, and 
Strothers.5 The depth of the magnet box is 3 cm, 
and its inside diameter is 18 cm. Compton elec- 
trons are ejected from the bottom of the cup C 
and in order to be counted must follow a circular 
path of radius 5.50 cm. The cup is so inclined that 
the electrons counted emerge at an angle of 
about 20° with the bottom of the cup, thus mak- 
ing the effective width of the source small. The 

$j. paieets Richardson and F. N. D. Kurie, Phys. Rev- 
50, 999 (1936). 

2]. Reginald Phys. 53, 124 (1938). 

3S. Kikuchi, Y. Watase, E. Takeda, and S. 
Yamaguchi, Proc. Phys. Math. Soc. Ja n 21, 260 (1939). 

*S. Kikuchi, Y. Watase, J . Itoh, pa Takeda, and S. 
Yamaguchi, Proc. Phys. Mtb: Soc. Japan 21, 381 (1939). 

5S. C. Curran, P. I. Dee, and J. E. Strothers, Proc. Roy. 
Soc. A175, 546 (1940). 

tJ. Itoh, Proc. Lig = Math. Soc. Japan 23, 605 (1941). 


G. Elliot, M. Deutsch, and A. Roberts, Phys. Rev. 
61, 99 (1942). 


Geiger-Mueller counters 7}, 7:2, 7; and the entire 
box are filled with an argon-alcohol mixture at 
10 cm Hg, the ratio of the partial pressures being 
about 9 : 1. The counters are 2.0 cm long and are 
mounted in holes which are bored in a block of 
Lucite. Their diameters are 1.0, 1.5, and 2.0 cm, 
respectively. Triple coincidences are recorded. 
The copper cases, of thicknesses 0.001 cm, are 
contiguous but the wire potentials are inde- 
pendently variable. The varying diameters allow 
for divergence of the beam of electrons after the 
focusing at the slit S;, just beneath the counters. 
The width of this slit is about 3.0 mm. The sides 
of this slit are strips of brass 2 cm in length, 1 cm 
in width, and sufficiently thick to stop an electron 
of 5-Mev energy, thus preventing the counting 
of particles not focused on the slit proper. Three 
other aluminum slits S2, S3, and S,are also placed 


Fic. 1. The spectrograph. 
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Fic. 2. ante, N, of Compton recoils from gamma-rays 
from Th(C’+C”) as a function of Hp. 


in the box for collimating purposes. The center 
slit S3, the defining one, is 1.8 cm wide, and the 
lower edge is 5.3 cm above the bottom of the box, 
thus making recesses into which electrons not 
entering the counters may drop. The lead block 
B employed to shield the counters from gamma- 
rays coming directly from the source, is 5.8 cm 
in thickness. The Hp spread allowed by the source 
width and the slit system is estimated to be 
about 5 percent. The maximum measurable ki- 
netic energy which a beta-particle may have in 
the field of the Weiss magnet used is about 12 
Mev. A typical Rossi circuit and the usual re- 
cording meter and thyratron were employed in 
counting the triple coincidences. 

The background of the counters was found to 
be constant for both Th(C’+C”) and Na* over 
a wide Hp interval above the end point of the 
hardest gamma-ray measured in each case. This 
background was also equal to the background 
at Hp=0. The background at intermediate points 
was therefore obtained by extrapolation. Since 
the background was small as compared to the 
total number of counts observed, any error in 
calculation of the relative intensities introduced 
by this procedure would be small. 

There are four factors for which corrections 


é 
He «10 OERSTED-CM 
Fic. 3. Momentum distribution of Compton recoils of 
gamma-rays from Th(C’+C”) 
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must be made: (1) increase with increasing H of 
the value of the H/p interval over which electrons 
are counted, (2) variation of the Compton scat- 
tering coefficient with energy, (3) dependence of 
the emission of an electron upon its range in the 
thick aluminum target, (4) variation with energy 
of the absorption of electrons in the walls of the 
counters. This last correction is small for electron 
energies greater than 1 Mev. Plotting N/Hp 
against Hp serves as a correction for (1). Factors 
(2) and (3) complicate the determination of the 
relative intensities of the various components of 
the gamma-ray spectra. They may be most satis- 
factorily dealt with, however, by evaluating the 
relative intensities with the equation 


Area under N/Hp plot 


(electron range in Al) (Compton coeff.) 


Rel. Int. = 


A rough correction for (4) is obtained from known 
absorption curves.® 


I. Gamma-rays from Th(C’+C”). 


Energy in Mev 0.73+.02 1.584+.03 1.77+.04 2.664.05 3.32+.10 
Rel. Int. 0.14 0.10 0.05 1.00 0.09 


8 E. Madgwick, Proc. Camb. Phil. Soc. 23, 970 (1927). 
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Figure 2 is a plot of N against Hp and Fig. 3 
one of N/Hp against Hp. The peaks of the out- 
standing components of the gamma-ray spectrum 
are clearly evident in the latter case. Their ener- 
gies and relative intensities are given in Table I. 
Ellis has reported gamma-rays from Th(C’+C”’) 
at 1.63, 1.80, and 2.62 Mev with intensities 0.08, 
0.04, and 1.00. Weak gamma-rays of quantum 
energy about 3.3 Mev have been reported by 
several authors. The intensity reported here is 
somewhat higher than the most recent report by 
Itoh and Watase.’® The intensity of the 0.73 
Mev quantum agrees closely with that given by 
these authors.'° 

The clearly separated peaks of the gamma-rays 
at 1.58 and 1.77 Mev demonstrate the high re- 
solving power of the spectrograph. The sharp 
rise to a peak and the equally sharp drop to a 
clearly defined end point from which the gamma- 
ray energy may be obtained is to be expected 
with a spectrograph of this type from the nature 
of the Compton process. It is further interesting 
to note the shape of the peaks of the gamma-rays 
of lower quantum energy and lower intensity in 
the neighborhood of the intense 2.66-Mev peak. 
This graph is very useful in analyzing the curves 
obtained from other radioactive sources. 


RADIO-SODIUM 


The gamma-rays of Na™ were initially investi- 
gated by Richardson and Kurie,' and again by 
Richardson,’ the expansion chamber and mag- 
netic field being employed in both instances. 
Quantum energies of 1.01, 2.04, and 3.00 Mev 
were reported, and the intensities of the two 
softer components were estimated to be about 
equal. Since the sum of their energies was ap- 
proximately equal to the energy of the third 
component, it was assumed that levels of 2 and 
3 Mev or 1 and 3 Mev existed in the Mg™ 
product nucleus. 


Il. Gamma-rays from Na™. 


Energy in Mev 0.84 +.02 1.31 +.03 +.03 ry +.06 


1 
Rel. Int. 0.28 0.41 0.45 00 


1935), D. Ellis, Rep. Internat. Conf. on Phys. (London, 
Py Itoh and Y. Watase, Proc. Phys. Math. Soc. Japan 
23, 142 (1941). 
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Fic. 4. Momentum distribution of Compton recoils of 
gamma-rays from Na*™. 


Upon the suggestion of Feather and Dun- 
worth," this level scheme was revised, the new 
assumptions being that gamma-rays of 1 and 3 
Mev were emitted in cascade and two rays of 
quantum energies about 2 Mev were also emitted 
in a cascade process. Kikuchi ef al.* reported 
gamma-rays of about equal intensities with 
quantum energies 1.49+.05 and 2.97+.07 Mev, 
thus establishing the 4-Mev level suggested by 
Feather and Dunworth at 4.46 Mev. Later work 
by Kikuchi et al.‘ resulted in the values 1.55+.05 
and 2.97+.05 Mev. In both experiments, weak 
radiation at about 0.8 Mev was also observed. 

Curran, Dee, and Strothers® have reported the 
energies 1.46, 2.00, and 3.03 Mev with relative 
intensities 1.17, 0.27, and 1.00. The most recent 
measurements®? have given indication of only 
two gamma-rays, having energies 1.38+.02 and 
about 2.80 Mev. These latter authors may have 
failed to detect any other radiation because of 
low intensities. 

In this experiment, a sodium chloride crystal, 
activated by a bombardment of several hours by 
a deuteron current of 0.3 microampere at 2 
million volts, was used as a source, -the half-life 


1 N. Feather and J. V. Dunworth, Proc. Camb. Phil. 
Soc. 34, 442 (1938). 
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Fic. 5. Level scheme for Mg*. 


being 14.8 hr. Radio-chlorine of 37-min. half-life 
is formed from the deuteron bombardment of 
Cl*? (relative abundance 25 percent). Although 
the Cl-d-p reaction is known to be very small as 
compared to the Na-d-p reaction at this voltage, 
and although measurements were not made until 
after any radio-chlorine present had decayed by 
about 5 half-periods, a potassium chloride crystal 
was bombarded under the same conditions for the 
sake of comparison. Contributions from radio- 
chlorine were found to lie well within the statis- 
tical probable error of the points of the curve for 
Na*™*. Above Hp=4000, about 800 coincidences 
were counted at each point. Below that value of 
Hp about 550 coincidences were recorded at each 
point. 

The curve for Na™*, V/Hp against Hp, obtained 
with the spectrograph just described, is shown in 
Fig. 4. The results are tabulated in Table II. 

These quantum energies and relative intensi- 
ties suggest that de-excitation of the Mg*‘ nucleus 
after beta-emission by Na*‘ may occur with the 
emission of a single quantum of energy about 2.90 
Mev or with the emission of quanta of energies 
1.66 and 1.31 Mev in cascade. It would seem then 
that the disintegration energy of Na* is 1 Mev 
less than that suggested by Feather and Dun- 
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worth." The level scheme is similar to that 
initially suggested by Richardson and Kurie,! 
though theirs was based upon different quantum 
energies. 

The quantum of energy 0.84 Mev is probably 
the weak line first observed by Kikuchi et al.,*4 
and is, as they have pointed out, indicative of 
the emission of a complex beta-ray spectrum by 
Na*™, contrary to the results of Feather and Dun- 
worth" and of Lawson,” who found the beta-ray 
spectrum to be simple and of maximum energy 
1.4 Mev. The results of Table II indicate that the 
soft component of the beta-ray spectrum is about 
one-fifth as probable as the more intense distri- 
bution of higher maximum energy. A complete 
level scheme is given in Fig. 5. By means of ex- 
periments on the inelastic scattering of protons 
by magnesium, Wilkins" has assigned excitation 
energies to Mg” of 1.37, 2.80, and 4.07 Mev in 
partial agreement with the results reported in 
this paper. 

It seems possible that other investigators have 
over-estimated the intensity of the gamma-ray 
at 1.31 Mev because of insufficient resolution and 
the unsuspected presence of the gamma-ray at 
1.66 Mev. Attributing all of the areas of the 0.84- 
and 1.31-Mev quanta and about one-half the 
area of the 1.66-Mev quantum to a hypothetical 
gamma-ray of energy about 1.45 Mev, and apply- 
ing all corrections for that energy, an intensity 
about equal to that of the 2.90-Mev quantum is 
obtained, thus offering a possible explanation for 
previous interpretations. 
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The multiple scattering of electrons from B" has been studied by means of G-M coincidence 
counters in two scatterers of equal N7Z* where N is number of atoms per cc of scatterer, 
T and Z are the thickness and atomic number of scatterer. The scatterers were of Pb and C, 
respectively, and had a theoretical scattering power which differed by 20 percent. Comparison 
of the arithmetic mean angles of scattering showed that their ratio was nearly the same as 


given by the Williams theory. 


INTRODUCTION 


REVIOUS measurements of multiple scat- 

tering of electrons of several Mev energy'* 
indicated that a discrepancy existed between 
experiment and theory*® which was greater for 
lead than for elements of low atomic number 
such as carbon and aluminum. Since this 
discrepancy seemed difficult to account for 
theoretically, and since previous measurements 
were made with a Wilson cloud chamber, it 
seemed advisable to endeavor to detect it by a 
different method. 


EXPERIMENTAL DETAILS 


Since the theory of multiple scattering has 
been given in terms of the projection of the 
scattering angles upon a plane, a two-dimensional 
experiment was chosen in which G-M coinci- 
dence counters were used. The experimental 
arrangement can be seen in Fig. 1: A boron 
target was bombarded with deuterons acceler- 
ated by an a.c. tube operating at 350-kv peak 
voltage. A portion of the B™ electrons emerging 
from the front of the target was collimated by 
two lead slits after they traversed the 15-mil 
aluminum target cup. The electrons were then 
scattered and recorded by two cylindrical 
counters placed in a light wooden frame which 
could be moved about an axis through the 


* Now at South Dakota State School of Mines, Rapid 
City, South Dakota. 

1W. A. Fowler, Phys. Rev. ” 773 (1938). 

Cc. W. Sheppard and W. A. Fowler, Phys. Rev. 57, 
273 (1940). 

*Oleson, Chao, Halpern, and Crane, Phys. Rev. 56, 
482 and 1171 (1939). 

*E. J. Williams, Proc. rg Soc. 169, 531 (1939). 

5S. Goudsmit and J. L . Saunderson, Phys. Rev. 57, 
24 (1940) ; 58, 36 (1940). 


scatterers. Behind the target was placed a single 
counter to be used’ as a monitor. This counter 
was heavily shielded but so situated that it 
received a part of the electrons which traversed 
the 10-mil copper target and emerged from the 
rear of the target cup. This arrangement gave a 
satisfactory check on the bombarding intensity 
during the experiment. The coincidence counters 
were of the tubular type. The glass walls had a 
surface density of 16 g/cm?*. The cylirfders were 
13” in diameter, 9” long, and were made of 
8-mil copper oxidized with nitrogen peroxide. 
They were filled with argon and alcohol, their 
pulses were found to be fairly uniform, and their 
plateaus quite broad. The singles counter was 
similar in construction except that it had a 
cylinder only 2” long. The natural background 
of the counters was measured before and after 


ree 


Leao 


Fic. 1. Experimental arrangement employed. 


the experiment. For the two large counters it 
was 111 and 163 counts per minute, respectively, 
at the beginning and 98 and 165 at the end. 
For the singles counter it was 59 before and 69 
after. Since high voltage stabilization was not 
employed it was necessary to check the counter 
voltages every two minutes during operation. 
The coincidence counters were coupled by a 
Rossi coincidence circuit to a pulse-lengthener 
and recorder. The counters were shielded by 3” 
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of lead and had a relatively low coincidence 
background when the tube was in operation. 
In all cases the counting rate was slow enough 
to minimize skipping effects. The singles counter 
was coupled through a Neher-Harper stage to a 

- scale-of-eight circuit. This was necessary because 
of the relatively high background to which it 
was exposed during operation which necessitated 
a proportionately higher counting rate. 


002 
= 12° 
& 


Fic. 2. Experimental results, arbitrary scale. 


SCATTERERS 


Two scatterers were measured of equal NTZ? 
(No. of atoms per cc, thickness and atomic 
number of scatterer), one of rolled lead and the 
other of Acheson graphite. The surface densities 
were determined by weighing before and after 
the scattering measurements were made. Both 
scatterers were 1” wide and 5” long. The 
important theoretical quantities for them are 
given in Table I. After correction of Z for extra- 
nuclear electron scattering and inclusion of the 
logarithmic factor the theoretical scattering 
differs by 20 percent. The notation which has 
been used in Table I is the same as has been 
used elsewhere.” 


EXPERIMENTAL PROCEDURE 


The counting equipment was arranged so that 
it could be stopped and started by means of a 
master switch. In addition a magnetic shutter 
was provided with which the bombarding beam 
could be interrupted in order to obtain back- 
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ground counting rates. The scatterers were 
interchanged constantly during the experiment. 
An estimate of the angle subtended by the 
sensitive region of the counters is rather difficult. 
In practice five angular intervals of 7° were 
taken in the region from 0° to 28°. Any over- 
lapping effect which might occur will affect both 
scatterers by the same factor. In order to 
compensate for the rather large difference in 
energy loss in the two scatterers end pieces were 
attached to them and each was placed between 
the counters during the measurements on the 
companion scatterer. The total number of 
coincidence counts recorded in any angular 
interval varied from 500 to 100 depending upon 
the scattered intensity. The monitor counts 
averaged about 15,000. The data were averaged 
and background rates subtracted. Mean errors 
were calculated and combined to obtain the 
mean error in the result. These errors were 
obtained both from the total number of counts 
in each case and also from the computed re- 
siduals, thus affording a good indication of 
consistent counter behavior. In only two cases 
was the latter error greater than the former. 
The final results were expressed as ratios of 
coincidence counter intensity to that of the 
monitor counter. These results are plotted in 
Fig. 2. A measure of the resolving power of the 
experiment was obtained by measuring the 
counting rates with no scatterer present. These 
data are shown in Fig. 3. 


TABLE I. Theoretical quantities for scatterers used. 


Ws Wa We: Wam 
in in in in 
Mev- Mev- Mev- Mev- 
Surface density ¢ eff deg. M deg. deg. deg. 
Pb 0.177 g/cm? 82.3 24.7 375 84 210 76 
= 1.28 g/cm? 6.49 27.0 2390 99 270 92 
RESULTS 


In Fig. 2 Gaussian curves, normalized to the 
same arbitrary scale, have been fitted to the 
data. The mean angle for lead is seen to be 
10° and for carbon 12°. The mean angle for the 
curve in Fig. 3 is approximately 7°, so it is 
necessary to compensate for the effect of poor 
resolving power. Since the mean angle for the 
resulting distribution obtained by measuring a 
Gaussian distribution with an instrument having 
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a Gaussian error is the square root of the sum of 
the squares of the mean angle of the distribution 
and of the error of the instrument, this compen- 
sation can be made. The resulting values of am 
are therefore closer to 10° and 7°, respectively. 
Their ratio is thus 0.7 whereas the theoretical 
quantities have a ratio of 0.8. 


DISCUSSION 


Energy loss effects are important in the 
experiment, and the question arises as to 
whether the placing of the carbon scatterer 
between the counters when the lead is scattering, 
and vice versa actually produces proper compen- 
sation for the reduction in intensity by the 
carbon due to energy loss when it is scattering. 
A consideration of the energy loss of electrons 
in the region above 5 Mev shows that the 
variation with energy loss by ionization is quite 
slow. Radiative losses in the carbon are still 
negligible. Straggling effects due to multiple 
scattering in the carbon must be considered. 
The best justification for the method of compen- 
sation used can be obtained by considering the 
total reduction in forward intensity from 0° to 
28° when the carbon is scattering compared to 
that when the lead is scattering. This reduction 
was found to be a factor of 0.54. If we compare 
this with the reduction in intensity in the case 
of the scattering of lead effected by placing the 
carbon scatterer between the counters we find 
that the latter reduction is by a factor of 0.53. 

The fitting of a Gaussian curve to the measured 
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INTENSITY 


Fic. 3. Intensity distribution of electrons with 
no scatterer present. 


i0 20 30° 40 
oc IN DEGREES 


Fic. 4. Distribution in scattered intensity to be expected 
on the basis of a sharp energy cut-off in the B® spectrum 
at 5 Mev. Dotted curve is Gaussian having a computed 
mean energy of 6.2 Mev. 


scattering distribution might appear to ignore 
the fact that the electrons from B™ have a wide 
energy distribution. In considering this effect it 
is necessary to make an estimate of the minimum 
energy with which electrons can traverse the 
system and be recorded in the last counter. 
This estimate is complicated by the effects of 
multiple scattering of the diffusion type and 
cannot be made with precision. For purposes of 
discussion an estimate of 5 Mev can be made for 
this cut-off energy. With this value, the electron 
spectrum® of B” was approximated by a poly- 
nomial. Figure 4 shows the distribution for 
carbon obtained by integrating according to this 
method. Since the mean angle is inversely 
proportional to the energy in the Williams theory 
the a», for the approximate Gaussian curve 
should be the r.m.s. value which can be obtained 
by computing a mean energy for the incident 
electrons. This can be obtained by integrating 
the mean square reciprocal of the energy. In the 
case at hand this was found to be 6.2 Mev. 
A Gaussian curve for carbon having a value of 
Qm computed from this mean energy is also 
shown in Fig. 4 and can be seen to fit within 
the accuracy of this experiment. The inclusion 
of the second term in a Gram-Charlier series 
expansion was found to give a very close fit 
indeed, but its use was not considered necessary 
in this case. 

The estimate of 6.2 Mev for the mean energy 
of the electrons is subject to too much uncer- 
tainty, and so the values of a» calculated from 


*D. S. Bayley and H. R. Crane, Phys. Rev. 52, 604 
(1937). 
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this cannot be depended upon. However, the 
ratio of a» for carbon and for lead should be 
correct within 10 percent at least. On the basis 
of this result it seems doubtful that any large 
systematic deviation from the Williams theory 
can occur with increasing atomic number. 
Recent measurements by others also support 
this conclusion.” * However, in order to examine 


7Oleson, Chao, and Crane, Phys. Rev. 60, 378 (1941). 
§L. A. Kulchitzky and G. D. Latyshev, Phys. Rev. 61, 
254 (1942). 
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the seeming discrepancy between the cloud- 
chamber and counter results it should be 
possible to perform an experiment in which 
measurements by the two methods could be 
made simultaneously. 

The author wishes to express his appreciation 
to Mr. Sylvan Rubin, Mr. Eldred Hough, and 
Mr. R. E. Logan for assistance with the experi- 
mental work, and to Professor W. A. Fowler 
and Professor C. C. Lauritsen for advice and 


assistance. 
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The errors in coincidence measurements with self- 
quenching counters, due to inefficiency inherent in the 
counters and due to apparent inefficiency arising from 
showers and scattering, have been measured. The ineffi- 
ciency inherent in the counters has been found to be 
almost entirely due to the dead time. The dead time, for 
the counters investigated and with the recording circuit 
employed, was 4X10-* second. This gives rise to an 
inefficiency of 0.2 percent in counters with a normal 
counting rate of 300 per minute. The apparent inefficiency 
of the counters and the error in coincidence measurements, 
for most of the experimental arrangements used were found 


INTRODUCTION 


S sources of error in coincidence measure- 
ments with a conventional cosmic-ray 
telescope we wish to consider the following 
effects: (1) chance coincidences, (2) inefficiency of 
the counters, (3). scattering in the absorber or in 
the counter walls, and (4) cosmic-ray showers. 
All four of these effects contribute to the “‘ap- 
parent inefficiency” of a counter; i.e., to the 
decrease in the counting rate when an additional 
counter is put into a cosmic-ray telescope. 

The existence of these sources of error has 
been recognized for many years. The inefficiency 
of a counter was determined with an anticoinci- 
dence method by Rossi in 1930,' and a counter 


1B. Rossi, Rend. Acc. dei Lincei [6] 11, 831 (1930). 


to be almost entirely due to showers. From the data 
presented below, we estimate that near sea level, with the 
counter telescope in the vertical direction and with 35 cm 
between the extreme counters, the coincidences due to 
side showers were about 15 percent of the normal coinci- 
dence rate when two counters were in coincidence, 7 
percent when three counters were in coincidence, and 3 or 
4 percent when the telescope contained five counters in 
coincidence. These percentages increased by a factor of 2 
when the apparatus was taken from 259- to 4300-m 
elevation. 


measurement of cosmic-ray showers was made 
by the same author in 1932.* In the years since 
then, there have been very numerous reports of 
measurements of these effects. However, in the 
determinations of the influence of cosmic-ray 
showers on telescope measurements, some effects 
have been included which have made the 
determination of the error inexact; and the 
measurements of the inefficiency of counters have 
been influenced by the recording of showers so 


2B. Rossi, Physik. Zeits. 33, 304 (1932). 

3 See, for example, T. H. Johnson, Rev. Mod. Phys. 10, 
193 (1938); Phys. Rev. 40, 638 (1932); D. K. Froman an 
J. C. Stearns, Can. J. Research A16, 29 (1938); M. E. Rose 
and W. E. Ramsey, Phys. Rev. 59, 616 (1941) ; M. Cosyns, 
Bull. Tech. de’l Assoc. Ing. sortis del Ecole Polytech. Brux. 
(1936) ; W. E. Danforth and W. E. Ramsey, Phys. Rev. 49, 
854 (1936). 
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that the results were only accurate for counters 
of rather low efficiency. 

The usual method of determining the ineffi- 
ciency of a counter has been to place three 
counters in line and to compare the threefold 
coincidence rate with the coincidence rate 
between the two extreme counters. However, 
part of the difference in the rates is due to the 
fact that side showers discharge the two extreme 
counters much more often than they discharge 
all three counters. If the efficiency is low, the 
effect of the showers is comparatively unimpor- 
tant. But for the highly efficient self-quenching 
counters, it is shown below that more than 
90 percent of the “apparent inefficiency”’ 
determined in this way is due to showers. 

The customary method of correcting for side 
showers has been to subtract the coincidence 
rate obtained with one counter placed just out 
of line from the rate obtained with all the 
counters in line. However, this correction 
includes not only the effect of side showers, 
but also the coincidences that occur when one 
particle traverses all the counters in line and an 
accompanying particle discharges the counter 
out of line. These cases should not be included 
in the correction, and give rise to a large error 
especially when electrons are being recorded. 

In the present experiments the method of 
measuring the inefficiency has been slightly 
refined so as to eliminate the effect of showers 
and allow a precise determination of the real 
inefficiency of the self-quenching counters. In 
addition, the apparent inefficiency has been 
measured by the usual methods under a wide 
variety of experimental conditions. By sub- 
tracting the real inefficiency and the effect of 
chance coincidences from the apparent ineff- 
ciency, the effect of side showers under these 
conditions has been found. 

This method of analysis gives accurately the 
apparent inefficiency due to side showers, but 
allows only an estimate of the number of side 
showers which produce coincidences in the 
counter telescope. The estimate is sufficiently 
accurate, however, to show that when the self- 
quenching counters are used, the coincidences 
due to side showers are a much more serious 
error under most experimental conditions than 
the error due to inefficiency. The precise determi- 


nation of the inefficiency makes it possible to 
remove a large part of this error by using a 
large number of counters in coincidence. Then 
the number of side showers recorded becomes 
extremely small and the correction may usually 
be neglected. The error due to inefficiency is 
thereby increased but is still very small and may 
be accurately known. 

The counters used in this investigation were 
made of brass cylinders 1 mm thick and 4.24 cm 
in diameter. The central wire was 4-mil tungsten, 
and contact was made to the wire by a thick 
Kovar rod (75-mils diameter) which emerged 
at one end of the counter through a glass seal. 
The length of the outer cylinder was 10 cm 
greater than the length of the central wire, and 
the length of the latter was 20 cm (except in a 
few cases specified below where the length was 
60 cm). The counters were filled with a mixture 
of 91 percent argon and 9 percent alcohol, to a 
total pressure of 11 cm of Hg. Each counter was 
operated at 90 volts above its starting potential, 
as determined by the recording circuit. The 
starting potentials were approximately 1000 
volts. 


I. INEFFICIENCY INHERENT IN THE 
SELF-QUENCHING COUNTERS 


A. Inefficiency Due to the Dead Time 


Part of the inefficiency of a counter is due to 
the fact that following each discharge there is a 


certain period of insensitivity, during which the - 


space charge which has accumulated around the 
central wire is removed from the immediate 
vicinity of the wire. Since this period depends 
on the motion of heavy ions, the order of 
magnitude is 10-* or 10-* second, which is much 
greater than the time required for an electron 
to traverse the counter; in fact, it is almost 
100 times as large as the time required for the 
discharge to be initiated, to rise to the maximum, 
and to be quenched. The dead-time effect does 
not depend on how far from the wire the first 
ions are created. 

During the course of this investigation, 
Stever! has published an analysis of the dead- 
time phenomenon in self-quenching counters, 
obtained by an oscillographic method. Stever’s 


4G. H. Stever, Phys. Rev. 61, 38 (1942). 
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HO cm 


Fic. 1. Experimental arrangement used in the measure- 
ment of the dead time. Counters 1,2,4,5,6 are in coinci- 
dence, counters 3 and S in anticoincidence. All of the 
anticoincidences are due to inefficiency of counter 3. Re- 
sults are in Table I. Results with counters 4 and 5 moved 
into line are given in Table II. When counters S are dis- 
connected, some of the anticoincidences are due to showers: 
see results in Table IX. 


photographs show that following a discharge 
there is a period of complete insensitivity (during 
which the counter acts as a proportional ampli- 
fier, and the pulses are too small to be seen on 
the oscilloscope), which is followed by a period 
of about the same length during which the pulse 
size gradually recovers to its full normal value. 
Obviously, regardless of the method used to 
detect the pulses, the period during which the 
counter is apparently ‘‘dead’’ will depend on the 
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sensitivity of the recording apparatus, just as 
the exact value of the “starting potential’’ of a 
counter also depends on the sensitivity of the 
detecting device. For our purposes, we wish to 
define the “‘dead time,’”’ ¢o, as the time interval 
following a discharge until the counter has 
recovered to the point where another pulse may 
be recorded by our circuit. This definition is 
made fairly precise by the fact that there is a 
small region of field strength during which the 
pulse size increases very rapidly with the field 
strength in a counter. Thus in the present 
experiments the counters were operated at 90 
volts above the recording circuit, and a decrease 
of 70 volts in the working voltage was found to 
produce no change in the measured dead time, 
within a statistical error of about 20 percent. 

In order to suppress any apparent inefficiency 
except that due to the dead time, the usual 
method of measuring inefficiency was modified 
as shown in Fig. 1. The counter for which the 
inefficiency was determined is in position 3, and 
was connected in anticoincidence. The four 
counters S were connected in parallel with 
counter 3 in order to eliminate anticoincidences 
due to side showers. Counters 4 and 5 were 
placed slightly out of line in opposite directions 
with displacements of } inch in the direction of 
the counter axes, and } inch at right angles to 
the axes. This displacement prevented the 
recording of anticoincidences due to any ineffi- 
ciency which occurs only near the edges of the 
counter, because a particle traveling very near 
the edge of counter 3 could not produce a 
coincidence. The lead was inserted between the 
counters so that practically all the coincidences 
recorded would be caused by mesotrons, which 
usually appear singly (at least at sea level), and 
which do not often undergo large-angle scatter- 
ing, such as occurs frequently in radiation 
processes of electrons. 

The criterion by which one can determine 
whether or not all the anticoincidences recorded 
with this arrangement are due to the dead time 
is the dependence of the apparent inefficiency on 
the discharge rate. If N is the discharge rate and 
to is the dead time, the inefficiency due to the 
dead time is (1—exp (— Ntéo)), or if Nto<1, it is 
simply Nt. Since the number of chance coinci- 
dences is completely negligible with five counters 
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in coincidence, and the other causes of apparent 
inefficiency do not depend on the discharge rate, 
the apparent inefficiency will be proportional to 
the discharge rate if, and only if, the anti- 
coincidences are all due to the dead-time effect. 
Accordingly, we have made the test with 
counters of different lengths in position 3, and 
during some of the measurements we have placed 
a radium source near the counters to stimulate 
the discharge rate. 

The results are shown in Table I. For each 
measurement, we have computed the dead time 
to with the formula 


Anticoincidences 
= Coincidences 


and thus have assumed that there is no ineffi- 
ciency except that which is due to the dead time. 
The consistency of the results obtained justifies 
the assumption.’ The weighted average of all 
the results is 


to=(3.8+0.2) X10~ sec. 


There is a small error in this determination of 
the dead time, because of the fact that occa- 
sionally when counter 3 fails to record the 
passage of a mesotron, a knock-on electron 
accompanying the mesotron discharges one of 
the side counters,S, and thus an anticoincidence 
is not registered. The magnitude of this effect 
was determined by repeating the above experi- 
ment with counter 3 disconnected. With this 
arrangement, we should record an anticoinci- 
dence whenever we record a coincidence, except 
for the cases when a knock-on electron discharges 
one of the side counters, and for the cases when 
a large shower from the side discharges all of the 
coincidence counters as well as one of the 
counters S. The ratio of the difference C—A 
between the numbers of coincidences and anti- 
coincidences to the number C of coincidences 
obtained in this measurement is 


(C—A)/C=(6.32+0.22) X10~. 


5 One measurement in the seven gave a result for fo which 
differed from the average by almost four times the sta- 
tistical error in the difference. If we assume the usual error 
curve for the probability of deviations, we find that the 
probability of obtaining a deviation this large or larger is 
1/5000. It is therefore unlikely, but not impossible, that the 
result is purely a statistical fluctuation. See footnote 6. 


Since any side shower which discharges all 
five coincidence counters is almost sure to 
discharge one of the counters S, the above ratio 
sets an upper limit for the error in the coincidence 
rate which is due to side showers discharging all 
of the counters in coincidence. Likewise it sets 
an upper limit for the error in t) due to the 
knock-on electrons discharging the side counters. 
The major part of the above effect is probably 
due to these knock-on electrons. If we assume 
that the effect is all due to the knock-on elec- 
trons, then the corrected value of to is 


to=(4.0+0.2) X10-* second. 


B. Spread of the Discharge 


One of the results in Table I is possibly in 
disagreement with the other results. The ques- 
tionable measurement was made with a long 
counter in position 3, so that a large part of the 
counter protruded outside of the counter tele- 
scope. The low value of the dead time obtained 


in this measurement led us to suppose that the. 


discharges which occurred in the “unused” part 
of the counter did not make the rest of the 
counter insensitive, or perhaps gave rise to a 
shorter dead time at the far end of the counter, 
i.e., that the discharge did not spread uniformly 
along the whole wire. 

In order to test this hypothesis, we repeated 
the measurement of the dead time under two 
conditions: (1) with the outer or ‘“‘unused’’ end 
of the counter stimulated by radium, and (2) 
with the inner or ‘‘useful’”’ end of the counter 
stimulated by radium. This was accomplished 
simply by shielding the end of the counter which 
was not to be stimulated with about 40 cm of 
lead. The result obtained was entirely negative ; 
the values found for the dead time were in 


TABLE I. Determination of the dead time to. The experi- 
mental arrangement is shown in Fig. 1. 


Coun- 

ter in Discharge Anticoinc. 

posi- Length rate {(1,2,4,5,6) to X104 

tion 3 (cm) (min.“) —(3+5S)] A/C X108 (sec.) 
a 20 299 46 1.84+0.27 3.7 
a 20 1090 64 8.0 +1.0 4.4 
b 20 251 27 1.87+0.36 4.5 
c 60 769 45 2.97+0.44 2.3 
c 60 2330 90 13.7 +14 3.5 
c 60 2406 92 16.0 +1.7 4.0 
d 60 794 58 5.53+0.73 4.2 
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Fic. 2. “‘Usefulness” of different sections of the counters 
as a function of distance from the wire. N(r)dr represents 
(in arbitrary units) the number of coincidence-producing 
rays which traverse a counter along paths with distances 
from the central wire lying between r and r+dr. The graphs 
are calculated for the counter arrangement shown in Fig. 1, 
with counters 4 and 5 moved into line. The assumption is 
made that over the small range of angle subtended by the 
diameters of the extreme counters the intensity of cosmic 
rays does not change. 


agreement with each other and with all of the 


_ other measurements of the dead time. The data 


have therefore been included in Table I, and are 
the two sets of data taken on counter c with high 
discharge rates. 

This result is also in agreement with Stever’s 
findings; i.e., that the discharge does spread 
uniformly along the whole counter wire. 


C. Inefficiency not Due to the Dead Time. 
Measurement of the Effective Diameter 


Rose and Ramsey* have shown the existence 
of a type of inefficiency which is not due to the 
dead time, and have investigated this inefficiency 
for counters filled with oxygen or with an 
argon-oxygen mixture. The effect which they 
have observed is greatest when the primary 
particle passes through the counter near the 
edge, and seems to depend critically on the type 
of gas used, being much greater for an oxygen- 
filled counter than for a counter containing an 
argon-oxygen mixture of which only 6 percent 
is oxygen. 

The causes for this phenomenon may be: 
(1) when the primary particle passes near the 
edge of a counter, the path length through the 
counter is short; therefore, it may happen 
occasionally that no ions are produced; or (2) 
when the particle passes near the edge, only a 


few ions are produced, and the field intensity is 
low; it is therefore likely that the electrons may 
be captured, and slowly-moving negative ions 
may be formed. This may either delay or prevent 
the discharge. The dependence of the inefficiency 
on the relative abundance of oxygen in the 
counters indicates that the formation of heavy 
ions is responsible for at least a part of this 
inefficiency. Because of the dependence on the 
number of ions produced, the probability for 
either of the processes (1) or (2) depends on the 
distance r from the wire to the path of the pri- 
mary ray in the form exp [—(a?—r’)!/k], where 
k is a constant and a is the radius of the counter. 
In addition, the probability for process (2) may 
have a term depending on the field strength, but 
this varies slowly with 7 except near the wire. 

It has been shown above that the inefficiency 
of self-quenching counters, except possibly near 
the edge, is very small (0.18 percent for the 
short counters) and can be accounted for entirely 
by the dead time. However, a function of the 
form given above may be negligible for small r 
and yet be very large for r~a. The net result 
would be to make the effective diameter of the 
counter slightly smaller than the geometrical 
diameter. 

The experiment described in Section A was 
repeated with counters 4 and 5 placed exactly 
in line. In this condition the ratio A/C of the 
anticoincidences to the coincidences gives a 
measure of the “‘over-all’’ inefficiency of counter 
3, which should be larger than the inefficiency 
measured in the first experiment if there is any 
additional inefficiency near the edges of the 


Taste II. Evidence for inefficiency near the edge of self- 
quenching counters. The edge effect is taken as the differ- 
ence between the “over-all” inefficiency (column 4) and the 
inefficiency due to the dead time (column 5). Weighted 
average of differences (discarding data on counter c):* 
Percent inefficiency due to edge effect =0.093+0.034 
percent. 


Anticoin- 
cidences A/C X10° A/C X108 
Coun- with with with 
ter in counters counters counters : 
Length 4 and 5 4 and 5 4 and 5 Differences 
3 (cm) _ inline in line out of line (A/C X10) 


20 49 2.66+0.38 1.84+0.27 0.82+0.47 
20 52 3.074043 1.87+0.36 1.20+0.56 
60 79 §.98+0.67 2.97+0.44 

60 87 6.10+0.65 5.53+0.73 0.57+0.98 


* See reference 6. 
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counter. The results are given in Table II and 
compared with the results of the first experiment.® 

The comparison in Table II indicates that the 
inefficiency of counter 3 is slightly greater when 
the coincidence-producing particles may traverse 
any part of the counter than when they cannot 
traverse the counter near the edge. The difference 
in efficiency, however, is very small, and is not 
quite definite, because the effect is only three 
times the statistical error. But this does not 
preclude the possibility that a small region near 
the edge of the counter may be very inefficient. 
In fact the effect of such an inefficiency on the 
“over-all” inefficiency of a counter depends on 
the position of the counter in the counter 
telescope, and is least for a counter near the 
center of the telescope. In Fig. 2 we have 
plotted the number of coincidence-producing 
rays traversing unit area of a counter, as a 
function of the distance r from the central wire 
to the path of the particle, for the counters in 
our telescope. It is seen that if the counter 
were completely insensitive within 1 mm of the 
outer walls, the resulting inefficiency would be 
only 0.28 percent for counter 3, but would be 
4.75 percent for one of the extreme counters. 
Thus the data in Table II may indicate that the 
outer mm of the counter has a low efficiency, 
which may be due to the failure of some particles 
to produce ions in the short length of path 
through the counter near the edge. This effect, 
if real, might be of considerable importance for 
the extreme counters in the array. 

In order to determine the effective diameter 
of the counters more exactly, the following 
experiment was performed. The counters were 
first arranged as in Fig. 3a, and by means of the 
coincidence-anticoincidence circuit we recorded 
simultaneously the coincidences (1,2,5,6) and 
the anticoincidences [(1,2,5,6)—(3+4+5) ]. 
Next, the counters were arranged as in Fig. 3b, 
and the same phenomena were recorded. In 


* Here again we have evidence that the inefficiency of 
counter c measured with counters 4 and 5 slightly out of 
line was subject either to large statistical fluctuation or to 
some other error, because the difference brought about by 
putting counters 4 and 5 in line is much larger for this 
counter than for the other three. If the inefficiency first 
measured were replaced by the expected result (calculated 
from the average dead time), the change brought about by 
ange counters 4 and 5 in line would agree with that found 
or the other counters. 


Fics. 3a AND 3b. Experimental arrangement used in 
measuring the effective diameter of a counter. Counters 
1,2,5,6 are in coincidence, counters 3,4, and S in anti- 
coincidence. 


arrangement (b), the anticoincidences were 
caused by: (1) particles traversing any inefficient 
region between counters 3 and 4, (2) inefficiency 
of counter 3 or counter 4 due to the dead time, 
or (3) shower particles discharging counters 1, 2, 
5, and 6 without discharging either of the 
counters 3, 4, or S. With arrangement (a), the 
effect of the inefficiency is negligible, so that the 
anticoincidences are due only to the last of the 
above-named causes. The shower effect should 
have been almost identical with arrangements 
(a) and (b), and was in any case small, because 
of the presence of the side counters S. The lead 
above counter 5 assured that we were recording 
mesotrons (which usually occur singly) rather 
than electron showers. 

The ratios of the anticoincidences to the 
coincidences obtained in this experiment were: 


A/C=8.2+0.30 percent with arrangement (b), 
A/C=0.4+0.10 percent with arrangement (a). 


The dead-time inefficiency of the counters 
accounts for an inefficiency of only 0.2 percent 
with arrangement (b), and zero with arrangement 
(a); therefore 7.6 percent of the coincidence- 
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producing rays traversed the inefficient area 
between counters 3 and 4 in arrangement (b). 
This figure is too low by about five percent (of 
itself) because of the mesotron secondaries 
generated above counters 3 and 4, which 
occasionally discharge these counters when a 
mesotron traverses the dead space ;’ hence the 
corrected result for the fraction of the particles 
which traverse the dead space is 8.0+0.4 percent. 

If we assume that this number is accounted 
for by the existence of a sharply defined and 
perfectly inactive dead space of width d between 
counters 3 and 4 in arrangement (b), we find 
(by reference to Fig. 2 and the graph for counter 
4) that d=1.94+0.10 mm. The thickness of the 
counter walls between the interior of counter 3 
and the interior of counter 4 was 2.0 mm. 
Therefore within 1/10 of a millimeter, the effective 
diameter of the counter is equal to the inside 
diameter of the counter wall. 

If the interior of the counters were perfectly 
efficient, however, the dead space found should 
have been less than the wall thickness, because 
the mesotrons traversing the counter walls 


" must often produce secondaries (which need not 


have very high energy) that discharge one of the 


Fic. 4. Experimental arrangement used to determine 
the effective length of a counter. Results are plotted in 


Fig. 5. 


7 See discussion and experimental results in’ Section IA. 


counters. Therefore these results are in agree- 
ment with the small difference found above 
between the “over-all” efficiency and the 
efficiency near the center of the counter. 

Finally, we may state that for self-quenching 
counters of 4-cm diameter and 20-cm length, 
filled with an argon-alcohol mixture, the effi- 
ciency near the center is 99.8 percent ;° but for 
rays traversing the counter within about } mm 
of the walls, the efficiency may be much smaller. 
As a result, the over-all inefficiency of the central 
counter in an array is about 0.3 percent. The 
edge effect in the extreme counters, however, is 
exactly compensated by the occasional discharge 
of the counters when a ray goes through the 
walls; therefore the effective inefficiency of these 
counters depends only on the dead time and is 
about 0.2 percent. 


D. Effective Length of the Counters 


The effective length of a counter was measured 
by using the arrangement shown in Fig. 4. 
The coincidence rate between the five counters 
was recorded as a function of the position of 
counter 4, as the latter counter was moved 
parallel to its own axis from a position completely 
out of the beam on one side to a position com- 
pletely out of the beam on the other side. 

The data are indicated by the open dots in 
Fig. 5, and the statistical errors are shown by 
the vertical lines. The solid curve is the theo- 
retical curve for a cylindrically-shaped counter 
of 4.2-cm diameter in our counter array, with 
an added background rate of 0.05 (min.~') 
which was obtained from the coincidence rate 
with the counter drawn completely out of the 
beam. The maximum height of the theoretical 
curve has also been taken from the experimental 
data, and the length has been chosen so as to 
give the best fit. 

The effective length of the counter, as obtained 
by this comparison between the theoretical and 
experimental curves, is /=18.8 cm. The actual 
length of the fine collecting wire in this counter 
was 20 cm. Since the outer cylinder extended 
5 cm beyond this wire at both ends, and since 
contact was made to the wire by means of a 

8 This value, of course, as explained above, is subject to 


variations depending upon the normal counting rate and 
upon the sensitivity of the recording circuit. 
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Fic. 5. Determination of the effective length. The open dots represent the 
coincidence rates obtained with the arrangement shown in y 4 for various 
or 


positions of the middle counter. The solid curve is calculated 


counter, with background rate added. 


thick rod, it had been expected that the effective 
length would be equal to the length of the fine 
central wire; however, it appears that the 
effective length is less than this by 1.2 cm. Asa 
result of this measurement, we find that the 
absolute intensities of cosmic rays reported 
recently by one of us® are all too small by 
10.3 percent. 


Il. APPARENT INEFFICIENCY DUE TO 
SHOWERS AND SCATTERING 


A. Apparent Inefficiency Due to Scattering 


With the apparatus shown in Fig. 1, the ratio 
of anticoincidences to coincidences recorded 
was increased by a factor of 5 when counters 4 
and 5 (those placed slightly out of line) were 
disconnected from the coincidence group. That 
is, the anticoincidences [(1,2,6)—(3+S)] were 
almost 1 percent of the coincidences (1,2,6) 
whereas the anticoincidences [1,2,4,5,6—(3+S) ] 
were only 0.2 percent of the coincidences 
(1,2,4,5,6). This result contrasts strongly with 
the small increase in anticoincidence rate, due 
to the edge effect, obtained when counters 4 and 
5 were not disconnected but placed exactly in 
line with the others. The result might easily 
have been explained had the scattering of 
mesotrons in the lead absorber been appreciable, 


°K. Greisen, Phys. Rev. 61, 212 (1942). 


a cylindrical 


since a fraction of the mesotrons traversing 
counter 6 might have been scattered by the lead 
so as to traverse counters 1 and 2 but miss 
counter 3. 

In order to test whether or not scattering 
was really responsible for the anticoincidences 
[(1,2,6)—(3+5S)] observed, the latter were 
recorded both as described above and also with 
the width of lead absorber reduced to exactly 
the width of the counters. In the latter condition, 
scattering could not possibly cause counters 1, 
2, and 6 to be discharged without counter 3 also 
being discharged. No difference was observed 
between the two measurements (within a 
statistical error of 12 percent of the effect). 
Later measurements, described below, showed 
that the entire increase in the apparent ineffi- 
ciency of counter 3 upon disconnecting counters 
4 and 5 was due to showers, in spite of the 
presence of the side counters S. 

The unimportance of the scattering effect in 
measurements on mesotrons at sea level is also 
predicted theoretically. Using the formula given 
by Rossi and Greisen,'® one finds that the root 
mean square angle of scattering for mesotrons 
of momentum 2X10° ev/c in 7 cm of lead is 
only 1.5 degrees. 


 B. Rossi and K. Greisen, Rev. Mod. Phys. 13, 240 
(1941); formula given on page 265. 
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Figure 7 


\\ 


Figure 8 


Fics. 6-9. Experimental arrangements used to study the 
effect of side showers on coincidence measurements. Re- 
sults corresponding to Fig. 6 are in Table III; results 
corresponding to Fig. 7, with varying numbers of counters 
in coincidence and counter 2 in anticoincidence, are in 
Tables IV, V, VIII, X, XI; results corresponding to Fig. 8 
ah a VI; and results corresponding to Fig. 9 are in 

able 


B. Apparent Inefficiency Due to Showers, and 
the Effect of Side Showers on Coincidence 
Measurements 


The effect of showers on an array of counters 
is strongly dependent on the number of counters 
in the array, their geometrical arrangement, the 
physical surroundings of the apparatus, and the 
altitude at which measurements are made, as 
well as on the type of absorber used and on its 
shape and position. In the experiments described 
below these factors have been varied in an 
attempt to determine the effect of showers on 
coincidence measurements made under various 
conditions. The experiments are grouped below 
according to the number of counters connected 
in coincidence. 


1. Experiments with Two Counters 
in Coincidence 
(a) Results with experimental arrangement indi- 
cated in Fig. 6.—The results are given in Table 
III. In order to distinguish the effect of showers, 
the number of chance coincidences had to be 
determined. Therefore measurements were made 


AND N. NERESON 


both with no artificial stimulation and with the 
counters stimulated by radium. The difference 
between the coincidence rates observed was 
assumed to be due to chance coincidences. The 
resolving time thus determined is r=7.7+0.6 
microseconds. With this value of 7, we have 
calculated the chance coincidence rates given in 
Table III, and by subtraction of the chance 
coincidences from the total we have obtained 
the shower rate. From these figures we see that 
when the counters are not stimulated with 
radium, 96 percent of the coincidences are due 
to showers and only 4 percent are chance 
coincidences. 

From the data in Table III we may also 
discover how many of the showers which 
discharge counters 1 and 2 also discharge one or 
more of the five anticoincidence counters 
surrounding counter 2. An _ anticoincidence 
(1,2-X) may occur when a chance coincidence 
(1,2) occurs, provided none of the counters X 
is also discharged. However, the anticoincidence 
rate with no artificial stimulation was 0.072 
(min.—'), while the chance coincidence rate was 
only 0.018. Therefore at least three-fourths of 
the anticoincidences are not due to chance 
coincidences. These anticoincidences are ac- 
counted for by showers discharging counters 1 
and 2 without discharging any of the counters X. 
Thus the anticoincidence rate due to showers is 
at least 0.054 (min.~'); i.e., at least 14 percent 
of the showers (0.40 per minute) which discharge 
counters 1 and 2 fail to discharge any of the 
five surrounding counters X. 

(6) Results with experimental arrangement indi- 
cated in Fig. 7—This experiment was performed 
outdoors, under a thin wooden roof at elevations 


TaBLE III. Measurements taken with the experimental 
arrangement shown in Fig. 6, in the basement of the physics 
building in Ithaca (elevation 259 m). 


Counters 
No artificial stimulated 
Quantity measured stimulation by radium 
Coinc. (1,2) per min. 0.417+0.015 0.908 +0.033 
Anticoinc. (1,2-X) per 0.072+0.006 0.470+0.024 
min. 
Counts per min. in coun- 276 1586 
ter 1, N; 
Counts per min. in coun- 257 1254 
ter 2, Nz 
Chance coinc. per min. 0.018+0.0013 0.51 +0.038 
Shower rate (min.~) 0.40 +0.015 0.40 +0.050 
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TaBLE IV. Measurements taken in the vertical direction in Ithaca with the apparatus shown in Fig. 7. Counters 3, 4, 
and 5 were not in use. Ni, N2, and Ng are the counting rates in the individual counters 1, 2, and 6. 


No artificial stimulation Counters stimulated by radium 
No lead 13-cm lead No lead 13-cm lead 
between between between between 
Quantity measured counters counters counters counters 
Coinc. (1,6) per min. 3.14 +0.022 2.46 +0.042 3.32 +0.048 2.60 +0.037 
Anticoinc. (1,6-2) per min. 0.197 +0.006 0.119+0.009 0.313+0.015 0.239+0.011 
Nz (min.~) 252 204 806 736 
N, (min.~) 252 220 515 475 
Ng (min.~) 268 218 1127 1086 
Anticoinc. rate due to inefficiency 0.005 0.003 0.018 0.013 
Chance coinc. per min. 0.015+0.002 0.011+0.001 0.126+0.013 0.113+0.012 
Anticoinc. per min. due to side showers 0.177 +0.006 0.105+0.010 0.169+0.020 0.113+0.016 


TABLE V. Measurements taken at various elevations and two zenith angles, with and without lead, with the apparatus 
shown in Fig. 7. Counters 3, 4, and 5 were not in use. C refers to the coincidence rate between counters 1 and 6; A is 
the anticoincidence rate [1,6-2]. The anticoincidence rates may be considered as proportional to the number of side 


showers which discharge counters 1 and 6. 


Elevation Absorber Zenith angle 0° Zenith angle 46° 

(meters) (cm of lead) A (min.~) A/C (percent) A (min.~) A/C (percent) 

259 0 0.197 +0.006 6.30.2 0.139+0.008 10.0+0.6 

13 0.119+0.009 4.9+0.4 0.110+0.009 9.5+0.8 

1616 0 0.35 +0.02 8.0+0.5 0.23 +0.02 11.9+1.2 

13 0.20 +0.02 6.1+0.5 0.23 +0.02 13.6+1.5 

3240 0 0.81 +0.03 10.6+0.5 0.52 +0.03 15.4+0.9 

13 0.39 +0.02 8.4+0.5 0.48 +0.03 19.0+1.1 

4300 0 1.31 +0.04 11.8+0.4 1.02 +0.04 20.5+0.8 

13 0.78 +0.03 12.5+0.6 0.77 +0.03 22.2+1.1 


TABLE VI. Measurements taken with the experimental 
arrangement shown in Fig. 8, in the basement of the 
physics building in Ithaca (elevation 259 m). 


Counters 
No artificial stimulated 
Quantity measured stimulation by radium 
Coinc. (1,2,3) permin. 0.125 +0.006 0.139+0.008 
Ratio of anticoinc. 9.9 +1.6 13.3 +2.2 
(1,2,3-X) to coinc. 
(1,2,3) in percent 
Coinc. (1,2) per min., 38 38 
12 
Counts per min. in 276 1586 
counter 3, N; 
Chance coinc. (1,2,3)  0.0027+0.0002 0.015+0.0011 
per min. 
Coinc. (1,2,3) permin. 0.122 +0.006 0.124+0.008 
due to showers 


259, 1616, 3240, and 4300 meters, both with and 
without lead absorbers between the counters. 
The thickness of the absorber when in use was 
13 cm of lead, which is sufficient to prevent 
electrons from traversing the counter telescope. 
The absorber had the same width as the counters, 
so that it could not possibly scatter particles in 
such a way as to produce a coincidence without 
discharging counter 2. 

An analysis of anticoincidence rates (1,6—2) 


obtained in Ithaca with the counter telescope in 
the vertical direction is given in Table IV. The 
anticoincidences are due to (1) inefficiency of 
counter 2, (2) chance coincidences between 
counters 1 and 6, and (3) side showers which 
discharge counters 1 and 6 without discharging 
counter 2. The effect due to inefficiency has been 
calculated by using the value of the dead time 
determined above (4X 10~‘ second). 

The increase in the anticoincidence rates 
which accompanied the stimulation with radium 
allows two independent determinations (from 
the measurements with and without lead) of 
the resolving time. The two values thus found 
are 6+1 and 7+1 microseconds. The average 
of these two values has been used to determine 
the total number of anticoincidences arising 
from chance coincidences. 

We see from the values in Table IV that 
inefficiency normally accounts for less than 3 
percent of the anticoincidences, chance coinci- 
dences account for about 8 percent, and side 
showers account for about 90 percent of the 
anticoincidences. 

The measurements at other zenith angles and 
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TABLE VII. Measurements taken with the experimental arrangement shown in Fig. 9, in the basement of the physics 
building in Ithaca (elevation 259 m). 


Counter a (20-cm length) 
in position X, no 
artificial stimulation 


Quantity measured 


Counter d in position X 
(60-cm length), no 


Counter a in position X, 
artificial stimulation 


stimulated with radium 


Coinc. (1,2,3) per min. 2.47 +0.018 2.55 +0.042 2.53 +0.032 
Anticoinc. [(1,2,3) —(X+.S)] per min. 0.026+0.0019 0.041 +0.0054 0.032 +0.0036 
Counts per min. in counter X, N, 299 1090 794 
Chance coinc. (1,2,3) per min. 0.003 0.011 0.003 
Anticoinc. per min. due to inefficiency 0.005 0.018 0.013 
Anticoinc. per min. due to side showers 0.018+0.002 0.012 +0.006 0.016+0.004 


other altitudes are presented in Table V. Here 
we have not separated the anticoincidences into 
their various components. However, over the 
range of altitudes included in our measurements 
the chance coincidences between two counters 
vary with altitude in very nearly the same way 
as do the showers; therefore the showers account 
for about 90 percent of the anticoincidences at 
all of the altitudes and zenith angles. 

Since an appreciable fraction of the side 
showers which discharge the two coincidence 
counters will also discharge the anticoincidence 
counter and hence fail to produce an anticoinci- 
dence, the number of side showers which produce 
coincidences must be even larger than these 
anticoincidence rates. Thus we see that the 
error due to side showers, which would be 
involved in a coincidence measurement with two 
counters, is a large error, completely over- 
shadowing the errors due to inefficiency and 
chance coincidences. The error becomes more 
serious at greater zenith angles of the counter 
telescope, where the coincidence rates decrease 
faster than the shower rate. And the error 
becomes more serious at high altitudes, because 
the shower rate increases with altitude faster 
than does the mesotron intensity. In this 
experiment the shower rate was 7 times as large 
at Mount Evans (elevation 4300 m) as at Ithaca, 
and the ratio to the intensity of single particles 
was twice as great at Mount Evans as at Ithaca. 


2. Experiments with Three Counters 
in Coincidence 
(a) Result with experimental arrangement shown 
in Fig. 8.—The results are given in Table VI. 


The chance coincidence rates have been calcu- 
lated by using the value of the resolving time 


determined above; the remainder of the coinci- 
dences are assumed to be due to showers. © 

In order to have a check on the reliability of 
the separation of the anticoincidence rate into 
its components, the experiment was performed 
both with no artificial stimulation and with the 
counting rates stimulated by radium. The 
results are seen to be in good agreement, since 
the shower rate should have been the same in 
the two conditions. 

The figures in Table VI indicate that when 
the counters are not stimulated by radium, 98 
percent of the coincidences are due to showers. 
The anticoincidence rate shows that 10 percent 
of the showers fail to discharge any of the five 
surrounding counters. 

(6b) Results with experimental arrangement 
shown in Fig. 9.—This experiment differs from 
the one described just above only in that 
counter 3 is placed in line with counters 1 and 2. 
This is the usual arrangement of three counters 
for intensity measurements. For our present 
investigation we are primarily interested in the 
anticoincidences [(1,2,3,) —(X+S) ]. These are 
due to (1) inefficiency of counter X, (2) chance 
coincidences between counter 3 and counters 1 
and 2, and (3) side showers. The chance coinci- 
dences and the effect of inefficiency have been 
calculated by using the values of the resolving 
time and the dead time given above, and the 
remainder of the anticoincidences is due to the 
side showers which fail to discharge any of the 
counters X or S. That the effect of scattering is 
negligible has already been shown (see Section 
IIA). 

The results are given in Table VII. In order 
to have a check on the determinations, some of 
the measurements were made with a long counter 
in position X (counter d, of length 60 cm) and 
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some with the counters stimulated by radium; 
for these two cases the inefficiency is larger, but 
the shower rates should be the same as when 
counter a is used (20-cm length) with no artificial 
stimulation. The results verify this prediction. 

If we compare these results with the coinci- 
dence rate itself, we see that for counter a in 
position X, the number of chance coincidences 
between the three counters is only 0.1 percent 
of the coincidence rate, and the anticoincidences 
due to inefficiency are only 0.2 percent of the 
coincidence rate, whereas the number of anti- 
coincidences due to showers is 0.7 percent of the 
coincidence rate. But the first experiment with 
three counters (Fig. 8, Table VI) showed that 
90 percent of the showers which produce a 
threefold coincidence also discharge one of the 
counters X or S and so do not produce anti- 
coincidences. Hence the number of coincidences 
produced by side showers is 10 times the number 
of anticoincidences produced by the showers, or 
7 percent of the coincidence rate. This indicates 
that even in the basement of a building at low 
altitude, the error due to side showers in coinci- 
dence measurements with three counters is very 
considerable. 

The experiment described just above (Fig. 9, 
Table VII) was repeated with the side counters 
disconnected and only counter X in anticoinci- 
dence. The ratio of the anticoincidences to the 
coincidences obtained in this condition was 
A/C=(3.3+0.6) percent. This is larger than 
the ratio obtained from the data in Table VII 
(i.e., 1.06 percent, for counter @ in position X), 
because when the side counters are disconnected, 
a larger fraction of the showers which produce 
coincidences also produce anticoincidences. 
Chance coincidences and inefficiency are re- 
sponsible for only 75 of the 3.3 percent. Since 
the shower rate, however, is 7 percent of the 
coincidence rate, we see that about half of the 
showers which discharge the three counters in 
coincidence also discharge the single counter X 
and thus fail to produce an anticoincidence. 
This fact is used in the analysis of the following 
experiments. It also indicates the way to elimi- 
nate the errors due to showers. If we had used 
four counters in coincidence instead of three, 
the fourth being placed in the position of the 
anticoincidence counter, the number of coinci- 


dences produced by showers would have been 
reduced by a factor of 2. 

(c) Results with experimental arrangement 
shown in Fig. 7.—The results are given in 


Table VIII. We know from the experiments 


described above that the anticoincidences arising 
from inefficiency and from chance coincidences 
are only a small fraction of the total number of 
anticoincidences, especially when no side coun- 
ters are used.'' Therefore the anticoincidence 
rates given in Table VIII may be considered as 
proportional to the number of showers dis- 
charging the coincidence counters. Moreover, 
since about half of the showers which discharge 
the three coincidence counters also discharge the 
anticoincidence counter (as shown above), we 
may consider the total number of showers which 
discharge the three counters in coincidence as 
about twice the number of anticoincidences. I? 
is seen that this implies a large error in coinci- 

TABLE VIII. Measurements taken at various elevations, 
outdoors, with the apparatus shown in Fig. 7. Counters 3 
and 4 were not in use and no absorber was used. C refers to 
the coincidence rate between counters 1, 5, and 6. A is the 
anticoincidence rate [1,5,6-2]. The anticoincidences are 


almost entirely due to showers, and are about half the 
number of showers which produce coincidences. 


Eleva- Zenith angle 0° Zenith angle 46° 
tion A/C Cc 
(meters) A (min.~!) (percent) A (min.~!) (percent) 


259 0.105+0.008 3.6+0.29 0.042+0.005 3.2+0.39 
1616 0.153+0.014 3.8+0.34 0.094+0.017 5.4+1.0 
3240 0.45 +0.035 6.340.51 0.161+0.021 5.9+0.80 
4300 0.69 +0.044 6940.46 0.29 +0.032 7.6+0.88 


TABLE IX. Evidence for side showers discharging an 
array of five counters in line. The data were taken with the 
apparatus shown in Fig. 1, in the basement of the physics 
building in Ithaca. The figures in the first three rows of the 
table represent anticoincidence rates in counts per 100 
minutes. 


Counters 4 and § 
displaced as Counters 4 and 5 
shown in Fig. 1 placed in line 


A X 100, with only counter 0.76+0.10 1.13+0.17 
3 in anticoinc. 

AxX100, with counter 3 0.40+0.06 0.64+0.09 
plus side counters S in 
anticoinc. 

Difference 0.36+0.12 0.49+0.19 

Percent ratio of difference 0.17+0.06 0.20+0.08 


to coinc. rate 


1 This has only been shown for the measurements in 
Ithaca, but we | be even more true at high altitudes, be- 
cause the shower rate is known to increase more rapidly 
with altitude than the total counting rates. 
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TaBLeEs X and XI. The following data were taken with the apparatus shown in Fig. 7, outdoors, under a thin wooden 
roof. The anticoincidences are mainly due to the side showers which discharge the coincidence array but miss the anti- 


coincidence counter. 


TABLE X. Anticoincidences (1,3,4,5,6 -2) per 100 min. 


Zenith angle 


Elevation Absorber in 
(meters) cm of lead o° 29° 46° 56° 
259 0 3.4+0.25 2.4 +0.27 1.6 +0.23 1.04+0.17 
13 2.2+0.20 0.94+0.18 0.79+0.17 0.62+0.15 
1616 0 6.2+0.27 3.6 +0.56 2.5 +0.42 1.5 +0.34 
13 2.8+0.51 2.1 +0.43 0.98+0.25 1.1 +0.26 
3240 0 14.6+1.1 7.6 +0.74 4.5 +0.65 4.0 +0.49 
13 5.6+0.72 44 +0.59 2.9 +0.52 1.9 +0.36 
4300 0 21.6+1.4 13.6 +1.1 10.2 +1.0 5.6 +0.71 
13 8.5+0.87 6.2 +0.74 5.2 +0.70 4.6 +0.65 
TABLE XI. Percent ratio of anticoincidences (1,3,4,5,6 —2) to coincidences (1,3,4,5,6). 
Elevation Absorber in Zenith angle 
(meters) cm of lead o° 29° 46° 56° 
259 0 1.19+0.09 1.14+0.13 1.28+0.18 1.40+0.24 
13 0.94+0.09 0.56+0.11 0.73+0.16 0.94+0.22 
° 1616 0 1.56+0.07 1.27+0.20 1.44+0.25 1.45+0.33 
13 1.01+0.18 0.93+0.19 0.76+0.20 1.28+0.31 
3240 0 2.19+0.17 1.65+0.16 1.73+0.25 2.66+0.34 
13 1.34+0.17 1.49+0.20 1.45+0.26 1.62+0.30 
4300 0 2.30+0.15 2.01+0.17 2.62+0.25 2.5340.33 
13 1.60+0.17 1.48+0.18 2.00+0.27 2.84+0.41 


dence measurements with 3 counters, the 
magnitude of which becomes more serious at 
high altitudes. 


3. Experiments with Five Counters 
in Coincidence 


The experiments with three counters in 
coincidence have shown that when one of the 
counters was surrounded by five counters in 
anticoincidence, 90 percent of the side showers 
which discharged the three counters also dis- 
charged one of the anticoincidence counters, 
and hence failed to produce an anticoincidence. 
When the side counters were not used but only 
one counter in line was connected in anticoinci- 
dence, about half of the side showers recorded 
discharged the single anticoincidence counter. 
If five counters are used in coincidence instead 
of three, the side showers which discharge the 
array will be much fewer and, in general, will 
contain more particles. Such a shower would 
have a larger probability of discharging the 
anticoincidence counters. Therefore when five 
counters are used in coincidence, the set of five 
anticoincidence counters surrounding one coinci- 
dence counter should be discharged by practi- 


cally all of the showers recorded, and we should 
obtain no anticoincidences due to showers. This 
conclusion is verified in the first experiment 
reported in this paper, where all of the anti- 
coincidences could be accounted for by the 
dead-time inefficiency. Moreover, if the side 
counters are disconnected and only one counter 
(in line with the coincidence array) is used in 
anticoincidence, more than half of the side 
showers recorded will discharge the single 
anticoincidence counter. Thus the number of 
side showers which discharge the five coincidence 
counters is several times the number of anti- 
coincidences recorded due to the showers. 

(a) Results with experimental arrangement 
shown in Fig. 1.—The results are given in 
Table IX. The difference between the two 
anticoincidence rates should represent the num- 
ber of showers which discharge the five counters 
in coincidence and fail to discharge counter 3. 

It is seen that the number of anticoincidences 
due to showers is about as large as the number 
due to the dead-time inefficiency ; the number of 
fivefold coincidences due to showers must be 
several times as large. 

(b) Results with experimental arrangement 


Aon 
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shown in Fig. 7.—The results are given in 
Tables X and XI. The anticoincidences occurring 
at Ithaca because of the dead-time inefficiency 
should be about 0.2 percent of the coincidence 
rates. There should be no anticoincidences due 
to scattering because the lead absorber had the 
same width as the counters. The anticoincidences 
arising from chance coincidences between the 
coincidence counters are quite negligible. The 
effect of inefficiency near the edge of the anti- 
coincidence counters might amount to 0.1 percent 
of the coincidence rate. The remainder of the 
anticoincidences should be due to showers 
discharging the anticoincidence counter. Thus 
we see that even at 259-m altitude the major 
part of the effect (about 70 percent) is due to 
the side showers. 

At higher elevations the effect of inefficiency 
near the edge of the anticoincidence counter 
should remain the same fraction of the coinci- 
dence rates as at Ithaca. The inefficiency due 
to the dead time should go up in the same ratio 
as the counting rates. But the effect of showers 
increases with altitude much faster than does 
the counting rate. Therefore at the higher 
elevations the fraction of the anticoincidences 
which is due to showers is even larger than at 
Ithaca. The anticoincidence rates may thus be 
taken as a proportional measure of the number 
of showers which discharge the counter array. 
However, the number of side showers which 
discharge the five counters in coincidence is 
several times as large as the number which 
produce anticoincidences. 

The variation of the anticoincidence rates with 
altitude, shown in Table X, should give a 
measure of the variation of the shower intensity 
with altitude, if the anticoincidences are caused 
mainly by showers. We may check this by 
computing the average ratio of the anticoinci- 
dence rates at the higher altitudes to the anti- 
coincidence rates at Ithaca, and comparing these 
ratios with the corresponding ratios of the 
intensities of the soft component, measured 
recently at the same altitudes by one of us.® 
This is done in Table XII. We have also given 
in Table XII similar ratios for the total intensity 


TABLE XII. Evidence that anticoincidences are mostly 
due to the soft component of cosmic rays. The second and 
third columns refer to data taken with the experimental 
arrangement shown in Fig. 7 (see Tables V, X, XI). The 
last two columns give the corresponding ratios for the in- 
tensity of the soft component and for the total intensity. 


Average Average 
ratio of ratio of 
anticoinc. anticoinc. 
rates, withS rates,with2 Ratio of 
counters intensities 
(1,6) of the soft Ratio of total 
in’ ties 


(meters) in coinc. incoinc. component 
259 1.0 1.0 1.0 1.0 
1616 1.6 1.8 1.9 1.4 
3240 3.5 3.9 4.4 ae 
4300 6.1 6.9 8.2 3.2 


of all particles which discharge a counter 
telescope of six counters (taken from the same 
reference). The comparison verifies the fact that 
the major part of the anticoincidences, both 
when five counters and when two counters were 
in coincidence, was due to some agent (showers) 
which increases in intensity with altitude in the 
same way as the electron intensity. 

An interesting fact, noted in Table XI, is that 
not only the anticoincidence rates but also the 
ratio of anticoincidences to coincidences de- 
creases when the lead absorber is put in place, 
by an average factor of 1.4. This may be at- 
tributed both to the absorption of shower 
particles and to the multiplication of the shower 
particles in the lead. That is, the lead prevents 
a shower of only a few low energy particles 
from producing a coincidence. Moreover, when 
a shower of high energy particles, or of a large 
number of particles, strikes the apparatus, the 
multiplication in the lead makes it very unlikely 
that the anticoincidence counter will not be 
discharged. Thus fewer showers produce coinci- 
dences, and a smaller fraction of these produce 
anticoincidences, when lead is used. The distribu- 
tion of the 13 cm of lead in layers 2 cm thick 
between the counters assists to make this 
multiplication effect pronounced, because 2 cm 
of lead is nearly the optimum thickness for 
shower production. 

The authors wish to express their thanks to 
Professor Bruno Rossi for suggesting these 
experiments and for his guidance in the analysis 
of the results. 


e anti- 
2 
3 
1 
) 
3 
l 

ould 
This 
nent 
anti- 

the 

side 

nter 
d in 
side 
ngle | 
r of 
once 
nti- 
nent 
in 

two 
um- 

ters 
1Ces 
ber 

r of 
be 
rent 


OCTOBER 1 AND 15, 1942 


PHYSICAL REVIEW 


VOLUME 62 


' A Comparison of X-Ray Line Breadth and Internal Friction for Alpha-Brass as 
Affected by Cold-Working and Annealing 
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X-ray line breadth measurements are obtained for samples of alpha-brass cold-stretched 
in differing amounts and then annealed at various temperatures. A comparison of these data 
with internal friction measurements made on the same samples reveals marked dissimilarities. 
For small amounts of cold-work, line breadths appear to be affected by the rate of cold-working, 
whereas internal friction seems more dependent upon the amount of cold-work. The internal 
friction introduced by cold-work is removed by low temperature annealing (far below the 
temperature of recrystallization). Line breadth changes little with annealing until just before 
recrystallization when it exhibits a comparatively sharp decrease. 


INTRODUCTION 


HE effects on x-ray line breadth and 
internal friction introduced by cold-work- 
ing metals have both been subjects of consider- 
able study. But little work has been done in 
relating these effects for the same samples of 
metal. One writer,! measuring internal friction 
at large strain amplitudes, has given some 
qualitative comparisons of line breadths and 
internal friction measurements as affected by 
cold-working and annealing. It is of interest to 
obtain a more definite quantitative analysis of 
these effects, especially for internal friction 
measurements at small strain amplitudes—less 
than those producing plastic flow. The effects of 
cold-working and annealing on the internal 
friction of alpha-brass have been studied by 
Zener’ and his collaborators. It is the purpose of 
this paper to compare their results with x-ray 
line breadth measurements obtained by the 
writers for specimens from the same samples of 
brass given heat treatment identical with that 
used for the internal friction measurements. 
Two possible causes of line broadening have 
been discussed at considerable length in the 
literature. These are broadening due to lattice 
distortion, and broadening due to fragmentation 
of the grains. 
A survey of the experimental work which 
attempts to distinguish between these effects in 
cold-worked metals has led the writers and 


1 John T. Norton, Trans. A. I. M. E. 137, 49 (1940). 
2C. Zener, H. Clarke, and C. S. Smith, Metals Tech. 
(Oct. 1941); (A. I. M. E. Tech. Pub. No. 1376). 


others** to adopt the view that the primary 
cause of line broadening in cold-worked metals 
is distortion of the crystal lattice. Reasons for 
this view may be summarized briefly as follows: 

1. There is a lower limiting grain size (indi- 
cated by maximum line breadth) which may be 
produced by cold-working a metal. In most cases 
the calculated broadening from grains of this 
size is negligible.5-7 

2. Two distortional effects which are inter- 
preted as warping of the reflecting planes* or 
changes in the shape of the unit cell®" have 
been demonstrated exclusive of any possible 
fragmentation effect. 

3. Certain metals, viz., aluminum and zinc, 
exhibit no line broadening when severely cold- 
worked at ordinary temperatures though con- 
siderable fragmentation must occur.” ” 

4. Line breadths, when broadening is intro- 
duced by cold-working the metal, are found not 
to be dependent on the wave-length of radiation 
used as those arising from small grain size are 
known to be.** 

5. Comparisons of intensity measurements 


3G. W. Brindley, Proc. Phys. Soc. 52, 117 (1940). 

4]. E. Wilson and L. Thomassen, Trans. A. S. M. E 
22, 769 (1934). 

5 K. Becker, Zeits. f. Physik 42, 226 (1927). 

6 F. E. Haworth, Phys. Rev. 52, 613 (1937). 

7™W. A. Wood, Proc. Roy. Soc. A172, 231 (1939). 

8 E. Orowan and K. J. Roscoe, Nature 148, 467 (1941). 

8W. A. Wood, Nature 131, 842 (1933). 

10 W. A. Wood, Phil. Mag. 19, 219 (1935). 

"1 V. Caglioti and G. Sachs, Zeits. f. Physik 74, 647 
(1932). 

2 U. Dehlinger, Zeits. f. Krist. 65, 615 (1927). 

3G. W. Brindley and P. Ridley, Proc. Phys. Soc. 50, 
501 (1939). 
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with line breadths indicate that it is impossible 
to reconcile the large primary extinction effect 
with a sufficiently small grain size to produce 
broadening.* 


METHOD 


The x-ray lines referred to in this paper are 
the familiar Debye-Scherrer diffraction circles 
obtained from a polycrystalline substance by the 
back-reflection method. The radiation used was 
the Ka doublet from the copper target of a 
metal Coolidge tube operated at 32,000 peak 
volts with full wave rectification. Exposures 
were of the order of 250 milliampere-hours. 

The line breadths were measured from traces 
of a microphotometer having a magnification of 
about 12.5 times. The measure adopted was the 
so-called full width at half maximum. This was 
taken to be the width of the a; line measured at 
the maximum of the a. Such a measure is 
independent of the sensitivity of the micro- 
photometer. All measurements were made on 
the second-order reflection from the (210) planes. 


MATERIAL 


The samples of alpha-brass used in these 
measurements were prepared and furnished by 
Dr. C. S. Smith of the American Brass Company. 
A complete description of their preparation is 
given elsewhere.’ Briefly, it is as follows: Bars 
of 70-30 alpha-brass were cold-drawn to size, 
annealed for two hours at 400°C, and cooled in 
air giving a grain size of 0.0012 cm. These bars 
were then cold-stretched in the amounts shown 
in Table I (reproduced in part from (2)). Small 
specimens were cut from these bars, faced in a 
lathe, ground flat with carborundum powder and 
etched to a depth of one or two millimeters to 
remove the surface effects of cold-working 
introduced in the preparation of the specimens. 
Care was taken to avoid excessive heating during 
the cutting and etching processes. The etching, 
done in a 30 percent solution of concentrated 
nitric acid at 15°C, was found to be quite 
important. A few samples exhibited a slight 
decrease in line breadth when given a second 
etch, so that in the preparation of subsequent 
specimens the amount of the etch was increased 
considerably. Finally all measurements were 
repeated for a second set of specimens prepared 


from the same lot and consistent results were 
obtained. 

The annealing process was nearly identical 
with that which had been used on these samples 
in making the internal friction measurements. 
Anneals (except that at room temperature) 
were for one hour at various temperatures from 
room temperature up to 375°C, which tempera- 
ture is well above that of recrystallization for all 
samples. The rate of cooling after annealing was 
found to be relatively unimportant. Samples 
cooled several minutes in air gave exactly the 
same results as those cooled several hours in an 
oil bath. 


RESULTS AND DISCUSSION 
(a) Results of the X-Ray Measurements 


Samples 5 and 6 were found to give results 
almost identical with those obtained from the 
unstrained control sample, number 2. None of 
these gave evidence of line broadening, indicating 
that there was little lattice distortion. Numbers 
7 and 8 produced broader and more continuous 
circles but with the doublet still well enough 
resolved to enable accurate measurement in the 
manner described above. Number 9 gave very 
broad and diffuse lines which could not be 
measured accurately in this way, but the effects 


TABLE I. The cold-working of the samples. Samples 5, 
6, and 7 were cold-stretched at a constant rate of 0.05 in. 
per minute; samples 8 and 9, at a rate 12.5 times as rapid. 


Maximum strain Residual per- Internal 
Bar No. underload % manent set % friction 
2 unstrained 0.000 1/Q=0.35 x 10° 
5 0.286 0.100 0.95 
6 0.700 0.475 2.30 
7 2.76 2.46 5.40 
8 2.78 2.47 4.90 
9 30.5 30.0 4.3 


of annealing on this sample could be observed 
qualitatively. 

Because of the negligible broadening observed 
for samples 5 and 6, only 2, 7, 8, and 9 were 
annealed. Sample 2 differed from the others only 
in that it had not been cold-stretched. Thus the 
unchanging appearance of this sample during 


4 These anneals do not properly belong on the graph in 
Fig. 1 since they were for a period of several months. 
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ANNEALING TEMPERATURE IN DEGREES CENTIGRADE 


Fic. 1. X-ray line breadth as a function of annealing temperature for samples 
7 and 8 cold-stretched to about 2.5 percent elongation. Number 8 stretched 12.5 


times as rapidly as number 7. 


annealing indicates that the decrease in line 
breadth observed with the annealing of the other 
samples is to be interpreted as recovery from 
lattice distortion. 

The results of annealing on samples 7 and 8 
are shown in Fig. 1. The line breadth for both 
exhibited first a slight increase as the tempera- 
ture of anneal was increased, then a sharp 
decrease just before evidence of recrystallization 
appeared. Number 7 had already been annealed 
at 225°C for a preliminary measurement. This 
explains the absence of data for lower tempera- 
ture anneals of this sample. The line breadths 
for sample 9 were observed to change in a similar 
way with annealing but because of the very 
diffuse appearance of these lines and the greatly 
increased amount of background scattering, 
accurate measurements could not be obtained. 

Evidence of recrystallization (if the appearance 
of sharp spots on the Debye-Scherrer circles 
may be taken as an indication of recrystalliza- 
tion) was first observed after an anneal at 350°C 
for sample 8 and one at 375°C for number 7. 
Figure 1 shows that for both samples this 


occurred after a large amount of the broadening 
due to cold-working had been removed by 
annealing. This decrease in line breadth before 
evidence of recrystallization occurs is thought 
by some writers to substantiate the view that 
distortion of crystal lattice is the primary cause 
of line broadening. 

The difference in the positions of these curves 
is also significant. Numbers 7 and 8 were both 
cold-stretched 2.5 percent but number 8 was 
strained 12.5 times as rapidly as was number 7. 
The data in Fig. 1 indicate that number 8 
exhibited greater line broadening and recovery 
at a lower annealing temperature. Both these 
effects are known to be associated with a greater 
amount of cold-work. Thus line breadth meas- 
urements appear to be sensitive to the rate at 
which the metal is deformed—at least, for small 
amounts of cold-work. 


(b) Comparison with Internal Friction Data 


The internal friction as measured in terms of 
the rate of damping of very small amplitude and 
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high frequency vibrations was found to increase 
rather rapidly with the amount of cold-work. 
Table I shows that for sample 6, which produced 
negligible line broadening, the internal friction 
was more than 40 percent of the maximum value 
observed which was that for number 7 unan- 
nealed. However, it was found that samples 7 
and 8, cold-stretched the same amount but at 
considerably different rates, gave very nearly 
the same internal friction measurements. This 
suggests the conclusion that for small amounts 
of cold-work the x-ray line breadths are the more 
sensitive to the rate of deformation whereas the 
internal friction is to a greater extent dependent 
upon the amount of cold-work. Since this result 
was not anticipated the evidence for it is some- 
what meager, and it must remain for the present 
a suggestion which the writers make with 
considerable reserve but one which should 
certainly be the subject of further study. 

Still another dissimilarity is seen in the fact 
that with increasing amounts of cold-work the 
internal friction first increases to a maximum 
and then begins to decrease whereas x-ray lines 
are found to increase more slowly to a maximum 
and to be unaffected by further cold-working. 
Also these maxima occur for very different 
amounts of cold-work. In this case a strain of 
about 2.5 percent produces a maximum internal 
friction and one of 30 percent actually produces 
less. The x-ray line breadths are much greater 
at 30 percent than at 2.5 percent deformation 
and are known to reach a maximum for this 
order of strain. It appears that there is no 
correlation between internal friction and x-ray 


line breadth measurements as affected by 
cold-working. 

With successively higher temperatures of 
anneal the internal friction was found to decrease 
and reach a minimum at 225°C for all samples. 
This minimum value is approximately equal to 
that of the unstrained sample 2. The internal 
friction introduced by cold-work thus has been 
removed by annealing at a temperature about 
100° below that of recrystallization. But as 
may be seen from Fig. 1 there is no corresponding 
change in line breadths at this temperature.® 

Finally, the internal friction is found to 
increase rapidly with higher temperature anneals 
above the recrystallization temperature. This is 
attributed to grain growth, and the internal 
friction is interpreted in terms of inter-crystalline 
thermal currents.'*!” There is no corresponding 
change possible in x-ray line breadths. These 
observations indicate that there is also a marked 
lack of correlation between internal friction and 
line breadth measurements as functions of the 
annealing temperature. 
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18 Preliminary work on this observation is referred to in 
C. Zener, Phys. Rev. 60, 455 (1941). 

16 R. H. Randall, F. C. Rose, and C. Zener, Phys. Rev. 
56, 343 (1939). 

17C, Zener and R. H. Randall, Trans. A. I. M. E. 137, 
41 (1940). 
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Intensities of continuous x-radiations from a thin target (199A) of metallic nickel under 
electron bombardment in the range 12 to 180 kv have been observed. Ross filters of Ag-Pd 
(485-508 x.u.) and Se-As (979-1044 x.u.) provided monochromatization, and measurements 
were performed with an argon-filled ion chamber and a quadrant electrometer. Sensitivities 
of the measuring system at various wave-lengths were determined and corrections were made 
for all absorptions in the path of the radiation. Small corrections for radiation from sources 
other than the target and for the several effects of the finite target thickness were applied. 
The corrected data show relative intensities J per unit wave-length interval for the selected 
wave-lengths and direction of observation as produced by electron bombardment of inde- 
pendent nickel atoms. The observed variation of J with bombardment energy is in excellent 
agreement with the Sommerfeld theory and the variation of J with wave-length is in approxi- 


mate agreement. 


INTRODUCTION 


S a field for quantitative measurement the 
continuous x-ray spectrum has been long 
neglected. Well known for several decades have 
been the qualitative aspects of the spectrum 
which proceeds from thick targets, including the 
wave-length-intensity distribution and its de- 
pendence upon bombardment energy, target 
atomic number, and (less certainly) direction of 
emission. We have also possessed one quantita- 
tive fact: the law of Duane and Hunt. Although 
more refined and quantitative thick-target data 
might have a certain technical interest, pure 
physicists have not been attracted to this field 
of measurement because the thick target is so 
special and complex a source as to defy exact 
theoretical analysis. Lacking the possibility of 
comparing measurements with theory, observers 
have regarded the field as unremunerative. 

The continuous x-radiation from thin, or vir- 
tually monatomic, targets under electron bom- 
bardment is susceptible of theoretical treatment,! 
being, indeed, the result of one of the most 
fundamental of possible collisions. There is no 
doubt here as to the utility of quantitative in- 
formation but experimental difficulties have stood 
in the way, and we are at this writing very far 
from the ability to describe the radiation from 
such sources in quantitative terms. Kulen- 


1A. Sommerfeld, Ann. d. Physik 11, 257 (1931); F. 
Sauter, Ann. d. Physik 18, 486 (1933); G. Elwert, Ann. d. 
(1942) 34, 178 (1939); R. Weinstock, Phys. Rev. 61, 584 


kampff,? Nicholas,* and others have observed 
continuous spectra from targets less thick than 
the productive’ layer of thick targets but not 
sufficiently thin to classify as atomic targets. 
These observations have been valuable in show- 
ing the conspicuous differences between thick- 
and thin-target spectra and in indicating qualita- 
tive agreement with thin-target theory. In none 
of these cases was an attempt made to deduce 
the precise shape of an atomic-target spectrum, 
nor would this have been easy, with the available 
data. To reduce thin-target spectra to atomic- 
target spectra, in the present state of our knowl- 
edge of the behavior of fast electrons in solids, 
it is necessary to have the target thin enough or 
the electrons fast enough so that the velocities 
of the latter (both speeds and directions) are 
modified by the target atoms only to that small 
extent for which theoretical corrections may be 
applied. 

To meet the conditions of theory the experi- 
menter must produce sufficiently thin targets 
and he must use observing means sensitive 
enough to measure the feeble radiation after 
spectral resolution. Targets of evaporated and 
condensed nickel satisfied the first condition in 
the present work. Because of their great ad- 
vantage in light-gathering power we have used 
Ross! filters for isolating the radiation bands to 


2H. Kulenkampff, Ann. d. Physik 87, 579 (1928). 

3 W. W. Nicholas, Bur. Stand. J. Research 2, 837 (1929). 

4P. A. Ross, Phys. Rev. 28, 425 (1926) ; J. Opt. Soc. Am. 
and Rev. Sci. Inst. 16, 433 (1928). 
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Fic. 1. Plan view of thin-target x-ray tube and measuring a 
mately to scale. J. we reading current to target. 2. 
Nickel-on-cellulose acetate target. 4. Clock-con- 
alanced ao and iron prefilter. 6. Argon-filled ion chamber. 
7. Electrometer. 8. Anode plate, eeted 
e filament. /0. Sylphon to permit aiming of 


cylinder target. 
trolled shutter. 5 


containing the target. 9. Cath 


tus approxi- 
ollow aluminum 


to pass electrons to field-free region 


electron stream. 11. Cathode Socatiion adjustment controlling shield at 9. Over- 


all length of x-ray tube = 130 cm. 


be studied. With these filters the radiation was 
intense enough for measurement by standardized 
ionization methods. Nickel was selected as the 
first target element in a projected series of 
studies covering a wide range of atomic numbers. 
This program as a whole would include intensity 
and polarization measurements with a variety of 
atomic numbers, electron energies, wave-lengths, 
and angles of observation. 


APPARATUS 


The nickel target was produced by a modifica- 
tion of the evaporation technique described by 
Pockman and Webster.’ A backing film of 
cellulose acetate 940A in thickness was mounted 
in vacuum 10 cm above a nickel ribbon which 
could be heated, to its melting point if desired, 
by a well-controlled alternating current of about 
100 amperes. It was possible to keep track of 
the progress of the deposit by measuring through 
photoelectric means the progressive darkening 


5L. T. Pockman and D. L. Webster, Rev. Sci. Inst. 12, 
389 (1941). 


of a glass window of the evaporation chamber. 
By previous calibration the constant relation 
between the window opacity and the surface 
density of the metallic deposit on the acetate 
film was known. All data to follow were obtained 
from a target whose calculated thickness was 
199A. 

During the evaporation the target was sup- 
ported on a light and flexible structure of 
aluminum wires as described by Pockman and 
Webster.’ The completed target on its supports 
was transferred for use to the x-ray tube shown 
in Fig. 1 where it is seen to be mounted across 
the axis of a deep aluminum cup. The cup 
served to catch that great majority of the bom- 
barding electrons which the target failed to 
stop; it was made of a substance of low atomic 
number to facilitate retention of the electrons 
and it was insulated, except for the microam- 
meter lead, because it was essential to measure 
the current to or through the target. 

It was an important feature of the design in 
Fig. 1 that the fragile target was inclosed in a 


E 62 
_| 
all] | 
| 
whe 
1 | 
rved Ef 
not 
zets. 
1OW- 
\ick- 
lita- 
1one 
luce 
um, 
able 
nic- 
lids, 
h or 
ities 
are 
nall 
be 
eri- 
rets 
tive 
fter 
and 
in 
ad- 
sed 
to 


336 


field-free metal box, protecting it at all times 
from the kind of destructive gas bursts which 
were the ruin of so many targets during 
earlier researches in this laboratory. One wall 
of this protective enclosure was penetrated by a 
4.8-mm hole through which electrons of the 
desired energy were injected, to fall almost 
normally upon the target surface. The cathode 
which opposed this hole at a distance of about 
10 cm (subject to adjustment) was equipped 
with a precisely controllable focusing shield and 
aiming mechanism which made it possible to 
direct a large fraction of the total emission 
through the hole. Those electrons which struck 
metal around the hole developed heat (a maxi- 
mum of 50 watts) which was removed by a 
current of oil flowing under pressure within the 
hollow wall. 

. Two pairs of Ross filters were used and data 
were thereby obtained concerning the intensities 
of two spectral bands, a wide variety of electron 
energies being used in each case. Filter elements 
and pass band wave-lengths were Ag-Pd (485-— 
508 x.u.), and Se-As (979-1044 x.u.) As the 
preparation of these and other filters is to be 
described in a separate publication, little need be 
recorded here except that the Ag-Pd filters were 
metal sheets and the Se-As pair were thin plates 
of the powdered metals in paraffin, all filters 
being practically of optimum thickness. The two 
members of each pair were balanced for equal 
transmission of radiation outside the pass band 
by rolling or pressing. The precision of balance 
was such that the difference between the in- 
tensities transmitted by two balanced filters at 
selected wave-lengths lying without the pass 
band was about } percent of the incident in- 
tensity at those wave-lengths. 

In most applications of balanced filters the 
effective thickness of a filter may be varied by 
tilting it so that the radiation passes through 
obliquely. In our case the size of the focal spot 
was such that the beam incident upon the filters 
was appreciably non-parallel. To tilt a filter 
transmitting such a beam is to vary the thickness 
in different ratios for the different incident rays, 
an operation almost completely incompatible 
with accurate balancing. For this reason the 
filters used in this work were brought to balance 


K. HARWORTH AND P. KIRKPATRICK 


with their planes parallel and were always used 
in this way. 

In using these filters an iron prefilter of known 
transmission was placed in the beam to minimize 
any effect caused by possible departure from 
perfect balance in the spectral region of the 
nickel K lines. 

The ionization chamber was designed to permit 
comparable intensity measurements at different 
wave-lengths. Specifically, the collectors were flat 
and parallel aluminum plates separated suffi- 
ciently to prevent interference with photoelec- 
trons, and provided with guard plates which in- 
sured precision in the location of the limits of 
the collecting region, 15.2 cm in length. Large 
Cellophane windows at both ends facilitated 
geometrical alignment, while contributing en- 
tirely negligible amounts of secondary radiation 
to the collecting region, 8.61 cm distant. The 
chamber was filled with dried tank argon to a 
pressure 18 mm below the existing atmospheric 
pressure, and mounted with its entrance window 
98 cm from the center of the target focal spot. 
The direction from target to ion chamber made 
an angle of 93.5° with the direction of bombard- 
ment. Ionization charges were measured with a 
quadrant electrometer operating at a moderate 
and continually checked charge sensitivity. 

The filters were mounted just before the en- 
trance window of the chamber where a hole in a 
lead plate defined the beam area in such a way 
that the calibrated regions of the filters were 
effective and that the beam traversing the 
chamber lay midway between the plates. A care- 
ful survey along the beam assured that the entire 
focal spot could irradiate the chamber window 
and that a minimum of other matter could con- 
tribute either primary or secondary x-radiation. 
A test of this latter condition will be described 
below. 


ADJUSTMENTS AND MEASUREMENTS 


By varying the adjustment of the focusing 
shield of the cathode and applying a small bias 
voltage between the shield and the filament the 
diameter of the focal spot was kept at 1 cm for 
all tube potentials used in this research, a range 
from 8 to 180 kv. This constancy of focal spot 
size was perhaps not essential but it avoided the 
extremes of a spot so large as to include the 
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RELATIVE SENSITIVITY 


WAVE-LENGTH ANGSTROMS 
02 04 06 08 10 12 16 


Fic. 2. Sensitivity curve of measuring system. Ordinates 
are proportional to collectible ionization charges resultin 
from emission toward the chamber of equal quantities 
radiation energy. 


target supports or other thick material and one 
so small as to produce, at the allowable current 
density, insufficient radiation for good measure- 
ments. As it was, the target was bombarded with 
a current of some 200 microamperes for the 
hundred hours or more required by these meas- 
urements without any observable reduction of 
the ample radiation output. The radiation meas- 
ured in this investigation proceeded from the 
target in a slightly diverging beam whose axis 
made an angle of 93.5° with the forward direction 
of motion of the bombarding electrons. 

In a typical run for the collection of intensity 
data the tube was operated on the pumps until 
a pressure of at most 10-* mm Hg could be 
maintained with an applied voltage of 180 kv 
and the normal target current flowing. At volt- 
ages above 80 kv a 5-megohm resistor was con- 
nected in series with the tube at the anode end 
to reduce the violence of the occasional gas 
burst. The precise value of this resistor was 
checked at the beginning and end of every run 
and its voltage drop deducted from readings of 
the high potential voltmeter to obtain the actual 
bombardment voltage. While one operator de- 
voted his attention exclusively to precise main- 
tenance of the selected tube voltage and current 
another observed the target current and operated 
the filters and electrometer, recording the ioniza- 
tion charges accumulated in clock-controlled 
exposure intervals with one and then the other 
of two balanced filters in the beam. 

Readings were not affected by electrometer 
drift since the electrometer was kept grounded 
and detached from the ion chamber until the 


exposure was concluded, at which time the 
chamber collector was connected to the elec- 
trometer, producing a ballistic deflection. As only 
a fraction of the tube current passed by way of 
the target, a special microammeter was mounted 
in a shielded case at the potential of the anode 
compartment so as to carry only the target 
current, and the mean reading of this instrument 
for every exposure was recorded. 

After five or six readings with each filter the 
tube voltage was lowered and observations re- 
peated at the new setting, continuing until the 
excitation voltage of the filter band was reached. 
There resulted an isochromat curve for the filter 
band wave-length, but a curve in need of the 
corrective treatment described below in order to 
bring out its fundamental significance. The 
curves to follow are the combined result of some 
fifteen runs of this kind. 


INTERPRETATION OF DATA 


Though no intensity measurements in absolute 
units were attempted it was desired to deduce 
from the data the isochromats of the free atom 
with all ordinates on the same arbitrary intensity 
scale. This required correcting the data for the 
unequal response of the ion chamber to different 
wave-lengths, for the presence of radiation from 
sources other than the target, and for the several 
effects of the finite bulk of the thin target upon 
speed and direction of the bombarding electrons, 
effects themselves speed-dependent. 

In an ion chamber with guard plates the 
radiation must undergo absorption in a pre- 
liminary gas volume before arriving at the col- 
lecting region. The preliminary absorption neces- 
sarily makes the chamber insensitive to radiations 
of long wave-length. On the other hand the 
finite length of the collecting region makes for a 
low sensitivity to short wave radiations. The 
result is a rather symmetrical sensitivity curve 
like Fig. 2, which represents the result of calcula- 
tions on the chamber of this investigation, calcu- 
lations taking account not only of the above 
considerations but also of the small effects of 
scattering and fluorescence within the chamber, 
and of the reduction of intensity by absorptions 
external to the chamber. We are involved with 
scattering and fluorescence effects because the 
fraction of the entering radiant energy which is 
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scattered out of the collecting region and fails 
therefore to contribute any ions to the collected 
charge is a function of wave-length, and because 
the escaping argon K radiation energy intro- 
duces a loss dependent upon wave-length. The 
method of making these corrections is essentially 
that of Allison and Andrew® with simplifications 
which our particular gas, wave-length range, and 
chamber design permitted. 

The absorbers outside the chamber included 
the target itself, thin windows of the tube and 
chamber, intervening air, prefilter, and balanced 
filters. By taking due account of the transmission 
ratios of these elements we arrive at the over-all 
sensitivity curve of Fig. 2, any two of whose 
ordinates stand in the ratio of the ionization 
charges which would be collected at the chamber 
as a result of primary emission toward the 
chamber of equal energies of radiation having, 
respectively, the wave-lengths denoted by the 
corresponding abscissas. It is improbable that in 


RELATIVE INTENSITIES 


60 80 100 
BOMBARDMENT ENERGY 


120 140 160 
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Fic. 3. Relative intensities of monochromatic x-rays from 
independent nickel atoms under electron bombardment. 
Ordinates are proportional to continuous spectrum radiant 
—_—~ per unit solid angle per unit wave-length interval 
or the selected direction of observation. Solid curves are 
experimental results, with probable errors of observation 
denoted by vertical line segments. Broken curves are 
theoretical (Sommerfeld-Weinstock). Theory and experi- 
ment have been brought into agreement at one point by 
arbitrary manipulation of the ordinate scale. 


Fig. 2 the ratio of any ordinate to the maximum 
ordinate is in error by more than 0.01. To put 
observations on the same intensity scale the 
value of each collected charge was divided by the 
appropriate ordinate of Fig. 2. 

Allowance for radiation from other sources 
than the target had to be deducted from the in- 
tensity data. The target was removed and the 


*S. K. Allison and V. J. Andrew, Phys. Rev. 38, 441 
(1931). 
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tube operated as before with the result that with 
all filters removed small but observable intensi- 
ties were produced at all voltages. The intensities 
being insufficient for monochromatic measure- 
ments, corrections were calculated by assuming 
for the radiation received at each voltage a 
Wagner and Kulenkampff’ thick-target spectral 
distribution and normalizing the total to agree 
with the observations. The background intensity 
contributions lying within the Ross filter pass 
bands were then readily deduced and applied to 
the thin-target intensity data as subtractive cor- 
rections, the magnitudes of which ranged from 
0.1 to 3 percent of the thin-target intensities. 
Although our target was so thin that nearly 
all the bombarding electrons emerged with devi- 
ations of only a few degrees due to multiple 
scattering, and with energy losses of the order of 
one percent there were nevertheless four separable 
effects of target thickness which affected the 
continuous spectrum intensities. Of these. the 
velocity loss has been mentioned. Its magnitudes 
for various bombardment energies were calcu- 
lated by Williams’s* empirical stopping power 
formula, effective bombardment voltages for the 
target as a whole being calculated and used in 
plotting the final results. The other three target 
effects are processes which cause the bombarding 
electron’s actual path in the target to exceed the 
normal target thickness. This may result from 
multiple scattering, from single scattering, or 
from return of the electrons into the target after 
large angle scattering by the cellulose acetate 
backing. The average actual path lengths in the 
target were calculated by a combination of the 
Rutherford® scattering formula and Bothe’s” 
formula for the most probable angle of scattering. 
At voltages above 50 kv the calculated path 
lengths exceeded the target thickness by only 
about one percent but at 15 kv the increased 
multiple scattering had raised the excess to about 
8 percent. The details of these calculations are 
elaborate and laborious and are omitted here 
because of the smallness of the resulting correc- 
tions and because they will be more extensively 


( 923) Wagner and H. Kulenkampff, Physik. Zeits. 23, 503 
1 
8 E. J. Williams, Proc. Roy. Soc. A130, 310 (1930). 

®E. Rutherford, Phil. Mag. 21, 669 (1911). 

10W. Bothe, Handbuch der Physik, ed. 1927, Vol. 24, 
15-18. 
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treated in a forthcoming paper to be submitted 
to this journal by Webster, Pockman, Harworth, 
and Kirkpatrick. 

After all data were obtained the target was 
examined and found in good condition. It was 
quite flat and exhibited a well-defined, darkened 
focal spot 1 cm in diameter. Under the micro- 
scope numerous minute holes could be seen, 
occupying a total of perhaps 2 percent of the 
focal area. Aside from this observation no at- 
tempt was made to evaluate the thickness dis- 
tribution over the focal spot. It is realized that 
if thin spots as well as perforations existed our 
target corrections are somewhat too large since 
they were based upon an assumed thickness of 
199A, the value obtained from the known loss 
of nickel from the filament in the evaporating 
process. We were watchful at all times, however, 
for decreases in the radiation output of the target 
and since none was observed, it is concluded that 
only negligible amounts of metal were lost during 
operation. 

RESULTS 


In Fig. 3 the corrected isochromats are plotted 
as solid lines with the statistical probable errors 
of the ordinates denoted by short vertical seg- 
ments. The errors of the abscissas are negligible. 
These curves are in striking contrast to the nearly 
straight lines of positive slope which are the 
isochromats of a thick target. Recently through 
the work of Weinstock! the Sommerfeld theory 
of the continuous x-ray spectrum has been stated 
in a form which makes comparison with experi- 
ment possible with only a moderate amount of 
straightforward calculation and in Fig. 3 the 
predictions of this theory are shown as broken 
lines, which have been normalized by arbitrary 
adjustment of the ordinate scale to produce a 
compromise fit of theory to experiment in the 
case of the 0.497A curves. In comparing the 
shapes of the curves we consider that the agree- 
ment is all that might have been hoped for in 
view of the approximations of the theory (neglect 
of screening and relativity) and the errors of the 
experiment. As to the dependence of intensity 
upon wave-length (the vertical separation of the 
isochromats) theory indicates a more pronounced 
wave-length variation than we observed and the 
discrepancy here is beyond the expectations and, 
for the present, is left unexplained. 


Seeking empirical formulas for representation 
of the experimental results one finds that the 
isochromats are approximately equilateral hyper- 
bolas, especially the 0.497A curve, and that all 
of the observations may be approximated by 
IV\!=constant where J is the observed intensity 
within a fixed small wave-length interval, V is 
the bombardment energy, and A the wave- 
length. The inaccuracies of fit exceed the errors 
of observation only near the high frequency 


Fic. 4. Relative intensities of monochromatic x-ra 
from independent nickel atoms, plotted to show the o 
served inverse proportionality of intensity to electron 
bombardment energy. The validity of the approximation is 
poy higher when frequency and bombardment energy 
are high. 


limits, where the intensities fall off somewhat. 
This simple descriptive formula may hardly hope 
for theoretical justification and it is advanced 
only tentatively and subject to adjustment in 
view of future results. Figure 4 presents the 
experimental data of Fig. 3 in such a way as to 
display this empirical inverse proportionality of 
intensity to bombardment energy. 

Naturally the theoretical isochromats may be 
represented by a similar abbreviated expression 
but in this case the exponent of A lies between } 
and 2. It is interesting that twenty-five years 
ago Webster" was able to deduce from thick- 
target data the conclusion that the intensity in 
the thin-target spectrum should vary inversely 
with \2, a conclusion here approximately sub- 
stantiated by both experiment and an adequate 
theory. 


1D. L. Webster, Phys. Rev. 9, 220 (1917). 
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An investigation of the range or absorption characteristics of secondary electrons in mag- 
nesium by means of examining (with high energy primary electrons) the secondary electron 
emission ratio 6 as a function of the thickness x; of magnesium layers evaporated on a carbon 
backing has led to an expression (/;,=/o:.*) for the ranges of secondary electrons as a function 
of the energy E of the secondary electrons at the point of production in the solid. An analysis 
of the results yields relations between 6 and x; (for the primary energies used) as well as the 
separate contributions of backing and layer to the total secondary emission. Values of /o1. 
obtained from this experiment on magnesium for two, four, six, and eight kilovolts are 1.3, 
1.7, 2.2, 2.3X10-* cm/(volt)?, respectively. The depth of origin of secondary electrons in 
the energy range from 10 volts to 200 volts will vary from 2X 10-7 cm to 4X 107° cm. 


INTRODUCTION 


N investigation of the range or absorption 
characteristics of secondary electrons in 
magnesium has led to an expression relating the 
ranges of secondary electrons to the energy of 
those electrons at the point where they were 
produced in the magnesium. The range-energy 
relation was obtained by examining the second- 
ary electron emission ratio as a function of the 
thickness of magnesium layers evaporated on a 
carbon backing. 

The method to be described for determining 
the range of the secondary electrons requires 
primary electrons of high energy. The energy 
of the primary electrons must be large enough to 
make their penetration sufficiently great so 
that the number of secondary electrons that 
they produce will be a slowly varying function of 
the primary electron energy. 

Magnesium was chosen for this experiment 
because it was hoped that comparisons could be 
made between the results obtained for magne- 
sium and the results obtained for magnesium 
evaporated through an atmosphere of oxygen 
because the latter is well known for its high 
secondary yield while the metal has a low yield. 
This comparison, however, was not actually 
made. 


APPARATUS AND TECHNIQUE 


The vacuum system (Fig. 1) contained the 
electron gun and collector for secondary elec- 


* At present with RCA Laboratories at Princeton, New 
Jersey. 


trons, and a target assembly which held seven 
polished carbon blocks upon which were evap- 
orated the magnesium layers of different thick- 
nesses. The assembly could be rotated by means 
of a ground glass joint to expose any of the seven 
surfaces to electron bombardment from the gun. 
An evaporation furnace was used for some of this 
work but was finally removed because a more 
satisfactory method was found for controlling 
and measuring layer thickness. The evaporation 
mechanism used in place of the furnace is shown 
in Fig. 2. This evaporation was done in a bell 
jar, and the samples were then transferred to the 
tube of Fig. 1 for the secondary emission meas- 
urements. The time required for transferring the 
samples was about an hour. Exposing the surfaces 
to air for this period of time produced no ap- 
preciable change in the secondary emission 
characteristics (with the primary energies used). 
This fact is demonstrated in the curves at Figs. 
5, 6, 7, and 8 because three of the points in each 
figure were taken on surfaces evaporated from 
the furnace while the others were taken from 
surfaces prepared in the bell jar arrangement 
and transferred to the tube for measurements. 
The pressure in the bell jar during the evapora- 
tion of these surfaces was less than 10-° mm of 
mercury. The surfaces used were evaporated 
quite slowly, over a period of twenty to thirty 
minutes. Surfaces produced by rapid evaporation 
do not have the density characteristic of metallic 
magnesium. Above a certain evaporation rate 
the density of the surface varies with the rate of 
evaporation. 
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The method of determining the thickness 
of the layers of magnesium was as follows: 
Three or four microscope cover glasses (care- 
fully dried and weighed) were placed at various 
points on the ladder arrangement shown in Fig. 2. 
The amount of material evaporated on each 
cover glass was determined by weighing before 
and after evaporation. If one assumes the inverse 
square law relation and the source of evaporation 
to be a point, the layer thickness ¢ is 


t=— cos 
r2 


K was evaluated from the amount of material 
found on each cover glass, and from the geometry 
of the evaporating arrangement cos@ and r 
were determined (i.e., measurement of x and d, 
Fig. 2). Values of K found by this method did not 
deviate from the mean of the three values for 
three cover glasses by more than about 12 


Fic. 1. Vacuum system with electron gun, target assembly, 
and evaporating furnace. 


1 J. Strong, Procedures in Experimental Physics (Prentice- 
Hall, 1938). 
? P. Bonet-Mourey, Ann. de Physique 11, 279 (1929). 


evaporated 
material 
on carbon blocks 


Fic. 2. Arrangement used for evaporating the surfaces in a 
bell jar. 


percent. An average value of 1.8X10-? cm* was 
used for K. 

By suitable design and arrangement of aper- 
tures the number of stray secondary electrons 
from these apertures was made negligible. One 
aperture was of carbon and so shaped as to 
reduce as far as possible stray secondaries so 
that the secondary electrons measured all came 
from the target. 

The values of 6 were determined by dividing 
the secondary electron current 7, by the primary 
electron current i,; 7, and 1, were measured by 
means of microammeters in the circuit shown in 
Fig. 3. Measurements were made for primary 
beam energies of 2, 4, 6, and 8 kilovolts. Plots of 
5 as a function of magnesium layer thickness 
x, are given in Figs. 4, 5, 6, and 7 for the primary 
energies mentioned. The values of 6 for x; =0 and 
different primary energies are those obtained 
for the carbon which formed the backing for 
the magnesium layers. (See Fig. 12 for’~ energy 
curve for carbon.) 

The data were taken for a range of layer thick- 
nesses from x,=0 to x;=2.5X10- cm. Addi- 
tional values of 6 were found for the four primary 
energies for x;=2.2X10- cm; this layer was 
sufficiently thick so that few of the primary 
electrons succeeded in traversing it. The value 
of 6 for x; =2.2 X10-* cm is shown on each curve 
by the short lines on the upper left and upper 
right. 

The unlabeled curve in each case is a plot of 
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c 


Ms Ma Ba 


= POWER PACK 0-800 VOLTS 
B2=POWER PACK 0-10000 VOLTS 

M, = MILLIAMMETER 

M2 M3Mg = MICROAMME TERS 

C = SECONDARY-ELECTRON COLLECTOR 
Sc=SCREEN AND CARBON APERTURE 

P, = FIRST ANODE 

T = TARGET ASSEMBLY 


Fic. 3. Circuit used for measurements. 


the data. The other two curves on each figure 
are calculated from the data by a method dis- 
cussed later. Of these two one is the calculated 
contribution of the carbon backing material to 
the secondary emission ratio 6, and the other is 
the calculated contribution of the magnesium 
layer to the ratio 6. The sum of these contribu- 
tions gives the experimental curve. 


ANALYSIS OF RESULTS 


An analysis of the problem has been made in 
the following manner : Consider that in the region 
of the magnesium layer (see Fig. 8) where 
0=x,;=x the number of secondary electrons 
formed in the energy range dE in the distance dx 
and going in the negative x direction is 


dny,= —Mipkifi(E)dEdx 


(in what follows the subscript 2 refers to backing). 
m1, is the number of primary electrons per 
square cm of path reaching dx. 

k, is a coefficient for the production of second- 
ary electrons (headed out of the layer in the 
negative x direction) by primary electrons. k; isa 
function of primary energy, but the primary 
energy has been made high in order to make k; 
almost independent of the primary energy. 

fi(E) is the energy distribution function giving 


the number of secondary electrons produced in 
dx in the energy range between E and E+dE. 
fi(E) is associated with the magnesium layer. 
There will be another distribution f2(Z) of this 
kind associated with the backing. f:(E) and 
fe(E) should be sufficiently alike so that no 
difficulties will arise from setting f;(E) =f2(E) 
=f(E). This has been done in what follows. 

The fraction of dn;, leaving the solid from a 
depth x may be expressed by dm, exp (—x/l,.), 
where /;, is assumed to be a function only of 
the energy E of the secondary electrons at the 
point of production. 

exp (—x/l;,) is the fraction of the 
incident primary electrons (mo) that penetrate 
to a depth x where 0Sx=x,. 

For the carbon backing (which is very thick 
compared with the magnesium layer), where 
x,:=x=0 there are the quantities similar to 
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Fic. 4. Secondary electron emission ratio 6 as a function 
of the layer thickness x; using an 8000-volt primary 
electron beam. 
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Fic. 5. Secondary electron emission ratio 6 as a function 
— layer thickness x; using a 6000-volt primary electron 
m. 
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those for the layer, namely, 
dn2,= — Nz pk2f2(E)dEdx 
and the fraction of dmz, which should get out 
will be 
exp (—x1/l1.) exp 


where /s, is also a function of E. 


N2p= No Exp exp [—(x—21)/lep ] 


is the fraction of m» incident primary electrons 
that penetrate to a depth x where SxS~. 

On the basis of the above formulation, the 
number of emitted secondary electrons per 
square cm should be 


exp (—x/l.)dm. 


r=0 


+f exp (—*1/l) 


Xexp [—(x—2x1)/l2s 
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Fic. 6. Secondary electron emission ratio 6 as a function 
hg layer thickness x; using a 4000-volt primary electron 
m. 
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Fic. 7. Secondary electron emission ratio 6 as a function 


fee layer thickness x; using a 2000-volt primary electron 
m. 


x carbon backing 


/ayer 


Fic. 8. Magnesium layer (0Sx;x) eva ted upon car- 
bon backing 


In order to secure the desired information 
about the behavior of the secondary electrons, 
it is necessary to make an assumption about the 
nature of the function /;,(Z). The assumption 
made tentatively is that elastic scattering occurs 
such that the mean distance between collisions 
depends on the energy E in the following way: 


hi 


where E is the energy of the secondary electrons 
at the point of production, and /»;, is a mean free 
path for secondaries of one-volt energy. 

The use of this E? relation probably weights 
or favors the high energy electrons because of 
screening effects in the collision process. This 
statement means that 1;,=/o:,E? results in 
values of /;, that are larger than they should be 
for values of E above several hundred volts. 
This fact is partially compensated for by the 
fact that the shape of the curve used for f(£) is 
such that it does not favor the high energies 
quite enough. 

On inserting 1,,=19:,E2 in the relation for m,, 
the number of secondary electrons leaving the 
emitter is 


Xexp (—x/louE?) exp (—x/hy)f(E)dxd E 


+ f —nok2 exp (—x/lip) 


E=0 
Xexp ] exp 

Xexp ]f(E)dxdE, 

where /.,=Jo2,E? as has been used; the argument 

for its use is the same as that for using /;,=/o1.E”. 


Because of the relatively high energy pri- 
maries used 1/l;, is considered negligible in 
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comparison with 1//;,, and 1/l2, is neglected in 
comparison with 1/l2,. Then 


$= J J ki exp (—x/lowE®)f(E)dEdx 


+f f ke exp (—x1/lo1.E?) 
Xexp ]f(E)dEdx. 
On integrating with respect to x 
5= f E*(1—exp (—21/louE?) 
0 


+halon FE? exp 
0 
and rearranging 


f 
0 


f E°f(E) exp (—x1/louE2)dE, 
0 
or 


f 


+B f “EY(E) exp (—%1/lo1,E*)dE. 


This expression relates the secondary emis- 
sion ratio 6 to the layer thickness x;. A reason- 
ably correct energy distribution of the secondary 
electrons at the point at which they are pro- 
duced in the solid together with data giving 6 


an 


10048 


200 100 400 
E-VOLTS 


Fic. 9. Energy distribution of secondary electrons produced 
in dx in energy range between E and E+dE. 


as a function of x; will permit the evaluation of 
lox, which in turn allows the calculation of /;,, 
the range of secondary electrons of energy E. 

The type of analysis undertaken is possible 
even approximately only because f(Z) for high 
primary energies changes very little for large 
changes in primary energy. Since f(E) cannot 
be measured directly, a calculated f(Z) was used. 
The best energy distribution at hand seemed to 
be that given by Bethe? for ionization of hydro- 
gen. From the calculated values given in refer- 
ence 4 a block distribution was drawn, and a 
curve was put through the block distribution. 
This was done for 1000-volt and 10,000-volt 
primary energies, and the point of interest is 
the small change in the distribution of energies 
when the primary energy is changed from 1000 
volts to 10,000 volts or even to 100,000 volts. 
The form used for the distribution function f(£) 
is shown in Fig. 9. 

By using the calculated f(£) mentioned and 
the following numerical methods, values of 6 are 
obtained from 


56=A f "EY 
0 


+B f E%f(E) exp (—C/E®)dE, 
0 


where 


C= X1/lots. 


The integrals are evaluated by means of 
plotting the integrand of each integral for a 
given value of the parameter C against the en- 
ergy E and using a planimeter to determine the 
area under each plot. The first integral can be 
evaluated from a plot made with the parameter 
C equal to zero. The second integral must be 
evaluated a number of times for a properly 
chosen range of values of C. The integral 


f " BY(E) exp (—C/E24E 


will be denoted by F(C); F(C) was evaluated for 
ten values of C ranging from 0 to 1.1105. The 
plots of E?f(E) exp (— C/E?) as a function of E 
are shown in Fig. 10. 


4H. Bethe, Handbuch de Physik 24, page 517, Table 31. 
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Fic. 10. The integrand used in the analysis as a function of the energy E of the 
secondary electrons. 


With F(C) evaluated in the desired range of C 
the constants A =kylo;, and B= —(kilo1s —Raloas) 
can be evaluated by using the values of 6 for 
x,=0 in which case only the backing material 
contributes to 6 and the value of 6 for x:=2.2 
X10-* cm where the value of 6 is determined 
almost entirely by the layer of magnesium. 

Values of F(C) for the range of C used are 
shown in Table I together with the values of 
Rilois, Relozs, the calculated 6, and the correspond- 
ing value of x, taken from the experimental curve 
for the particular primary beam voltage con- 
sidered. Calculated values of 6, (backing) and 
6, (layer) are also included in Table I. 

With Rilois and Roloss determined, the 6 equa- 
tion in the particular case of a 2000-volt primary 
beam turns out to be 


5=1.08—2.3010-*F(C), 
where 


F(C)= f “EY(E) exp —C/E*dE. 


Since F(C) has been evaluated for a desirable 
range of C, 6 can be calculated for this range of 
C. From the 6~x; curve for 2000 volts the value 
of x; can be found which corresponds to this 
particular value of 6 just calculated. With the 
value of C corresponding to the value of F(C) 
used to calculate this 5, a value of Jo, is found 
thus: 

C=xi/lo. or low=x1/C. 


6, and 6, in Table I were obtained by consider- 


ing the 6 equation 


-f exp (—Ca/ 


+ (talon) EY(E) exp 
0 
Now as C,— (i.e., x1) it is apparent that 
f E*f(E)dE, 
0 


which is the 6; due to the layer alone, since the 
backing cannot contribute for an_ infinitely 
thick layer. When C,- 0 (i.e., 


5—>(kalons) f EY(E)dE, 


which is the 6, due to the backing only since 
there is not a magnesium layer when x;=0. 

The values for the curves marked “‘layer”’ and 
“backing” in Figs. 4, 5, 6, and 7 were calculated 
by using 


= (Relozs) f exp (—C,/E*)dE; 


then we set 6;=5—6, to calculate 6; which can 
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TABLE I. Values of F(cn), Rilois, 8, Sjayer for the four primary energies used. 


2000 4000 6000 8000 

0 0 2.05 X105 
1 2.5 X10? 2.03 0.61 0.59 
2 15 X10% 1.92 0.64 0.56 0.08 2.5X10-*|0.51 0.44 0.07 2.5X10-*/0.47 0.41 0.06 3 X10-*| 0.44 0.38 0.06 4 X10-6 
3 4 1.80 0.67 0.52 0.15 5 X10*/0.54 0.41 0.13 5 X10-*/0.50 0.39 0.12 8 X10-*| 046 0.35 0.11 8 
4 X108 1.69 0.69 0.49 0.20 7 X10*/0.57 0.39 0.18 8.5 X10-*/0.53 0.36 0.17 1.4X107| 0.49 0.33 0.16 1.5107 
5 1.25X104 1.50 0.73 0.44 0.29 1.1X10-§|/0.62 0.34 0.28 1.6X1075/0.58 0.32 0.26 3 0.53 0.29 0.24 3 
6 2 1.33 0.77 0.39 0.38 1.8X10-§|0.67 0.31 0.36 2.6X10|0.62 0.29 0.34 0.58 0.26 0.32 5.7X107 
7 4 X10* 1.0 0.85 0.29 0.56 4 X10-§/0.75 0.23 0.52 6.8X10-§|0.70 0.21 0.49 1.75 0.66 0.20 0.46 1.8X1076 
8 6 X104 0.79 0.90 0.23 0.67 1.4X10*/0.81 0.18 0.63 2 X1074 
9 1.1 X105 0.48 

ki lows 5.2 X10-¢ 4.96 4.64 X10-* 4.46 X107* 

2 lore 2.9 X10-* 2.29 2.14 1.95 X10-* 


also be obtained from 
6 = Rulois E dE 
1= (kilou) J f(E) 
f E¥f(E) exp | 
0 


The constancy of 1);, over the range of C or x 
is considered to be a justification of the assump- 
tions made heretofor, namely, (1) /:, independent 
of x; (2) lis (3) the form of f(EZ) shown in 
Fig. 9. 

For the range of C from C; to C; the variation 
of Jo, for 2000-, 4000-, 6000-, and 8000-volt 
primary beams is shown in Table II. The reason 
for not giving values of Jo, for C; is that there is a 
large uncertainty in attempting to pick x; from 
the experimental curve because x; in this case 
is very small. When C=Co, x1=0 and only the 
carbon backing is present. The determination of 
lors beyond C; requires more data (except for 
2000 volts). 

The range of x; covered is that from 310-° 
cm to 1X10-‘ cm corresponding to the range of 
C from C, to 

Figure 11 shows how /9;, varies with x;. The 
6~x1 curves have been plotted for the four 
primary energies to compare the data (curves) 
of Figs. 4, 5, 6, and 7 with the values of 6 calcu- 
lated (points) from the equation 


6=AF(Co)+BF(C), 


which has been discussed in detail. 

The four values of Jo:, used for this calculation 
are the average values (averaged over the range 
of x; for the particular primary energy) given in 
Table II. 


The 6~x; curve of Fig. 11 for 2000-volts 
primary energy shows that the calculated points 
fall below the curve for the thinner layers while 
for thicker layers the calculated points will be 
above the curve. The fact that some points lie 
below the curve indicates that in this range of x; 
the theory does not account for enough secondary 
electrons, which is another way of saying that /;, 
is too short in this range of x1. 11,=/01,/? does not 
favor the electron paths enough in this range of 
x, which may mean that a modified form of the 
relation between /;, and E would represent the 
facts better than the relation we have used or it 
may mean that /;, is not independent of x as we 
have assumed it is. The relation 1,,=/9:.E? 
appears to yield paths somewhat too long for 
larger x; which means that /;, is too long—that 
the energy variation lies between the square of 
E and the first power of E. 

When the primary energy is increased to 4000, 
6000, and 8000 volts the agreement between 
curves and calculated points becomes increas- 
ingly better at smaller x; and somewhat worse 
for larger x;. This better agreement at small x; 
for increasing primary energy indicates that 
l1s=1:,E2 is about the right value of /;, for the 


TaBLe II. Values of /:, for the four primary energies used. 


n 2000 4000 6000 8000 
2 1.7X 107° 1.7X10-° 2 x10-° 10-° 
3 1.3 1.3 2 2 
4 1.0 1.2 2 2.1 
5 0.9 1.3 2.4 2.4 
6 0.9 1.3 2.6 2.8 
7 1.0 1.7 
8 2.3 3.3 
1.3 1.7 2.2 2.3 
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higher primary energies. Now for small x; and 
at 2000 volts /;,(E) is too small and at 8000 
volts about right so that it might be said that 
there are more slow secondary electrons at 
2000 volts than at 8000 volts and that /;,(£) does 
not allow sufficiently large values of /,, for the 
slow secondary electrons. 

Another factor involved is the shape of the 
distribution f(Z) which probably falls off too 
rapidly at high energies. f(Z) is multiplied by 
E? which enhances considerably the effect of the 
high energy part of f(Z) in the numerical integra- 
tion. /,, should be a measure of the depth from 
which the secondary electrons can on the average 
just get out of the surface. The depth of origin 
of secondary electrons in the energy range from 
10 volts to 200 volts will vary from 210-7 cm 
to 4X10-* cm, a factor of 200 in the range for a 
factor of 20 in the secondary energy. 

For a surface having a relatively large yield 
lors might be larger than J»), for a low yield sur- 
face such as the one being dealt with here, but it 
is more likely that the dependence on secondary 
electron energy would be different. No attempt 
was made to secure data for a high yield surface, 
but it would appear that in a good secondary 
emitter the dependence of range on the energy 
must be such as to allow abnormally long 
range at low energies—an effect equivalent to a 
kind of Ramsauer effect. Such an effect, if 
present, could be measured by measuring /o., 
and comparing with J, for a relatively poor 
emitter such as magnesium. It is likely that 
with a pronounced Ramsauer effect the present 
analysis would break down because /;,=101,E? 


VOLTS 
ENERGY 
0.8 > 
4000 
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000, 
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Fic. 11. Comparison of data and calculated values of the 
secondary emission ratio 6 as a function of layer thickness 
2%; for the four primary energies used. 


would no longer apply even approximately. 
lox. would fail to be constant as soon as /;,(Z) 
failed; so measurement of Jo:, for a high yield 
surface would tell something about the manner 
in which the slow electrons are favored in getting 
out of the surface. 


Electron Ranges in High Yield Surfaces 


It is believed that the effect of oxygen in 
producing a good secondary emitter is to in- 


CARBON 
0.7 
0.6 
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Fic. 12. Secondary electron emission ratio 6 as a function 
of primary electron energy V, for the carbon backing 
material used. 


crease the range of the slow electrons. One 
possible qualitative view has been given by 
Bethe.® It is known that magnesium alone is a 
very poor emitter. It is also known that the 
straight magnesium oxide without occluded 
oxygen is a poor emitter. When, however, mag- 
nesium or magnesium oxide is evaporated through 
an atmosphere of oxygen in the proper way a 
mixture results which is not MgO alone but 
MgO with perhaps some free magnesium and 
probably with some free oxygen. MgO is a 
ceramic and a good insulator, and, if bombarded 
with a primary beam, will simply charge up and 
stop emitting (particularly with low energy 
primary beam). If, however, some* impurity 
such as oxygen is introduced, impurity energy 
levels are introduced between the allowed 
empty band and. the allowed filled band of the 
insulator, thereby permitting transitions from 
the filled to the empty band. If few impurity 
levels are introduced, the insulator becomes a 
poor semiconductor. The presence of these few 
impurity levels permits poor conduction so that 
the bombarded surface cannot charge up to any 
great extent. At the same time the few impurity 


5 In a conversation with Professor H. Bethe. 
6 The actual amount required may be extremely small. 


10-8 
10-8 
10-8 
10-4 
Its 
nts 
rile 
lie 
ary 
t Lis 
of 
the 
the 
r it | 
we | 
for 
hat 
of | 
00, 
2en 
as- 
rse | 
hat 
the 
sed. 
j H 
| 
‘ 


348 W. F. BRUCKSCH, JR., 


levels do not provide many levels for electrons 
to fall into in transition from the filled band 
to the empty band, thereby improving the 
chances of completion of the transition. This 
state of affairs will permit a continuous supply 
of electrons to the conduction band for emission, 
and at the same time provide a minimum number 
of levels to trap possible secondary electrons. 
The desirability of having the bottom of the 
allowed empty band as high as or higher than the 
potential energy of an electron in vacuum out- 
side the emitter has been pointed out by Bruin- 
ing and De Boer.’ 

Many investigators have found that the so- 
called oxides of barium and magnesium (pre- 


7H. Bruining and J. H. De Boer, Physica 6, 836 (1939). 


AND W. T. ZIEGLER 


pared by raising the temperature of the metal 


in the presence of oxygen) have maximum 
secondary yields in the neighborhood of three 
or four. Strictly speaking, these are not simple 
oxides; they are, in the writer’s opinion, oxides 
with too little or too much impurity to allow 
them to be really good secondary emitters. 
They are, however, better emitters than con- 
ductors are. 

The author expresses here his appreciation for 
the encouragement and advice of Professor 
Lloyd P. Smith who suggested this problem. 
The writer is also indebted to Professor Hans 
Bethe for valuable discussions, and to Professor 
D. H. Tomboulian for assistance with evapora- 
tion problems. 
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Studies in Superconductivity * 
II. Evaporated Lead Films 


W. F. Bruckscn, Jr.t AND WALDEMAR T. ZIEGLER 
Cryogenic Laboratory, Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 


(Received June 29, 1942) 


The electrical and magnetic properties of evaporated lead films of several thicknesses 
(1000-—3000A) have been studied in the superconducting state. The films exhibited transition 
temperatures (in zero magnetic field) of 7.23+0.03°K irrespective of thickness. The disap- 
pearance of resistance (99 to 1 percent) usually occurred within a range of 0.10°. The depression 
of the transition temperature with increase in measuring current J, was found to be greater 
the thinner the film, all films showing larger depressions than a 254-micron lead wire 
(T.=7.20°K). The depression of the transition temperature in several magnetic fields (0-80 
oersteds) perpendicular to the plane of the film was found to depend upon thickness, the 
thinner films requiring larger external fields to produce a unit depression. The J, versus T-. 
curves obtained exhibit inflection points not heretofore observed, while the H, versus T. 
curves obtained suggest the existence of a similar phenomenon. The effect of measuring 
current on the depression of the transition temperature in constant magnetic fields also has 


been studied. 


N the course of the development of a super- 
conducting bolometer suitable for the meas- 
urement of small amounts of infra-red radiation! 
the properties of evaporated lead films of several 
thicknesses and tin, lead, columbium, and 


*From a dissertation submitted to the Board of Uni- 
versity Studies of The Johns Hopkins University by 
W. F. Brucksch, Jr. in conformity with the requirements 
for the degree of Doctor of Philosophy. 

mi Present address: U. S. Rubber Company, Detroit, 
ichigan. 

1 Andrews, Brucksch, Jr., Ziegler, and Blanchard, Rev. 
Sci. Inst. 13, 281 (1942), hereafter referred to as Paper I. 


tantalum wires of various sizes were investigated 
in the superconducting state. It is the purpose 
of this and a subsequent paper to present some 
of the results obtained in these studies, since we 
believe them to have an interest quite apart from 
their original intent.? In general our studies were 


2 Some of the material to be discussed has already been 
ang before meetings of the American Physical 

iety. See Brucksch, Jr., Ziegler, Blanchard, and 
Andrews, Phys. Rev. 59, 688 (1941); Brucksch, Jr., 
Ziegler, Horn, and Andrews, ibid., 60, 170 (1941) 
(Abstracts). 
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confined to the measurement of the change of 
the electrical resistance of the test specimen 
with temperature in several magnetic fields 
(0-80 oersteds) and with various measuring 
currents. From these measurements we were 
able to deduce the change of the superconducting 
“transition temperature” 7 with the measuring 
current J. (the critical current) and the applied 
external magnetic field H, (the critical field). 
Throughout this paper the transition tempera- 
ture J, is taken to be that temperature in the 
transition range at which the resistance of the 
sample is equal to one-half that which it is in 
the normal state just above the transition range. 


EXPERIMENTAL PART 


The lead films were prepared by evaporation 
of the metal from a heated tungsten filament in 
a high vacuum onto soft glass bases at room 
temperature. These glass bases were provided 
with two tiny sealed-in platinum leads to which 
current and potential leads could be attached. 
The resistance of the film was followed during 
deposition and various thicknesses obtained by 
varying the time of deposition. 

The soft glass bases were of two sorts: (a) flat 
ribbons and (b) glass ‘“‘bows.’’ They were made 
by drawing an ordinary microscope slide into a 
thin ribbon approximately 1 mm wide and 0.5 
mm thick. The “bow” type of film support 
was designed for the aforementioned radiation 
studies,’ and arranged so that the glass ribbon 
(on which the film was to be deposited) was 
sealed to the glass base only at the ends. 

After the films were deposited at room 
temperature they were inserted in the cryostat 
and usually cooled to about 80°K the day after 
they had been prepared. 

The low temperature cryostat used in these 
studies has been described in Paper I and will 
not be discussed in detail here. The sample to 
be studied was mounted on the top of the copper 
helium vessel in a horizontal position. The 
current and potential leads were soldered to 
wires passing around the massive helium vessel 
and thence to the outside. 

Three methods of mounting were used. 
Method I: The glass base on which the film was 


3 Andrews, Brucksch, Jr., Ziegler, and Blanchard, Phys. 
Rev. 59, 1045 (1941). 


deposited was suspended between binding posts 
by means of the electrical leads. These binding 
posts were attached to (but insulated electrically 
from) a copper plate which was then soldered to 
the helium vessel. Method II: The glass base on 
which the film was deposited was cemented 
directly to the copper plate which was then 
inserted as in Method I. Method III: The film 
was mounted as described in Method I. However, 
the copper plate was supported on glass rods in 
such a manner as to afford poor thermal contact 
with the helium vessel. (This arrangement was 
later used for radiation studies.) 

We believe that the thermal contact afforded 
by Method II may be described as very good, 
by Method I as moderately good, and by 
Method III as poor. 

All electrical measurements were made with 
White and Type K potentiometers and sensitive 
galvanometers manufactured by the Leeds and 
Northrup Company. Potential measurements 
were made to +0.2 microvolt. 

The measurement of the temperature of the 
test sample was carried out by means of a 
constantan resistance thermometer (R= 6242 at 
4.2°K) wound on the outer surface of the helium 
vessel. Measurements of potential permitted a 
relative temperature‘ measurement to be made 
with an accuracy of +0.008°. The resistance 
thermometer was calibrated with a helium gas 
thermometer of the type described by Wood- 
cock. This acted as a helium vapor pressure 
thermometer below 2.5°K. Using the normal 
boiling point of helium (4.22°K) as a fixed point 
we obtained the following results® for the 
transition temperatures of tin, lead, and colum- 
bium metals: 7.(Sn)=3.60°; .(Pb)=7.20°; 
T (Cb) =9.58°. These may be compared with the 
generally accepted values for these transitions: 
T (Sn) =3.73°;7 T.(Pb)=7.2°;8§ and T7.(Cb) 
=9.3°.° Inasmuch as we did not correct for gas 


‘ All temperatures reported in this paper refer to these 
relative temperatures, since their accuracy relative to one 
another is much greater than their absolute accuracy. 
If absolute temperatures are desired cognizance must be 
taken of the errors discussed below. 

&’ A. H. Woodcock, Can. J. Research A16, 133 (1938). 

6 These are discussed more fully in Paper III of this series. 

7 Silsbee, Scott, and Brickwedde, J. Research Nat. Bur. 
Stand. 18, 295 (1937). W. Tuyn and H. K. Onnes, Leiden 
Comm. 181 (1926). 

8H. K. Onnesand W. Tuyn, Leiden Comm. b160 (1922). 

® Meissner, Franz, and Westerhoff, Ann. d. Physik 17, 
593 (1933). 
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imperfection (which would tend to make our 
temperature scale too high above 4.2° and too 
low below 4.2°) we have taken the above results 
as indicating the error in our absolute tempera- 
ture scale. 

The magnetic field was produced by a solenoid 
(22 cm long; 12 cm mean diameter) immersed in 
liquid nitrogen. The solenoid was wound in such 
a manner that the field produced was perpen- 
dicular to the plane of the film, the film being 
located approximately at the center of the 
solenoid. By varying the exciting current from 
0 to 200 milliamp. fields of 0 to 80 oersteds 
could be produced.” 


RESISTANCE VS. TEMPERATURE MEASUREMENTS 


After the film had been mounted but prior to 
cooling the cryostat a measurement of the 
resistance of the sample was made at room 
temperature (~25°C). This has been called Ro. 
Additional values of the resistance were obtained 
at the boiling point of liquid N2 (77°K), liquid 
Hz (20.4°K), and the lowest solid hydrogen 
temperature (~ 14°K) reached prior to expansion 
of the helium gas." 

In the temperature range 12.5-4.2°K the 
helium vessel was cooled in steps by partial 
expansion of the enclosed gas. When the sample 
being studied was found to have lost its re- 
sistance (i.e., become superconducting) cooling 
by expansion was stopped. The helium vessel 
was then allowed to warm up gradually (warming 
rate~0.05° per minute) which it did due to the 
heat leaks into the apparatus from the environ- 
ment (12-13°K). Measurements of temperature 
(resistance thermometer) vs. time and resistance 
of film vs. time permitted a determination of the 
resistance of the film as a function of temperature. 

Critical field studies were made by adjusting 
the external field to the desired value, cooling 
the sample until it became superconducting, and 
then following the above procedure for obtaining 
the R vs. T curve. 

Critical current studies were carried out by 
first cooling the sample well below its normal 
transition temperature and then adjusting the 

10 These fields were calculated from the expression for 


the field at the center of a solenoid whose length is several 


times its mean diameter. . 
1 Liquid helium was produced by the Simon expansion 
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current through the film to some value below 
that required to destroy superconductivity. The 
helium vessel was then allowed to warm slowly 
(~0.03-0.05° per min.). As the temperature of 
the sample approached the transition tempera- 
ture corresponding to this current the galva- 
nometer beam began to show the restoration of 
a small fraction of the total resistance (perhaps 
0.1-1 percent), and then within an interval of 
~0.01° the remainder of the resistance appeared. 
In order to prevent damage to the film as well 
as to the galvanometer the restoration of the 
normal resistance was used to operate a relay 
which cut off the current within the space of 
0.1-1 second. 

After superconductivity had been destroyed 
by a given current, the current was decreased, 
the circuit reestablished, and the cycle repeated 
on the same ‘“warm-up.”’ In this way the 
variation of transition temperature with critical 
current was obtained. 


EXPERIMENTAL RESULTS 


Table I summarizes the properties of and the 
results obtained for lead films of various thick- 
nesses, with small measuring currents (~10~ 
amp.) in zero applied magnetic field. 

Film thicknesses were obtained in two ways. 
An electrical thickness was calculated from the 
length and width (average dimensions: 16 mm 
and 1.5 mm, respectively) and electrical re- 
sistance of the film at room temperature. In this 
we assumed the metal to have normal resistivity 
and density. A chemical thickness was determined 
after the electrical experiments were completed 
by chemical analysis” for the mass of the film 
studied. By assuming a normal density this 
method gave thicknesses with an estimated 
maximum error of 20-25 percent. 

On the basis of the chemical thickness studies 
the films fall into three groups: 3000A, 1400A, 
900A. An examination of the electrical thick- 
nesses shows three corresponding groups: 2300A, 
400A, and 70A. With one exception the R/Ro 
ratio just above the superconducting region 
(column 7, Table I) shows corresponding groups." 

2 The method employed was a variation of that de- 


scribed by J. H. Yoe, Photometric Chemical Analysis, Vol. 1 


(Wiley, 1928), p. 250. 
13 Tt seems reasonable to assume that the thinner the film 


the larger would be the R/Ro ratio just above the transition, 


SUPERCONDUCTIVITY LEAD FILMS 351 


TABLE I. Summary of data on lead films. 


Method of R/R 
Film support. Re(ohm) 80°K 14°K 7.3% arial 10-00% 

7A* I 12.98 0.2506 0.0486 0.0304 0.0198 2100 3300 7.24 0.08 0.02 

7B* I 10.81 0.2556 0.0534 0.0347 0.0237 2600 2700 7.25 0.10 0.03 

9A*® I 84.51 0.2747 0.0754 0.0553 0.0459 380 1600 7.24 0.08 0.06 

9B* I 59.47 0.2692 0.0711 0.0501 0.0406 550 1200 7.20 0.03 0.01 
10A* II 405.0 0.4492 0.2933 0.2789 0.268 70 1000 7.21 0.11 0.05 
10B* II 451.3 0.4601 0.3052 0.2909 0.279 70 900 7.22 0.18 0.05 
13A% III 42.41 0.2664 0.0684 0.0509 0.0419 460 1400 6.93 0.11 0.08 
13B> III 51.24 0.2621 0.0635 0.0407 0.0358 380 1200 7.02 0.10 0.06 
14A> III 109.6 0.2566 0.0656 0.0473 0.0367 140 1500 7.11 0.10 0.07 
14B> III 499.7 0.2593 0.0681 0.0513 0.0411 50 1100 6.94 0.13 0.06 
Bulk Pb? 0.221 0.0269 0.0101 7.24 

* Ribbon base. 2 The tem ture ran 


1 Temperature at top of transition range. 


We conclude, therefore, that regardless of which 
criterion of thickness is used the relative thickness 
of the three groups is approximately correctly 
given. 

In zero magnetic fields the transition tempera- 
tures T. of films 7A-10B were 7.23+0.03°K 
(Table I). The R vs. T curves for these films 
were smooth curves with a transition range of 
the order of 0.05—0.15°. These observations are 
in accord with those of Onnes and Tuyn"™ as 
well as measurements made by us for lead wires'® 
and indicate that our films behaved like bulk 
metal in these respects. 

On the other hand, films 13A-14B, which 
were placed in the cryostat immediately after 
preparation and were cooled to 80°K, within 
4-5 hours thereafter exhibited transition temper- 
atures definitely below 7.2° as well as transition 
curves of a somewhat more irregular form. The 
poorer thermal contact between the films and 
the helium vessel afforded by the method of 
mounting III) may be partly responsible, but 
since the transitions were always determined 
during warming, the values for 7. reported in 
Table I are to be regarded as upper limits. 
However, we believe that the anomalies observed 
with these films are due primarily to the fact 
that at room temperature annealing of the film 
is a rather slow process,'* and hence the films 


because of the effect of strains within the metal lattice 
arising (a) from the surface forces of the substrate and (b) 
from imperfect annealing. These effects are difficult to 
disentangle. See Appleyard and Bristow, Proc. Roy. Soc. 
172A, 530 (1939). 
4H. K. Onnes and W. Tuyn, Leiden Comm. b160 (1922). 
18 See Paper III of this series. 


«Onnes and Tuyn, Leiden Comm. b160 (1922). 


13A-14B had not been allowed sufficient time 
to come to equilibrium." 

Figure 1 presents a summary of H, vs. T, for 
various film thicknesses (chemical), including 
two of the ‘‘unannealed”’ films. Figures 2 and 3 
present data for J, vs. T, in the absence and 
presence of a magnetic field, respectively. The 
magnetic fields were usually limited to the 
values 0, 40, and 80 oersteds, the critical currents 
to 100 milliamperes or less. In each of these 
figures experimental data for a 25-micron lead 
wire is included for comparison. 


DISCUSSION OF RESULTS 


An inspection of the curves in Fig. 1 shows, 
if the unannealed films be excluded, that the 
thinner the film the greater is dH./dT,. As the 
film thickness becomes larger (3300A) the slope 
much more nearly approaches that for 25-micron 
wire (dH./dT.2170 oersteds per deg.) which 
approximates that for bulk lead (dH ./dT.2120 
oersteds per deg.).? Similar results have been 
observed for annealed evaporated mercury film 
by Appleyard, Bristow, London, and Misener."” 

The experimental results with critical currents 
in zero magnetic field (Fig. 2) show that the 
thinner the film the smaller the critical current 
required to produce a given lowering of the 


16 See, for example, R. Suhrman and H. Schnackenberg, 
Zeits. f. Elektrochemie 47, 277 (1941). 

17Similar results have been obtained for unannealed 
mercury films by Appleyard, Bristow, London, and 
Misener, Proc. Roy. Soc. A172, 540 (1939). On the contrary 
Shalnikov (Nature 142, 74 (1938)) has found an increase 
in transition temperature for unannealed lead films de- 
posited at 4.2 K. 
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70 
CRITICAL TEMPERATURE “K 


Fic. 1. Variation of transition temperature 7, with 
applied magnetic field for lead films of several thicknesses. 
The 25-Pb wire curve has a third point (not shown) at 


6.8°K for 84 oersteds. 


transition temperature. By comparison with the 
lead wires it is seen that the films are much more 
sensitive to current in the superconducting state 
than is a wire of the same metal. The resistance 
was always restored in a very abrupt fashion, 
the transition range being of the order of 0.01°. 
This effect was independent of the method of 
mounting the film in the cryostat. Shalnikov'’ 
has observed similar effects for evaporated lead 
and tin films. 

Figure 3 shows the results obtained when 
critical current studies were made in constant 
magnetic fields. An inspection of the curves 
shows that the combined effect of current and 
field results in a greater lowering than one would 
expect from a combination of the data given in 
Figs. 1 and 2. Moreover, the excess lowering of 
T. is larger the greater the magnetic field being 
~0.06° and ~0.14° for fields of 42 and 84 
oersteds, respectively. Hence it is seen that a 
much smaller current is required to destroy 
superconductivity in a film which is in a magnetic 
field than would be predicted on the basis of a 


AND W. T. ZIEGLER 


combination of H, vs. T, and I, vs. T. curves of 
the sort given in Figs. 1 and 2. 

A consideration of the curves in Fig. 1 shows 
them to be concave upward, although the 
thickest film seems to have an inflection point in 
its HT, vs. T, curve. Likewise, with low critical 
currents the J.—7. curves (Fig. 2) are all 
concave upward, but show inflection points. 

We have been unable to find any published 
data for films which indicate the existence of an 
inflection point in the H, vs. T. and I, vs. T, 
curves. Extensive studies of mercury films in 
varying fields down to 100 oersteds by Apple- 
yard, Bristow, London, and Misener gave curves 
which were concave downward.'® This is true 
generally for bulk metal samples,'® although a 
slight upward concavity of H, vs. T. has been 
found for some samples of tantalum wire in 
low fields.*° 

Shalnikov" has made critical current studies 
for very thin lead and tin films. The J, vs. T, 
curves reported by him are concave downward, 
but the measurements are confined to tempera- 
tures below 4.2° where large currents (0.15 amp.) 
are required to destroy superconductivity. Below 

=7°K our curves are also concave downward. ¢ 


CRITICAL TEMPERATURE 


Fic. 2. Variation of transition temperature 7, with 
measuring current (in zero applied magnetic field) for 
lead films of several thicknesses. 


18 Reference 17. These investigators give data for a 
single Hg film 2410A thick which permits one to deduce 
the shape of the H. vs. T. curve for fields from 60-130 
oersteds. Within this range the curvature is still concave 
downward. 

19See D. Schoenberg, ad (Cambridge 
University Press, 1938), p. 107. 

20K, Mendelssohn and J. R. Moore, Phil. Mag. 21, 532 
(1936); Silsbee, Scott, an Brickwedde, reference 7. 
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The simultaneous occurrence of the upward 
concave nature of the H, vs. T, and IJ, vs. T, 
for low magnetic fields and critical currents is 
consistent with the generally accepted view that 
the action of the critical current in destroying 
superconductivity is to be associated with the 
magnetic field produced by the current. Hence 
the occurrence of the effect in the two sets of 
measurements arrived at in entirely different 
ways leads us to believe the effect to be real. 

Several possible objections to this conclusion 
may be raised. Conceivably the constantan 
resistance thermometer used for temperature 
measurements might have exhibited an anoma- 
lous resistance in the range 7.0-7.25°K.?! How- 
ever, a plot of potential across the thermometer 
vs. time during a slow ‘‘warm-up” was always 
perfectly smooth in this range. Hence we 
conclude that no such anomaly exists here. 

Also heating effects on the one hand might be 
assumed to invalidate the measurement of J, 
vs. T, while the assumption that the applied 
field varied linearly with applied solenoid current 
might also be incorrect. The fact that the 
observed slope of the H, vs. T, curve for Pb 
and Sn wires (unpublished results) is constant 
(i.e., linear curve) shows that the effect is 
associated with the films and not with an 
erroneous temperature or magnetic field scale. 
In addition since the method of mounting the 
films (Methods I, II) did not affect the presence 
of the anomalous curvature the argument that 
heating effects are responsible loses much of its 
force. We conclude, therefore, that the effects 
observed are to be associated with the films 
themselves and not with the method of making 
the measurements. 


%1 See Onnes, Leiden Comm. b219, 13 (1932). 


CRITICAL TEMPERATURE 


Fic. 3. Variation of transition temperature 7, with 
measuring current for Pb film 7A in several constant 
magnetic fields. 
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Studies in Superconductivity 
III. Sn, Cb, Ta, and Pb Wires 


W. T. ZreEGLER, W. F. JR.,* AND J. W. HICKMAN 
Cryogenic Laboratory, Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 
(Received June 29, 1942) 


The electrical and magnetic properties (in fields of 0-80 oersteds) of Sn, Cb, Ta, and Pb 
wires of several diameters (25-250 microns) have been studied in the superconducting state. 
In addition x-ray powder photographs (made with copper Ka radiation) have been taken 
of the wires studied. The experimental results obtained showed that a low resistance ratio, 
R/Ro (Ro=room temperature resistance), high transition temperature, and small transition 
range were always found together for a given metal. For lead the 25,127- and 254-micron 
wires studied, all exhibited a transition temperature of 7.20+0.01°K, a fact which is attributed 
to the similar crystalline character of these wires. 


N the course of an investigation of the proper- 
ties of attenuated superconductors with a 
view to their use as receiving elements in low 
temperature radiometry'? we studied several 
different metals in the form of very fine wires. 
In our search for a metal wire having the de- 
sirable characteristics of high resistance just 
above the transition region, a narrow transition 
range (i.e., large value of dR/dT) and a transition 
temperature in a temperature region easily ac- 
cessible to us (above 3°K) we tested the electrical 
and magnetic properties of fine filaments of tin, 
lead, columbium, and tantalum wires in the 
superconducting state before selecting the latter 
as the metal from which a bolometer was con- 
structed.” In general our studies were confined to 
the measurement of the change of the electrical 
resistance with temperature in several magnetic 
fields (0-80 oersteds) and with small (~10-‘ 
amp.) measuring currents. Certain of the results 
obtained are new, while others constitute a sub- 
stantiation of certain isolated experiments re- 
ported by other investigators. It is the purpose 
of this paper to treat of the results obtained for 
the above-mentioned wires. 


EXPERIMENTAL 


All the wires studied were obtained from com- 
mercial sources and were very pure according to 


* Present address: U. S. Rubber Company, Detroit, 
Michigan. 

1 Andrews, Brucksch, Jr., Ziegler, and Blanchard, Phys. 
Rev. 59, 1045 (1941). 

2 Andrews, Brucksch, Jr., Ziegler, and Blanchard, Rev. 
Sci. Inst. 13, 281 (1942). Brucksch, Jr., and Ziegler, Phys. 
Rev. 62, 348 (1942). These papers will be referred to as 
Papers I and II, respectively. 


the statements of these manufacturers.’ These 
wires varied in diameter from 25 to 250 microns. 
X-ray powder photographs were taken (with 
copper Ka radiation ; sample stationary) of each 
wire to determine its crystalline character as well 
as the presence of impurities. The results ob- 
tained are listed in Tables I, II, and III. Only the 
254-micron columbium and 230-micron tantalum 
wires showed the presence of impurities, faint 
lines for SiC appearing on each photograph. 

A 30-32 mm length of each specimen was 
mounted in the low temperature cryostat, pre- 
viously described, which could accommodate as 
many as four samples in a single experiment.‘ 
Each wire was suspended between massive bind- 
ing posts, to which potential and current leads 
were affixed.’ These supports were mounted on 
a copper plate, but were electrically insulated 
from it by thin varnished mica sheets. This plate 
in turn was soldered to the top of a heavy-walled 
copper helium vessel, thus bringing the wires 
into excellent thermal contact with the vessel. 

The temperature of this helium vessel was 


Extruded wires of Sn and Pb were obtained from 
Baker and Company, Newark, New Jersey. Hard-drawn 
wires of Ta and Cb were purchased from the Fansteel 
Metallurgical Corporation, North Chicago, Illinois. A 
special sample of 0.001” Cb sheet was procured from 
A. D. Mackay, 198 Broadway, New York City. 

‘For a description of the cryostat see Papers I and II. 
The operation, temperature control, electrical measuring 
devices, and temperature scale used in this research are 
discussed in these papers. o- 

5 Sn and Pb wires were soldered directly to the binding 
posts by means of woods alloy. Columbium and tantalum 
wires were mounted by first spot welding a small tab of 
iron ribbon to each of the wire extremities [see J. Strong, 
Procedures in Experimental Physics (Prentice-Hall, 1939), 
p. 145]. It was then possible to solder the wires into place. 
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TABLE I. Tin wires. 


K 
microns Crystallinity field 10% 1% 99% oersteds 10-90% 1-99% 
25 Small crystals 31.5 3.46 0.090 0.440 2.69 0.215 0.575 107 
254 Small crystals 26.9 3.60 .035 .070 2.81 .220 615 109 
220' 0.98 3.724 .023 116 
190? 0.96 3.735 .020 .045 1398 


1 Silsbee, Scott, and Brickwedde, J. Research Nat. Bur. Stand. 18, 295 (1937). 


2? Tuyn and Onnes, Leiden Comm. 181 (1922). 
3de Haas and Engelkes, Physica 4, 325 (1937) (for a single crystal). 


TABLE II. Columbium wires. 


R/RoX10° T°. K AT? AT° 
Wire R/RoX108 at bottom zero Transition range Transition range dH./dT 

diameter at top of of transi- R/Ro (H =84 R/Ro oersteds 

microns Crystallinity transition tion eld 10%-90%  1%-99%  oersteds) 10%-90% 1%-99% deg.~ 

28 Small crystals* 456 5.2 5.09 0.020 0.030 
254 Small crystals* 270 1.0 8.64 .005 .035 8.57 0.015 0.055 1200 
300X 25 Polycrystalline 1.6 .005 9.50 025 -120 1050 


* Powder photographs showed extra lines indentified as SiC. These lines disappeared when wire diameter was decreased 20 percent by etching, 


indicating a surface impurity. 


! Meissner, Franz, and Westerhoff, Ann. d. Physik, 17, 593 (1933). This sample was 99.9 percent Cb. Meissner and Franz, Zeits. f. Physik 63, 558 
(1930), also studied a columbium sample of P—~ 98.43 percent Cb. This specimen was a strip 600 X90u and had a transition range of 8.54-—8.69°K, 


The data for our 254-micron wire com 


varied in two ways. Below 5°K this was accom- 
plished by regulation of the vapor pressure of 
the boiling helium (produced by the Simon ex- 
pansion method). In the range 5-10°K we em- 
ployed the technique (described in Paper II) of 
expanding only a part of the highly compressed 
helium gas present in the above-mentioned 
vessel. This partial expansion produced partial 
cooling, but no actual liquefaction of helium. 
When a sufficiently low temperature had been 
reached such that the sample became supercon- 
ducting, expansion was stopped and the helium 
vessel allowed to warm up gradually, resistance 
vs. temperature measurements being carried out. 
Because of the large heat capacity of the helium 
gas usually remaining in the vessel, this rate was 
quite slow (~0.05°/min.). 

Temperatures were measured by means of a 
constantan resistance thermometer coiled about 
the helium vessel. The temperature so measured 
has been called the relative temperature and 
could be measured with~a precision of +0.008°. 
These relative temperatures, in turn were cor- 
rected by the aid of a calibrated helium gas 
thermometer of the type described by Wood- 
cock.® It is estimated that these relative tempera- 


* A. H. Woodcock, Can. J. Research A16, 133 (1938). 


very well with this sample as regards 7: and R/Ro. 


tures are correct to +0.05° at 3°K, +0.1° at 
7°K, and +0.25° at 9.5°K on the absolute tem- 
perature scale.’ 


EXPERIMENTAL RESULTS 


The experimental observations for tin, colum- 
bium and tantalum are summarized in Tables I, 
II, and III. The symbols used require a few 
words of explanation. Ro represents the resist- 
ance of the specimen at room temperature. The 
ratio R/Ry represents the fraction of this re- 
sistance which remains at+a given temperature 
in the transition range. The columns headed 
10-90 percent and 1-99 percent represent the 
temperature interval in which the ratio R/Ro 
for the sample changed by this percentage, 100 
percent representing the value of the ratio in the 
normal state just above the transition. Through- 
out this paper the critical temperature 7, is 
taken to be that temperature at which the re- 
sistance of the sample is one-half the value in 
the normal state just above the superconducting 
region. 

The variation of resistance with temperature 
was taken in several constant transverse mag- 


7For a more complete discussion of the temperature 
scale see Paper II. 
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TABLE III. Tantalum wires. 


Wire R/Ro X108 Te range dH-/dT- 
diameter at top of (H=0) 10-90% oersteds 
microns Crystallinity transition °K R/Ro deg.“ 
28 Polycrystalline 389 <3.19° -- 
127* Polycrystalline 163 4.11 0.030 770 
230 Polycrystalline? 91 4.28 010 333 
50! 35.3 4.38 .030 — 
75? 260 3.96 .050 1250 
127? 260 4.07 .055 1020 


* This sample was also tested when coiled as a helix on a mandrel 3 
times the wire diameter. Measurements showed 7. =3.49°, R/Ro X10° 
=124, 10 percent to 90 percent R/Ro =0.040°. The displacement. of Te 
was considered directly the result of introducing additional strain in the 
wire by cold working. 

> Powder photographs showed extra lines identified as SiC. These lines 
disappeared when wire diameter was reduced 20 percent by etching, 
indicating surface impurity. 

¢ This wire not superconducting at 3.19°K. 

1 Meissner and Franz, Zeits. f. Physik 63, 558 (1930). 
on Woy Scott, and Brickwedde, Jr. Research Nat. Bur. Stand. 18, 


netic fields (0-80 oersteds) with constant current 
flow in the wire. The average slope of the critical 
magnetic field (H.) function of temperature 
(dH./dT.) was taken graphically from the plot 


of H, versus 


The crystallinity of the samples is described as 
large crystals, small crystals, or polycrystalline 
depending upon the x-ray photographs obtained. 
The first term implies that the photograph con- 
sisted largely of spots, the second spots with 
some diffraction rings, the third typical diffrac- 
tion rings. 

Three lead wires 25, 127, and 254 microns in 
diameter were studied; x-ray powder diagrams 
(copper Ka radiation) showed all three wires to 
consist of large crystals. All three wires showed 
R/Ro ratios of 1.8X10-* just above the super- 
conducting region. The transition interval was 
very small for all wires (~0.04° for 1-99 percent 
R/R>) and the transition temperatures were all 
7.20+0.01°K in zero magnetic field. These results 
may be compared with those of Onnes and 
Tuyn.® These investigators found R/R,)=1.0 
X10-* and 7,=7.2°K for a very pure annealed 
lead crystal. In a field of 84 oersteds the 25- 
micron wire exhibited a somewhat broader 
transition (~0.135°). For magnetic fields varying 


8H. K. Onnes and W. Tuyn, Leiden Comm. b160 (1922). 


from 0-84 oersteds dH./dT.=170 oersteds per 
deg.® 
DISCUSSION OF RESULTS 

The experimental results obtained showed that 
a low resistance ratio, high transition tempera- 
ture, and small transition range were always 
found together for a given sample of the metal. 
The fact that the smallest wire always exhibited 
the highest R/Ro ratio just above the transition 
as well as the lowest transition temperature and 
broadest transition is to be attributed, we believe, 
to the fact that the strains set up in drawing the 
wires are probably largest in these wires. The 
lead wires are exceptions to this, but here all 
wires consisted of large crystals which probably 
accounts for the sharpness of the transitions and 
uniform value. 

The effect of small magnetic fields (0-80 
oersteds) was to depress the transition tempera- 
ture, the H, vs. T. curves being essentially straight 
lines. For Sn, Pb, and Ta the observed de- 
pressions are in general agreement with results 
obtained by other investigators. For columbium 
the only data available for comparison are the 
observations of Daunt and Mendelssohn" in the 
temperature range 1.5—4.2°K, where rather high 
fields were necessary. When their H, vs. T, data 
are extrapolated to 9.3°K the value of dH./dT, so 
obtained is approximately the same as that found 
for our samples. 
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corresponding data for the 127- and 254-micron wires 
were not completed. The limited measurements made 


fort a dH./dT. of approximately the same value as that 


or the 25-micron wire. 
10J. G. Daunt and K. Mendelssohn, Proc. Roy. Soc. 


A160, 127 (1937). 
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Fluorescence of Solutions and Dielectrical Properties of Solvents 
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Because of different red shifts of the absorption and fluorescence spectrum of anthracene in 
solution, the longest wave absorption band is separated from the shortest wave fluorescence 
band. This ‘‘(0-0) separation” and the red shift are explained as effects of interaction between 
molecules in the liquid state and are discussed in the light of Onsager’s theory of liquids. 


1 


T has been shown in a previous paper! that 
there is a disparity between the red shift of 
the absorption and fluorescence spectra of sub- 
stances in solutions as compared with these 
spectra in the vapor state. It was found that the 
shift of the fluorescence spectrum is larger than 
that of the absorption spectrum. Consequently 
the longest wave absorption band no longer co- 
incides with the shortest wave fluorescence band, 
as is the case when these bands are measured in 
the vapor state. The discrepancy which will be 
referred to below as the (0-0) separation in- 
creases with the general red shift and can reach 
values up to several hundred cm. The red shift 
is accompanied by a broadening of the bands 
both in absorption and emission. The broadening 
and red shift of bands in solutions are analogous 
to effects produced in spectral lines by the 
pressure of foreign gases. No analogy to the (0-0) 
separation has, however, been observed in line 
spectra. 

Several attempts have been made to explain 
the red shift but these have usually been re- 
stricted to primary correlation with the dielectric 
properties of solvents. An effort by Henri? to 
connect the dielectric constant of the solvent 
with the red shift was not very satisfactory. 
Seshan* attempted to set up a relationship be- 
tween the red shift and the refractive index. 
Discussion has also been directed to the influence 
of the dipole moment on the shift.‘ All these 
explanations have been qualitative. It appears 


1S. Sambursky and G. Wolfsohn, Trans. Faraday Soc. 


36, 427 (1940). 

2V. Henri, J. de bn ing e 3, 181 (1922). 
oP, K. dy Trans araday Soc. 32, 689 (1936). 
*G. Scheibe and W. Frémel, Hand- u. Jahrb. Chem. 


Phys. 9, 174 (1936). 


likely that a quantitative explanation of the shift 
as well as of the correlated (0-0) separation can 
be developed only on the basis of recent theories 
of the liquid state. Such an explanation when 
found would also provide an optical method for 
the study of the liquid state through measure- 
ments of the fluorescence and absorption spectra 
in solutions. 

If we assume that the observed effects are 
caused by the interaction of the fluorescing 
molecule with the molecules of the solvent, we 
shall have three distinct types of interaction to 
consider: the orientation effect, the induction 
effect, and the dispersion effect. The total energy 
arising between two molecules f and s with the 
dipole moments y; and y, and the polarizabilities 
a, and a, is given by® 


12 
a; 


2 2 
Ger 
— 1 


In our case index f will represent the fluorescent 
molecule and the index s the molecule of the 
solvent. The sum in the last term has to be ex- 
tended over all transition frequencies v,, and v., 
of both molecules and their corresponding mo- 
ments of transition a,, and a,,. The ground 
states of the molecules f and s are g and o and 


ex = Eg) /h= — (la) 


The orientation effect given by the first term pre- 
supposes that both molecules have permanent 
dipole moments. The second and third term 
represent the induction effect, i.e., the inter- 
action by polarization of the dipole molecule 


5 F. London, Trans. Faraday Soc. 33, 8 (1937). 
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with its neighboring molecules, while the dis- 
persion effect (last term) describes the syn- 
chronized vibrations induced in one molecule by 
the rotating charges of another. This interaction 
occurs in both polar and non-polar molecules and 
can reach considerable values in cases of reso- 
nance. 

In considering the (0-0) separation between the 
fluorescence and absorption spectrum and the 
total red shift of both, it will be found convenient 
to choose solvents whose molecules possess no 
dipole moment so that y, is zero and there is no 
orientation effect. If we also choose as fluorescent 
substance a dipole-free molecule such as anthra- 
cene, then the interaction with the surrounding 
molecules will be restricted to the dispersion 
effect, as long as the fluorescent molecule is in 
its ground state, and consequently the energy of 
the ground state will be lowered by a certain 
amount L; compared with its value in the 
gaseous state. When, however, this molecule has 
been excited by absorption, two effects occur, as 
follows: 

First a dipole is generally created, polarizing 
the solvent by induction. It should be empha- 
sized that the two accompanying phenomena of 
absorption and induction in solution should be 
viewed as a single process. Therefore the energy 
of the excited anthracene molecule in solution 
will be larger than that of the same molecule in 
the vapor state by the amount of energy in- 
volved in the induction process. This induction 
energy is distributed over the excited molecule 
and the polarized molecules of the surrounding 
solvent and all or part of it is dissipated as 
thermal energy and therefore does not return, 
when the excited molecule returns to its ground 
state. The difference in energy between the ab- 
sorption and emission processes can thus at least 
partially be identified with the dissipated energy. 

There is a second effect which will operate in 
the same direction: the dispersion energy will 
also change by excitation. It will further be 
shown that the binding of the fluorescent mole- 
cule with its surroundings will generally increase 
by excitation, causing an increased lowering of 
the excited state as compared with its value in 
the gaseous state. It can be assumed that a part of 
this increase of energy too will be dissipated as 
thermal energy before fluorescence occurs. 


Let E; and E, be the energies of the ground 
level and the excited level in the gaseous state, 
I, the dispersion energy in the ground state, 
L, the dispersion energy in the excited state, and 
D the induction energy accompanying the crea- 
tion of the dipole » by excitation. Then the 
energy of the ground state of anthracene in the 
surrounding solvent will be given by E’=E£,+1). 
After absorption the molecule reaches an excita- 
tion energy 


—(L2—L1)5, 


where y and 6 are the fractions of D and of 
I2—IL, dissipated before fluorescence occurs. 
Hence the excited state before fluorescence has 


the value 
;=E2+L24+ D. 


The (0-0) separation is given by 
(2a) 


while the total red shift of the fluorescence 
spectrum is 


(E2—E” ;) (E,-—E’) (D+L2—L)). (2b) 


Quantitative results may be obtained with the 
aid of Onsager’s theory of polar liquids which is 
in good agreement with experiment despite the 
simplicity of his assumptions. According to 
Onsager® we may picture the anthracene molecule 
in the center of a spherical cavity with a dielectric 
constant 1 and a radius a. In a sufficiently dilute 
solution the sphere will be surrounded only by 
molecules of the solvent whose dielectric constant 
is «. Let us suppose that the excited anthracene 
molecule in the vapor state has a dipole mo- 
ment yu. In solution it will induce a reaction field 
in the cavity whose field strength is given by 


2(e—1) 
Qet1 a?’ 
(n?+2)(2e+ 1) 
3(2e+n?) 
is the dipole moment of anthracene in the reac- 
tion field. yu’ is of course larger than u» owing to 
the fact that both n, the refractive index of 


anthracene, and e are larger than unity. The 
interaction energy D between the excited mole- 


where 


® L. Onsager, J. Am. Chem. Soc. 58, 1486 (1936). 
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TABLE I. Summary of data. The frequencies 
are given in cm™, 


eo. 
Solvent é A Yabs vst “ton. (fluor.) 
CO: 1.54 0.325 26990 26900 90 415 
ntane 1.83 0.470 26716 26600 116 715 
xane 1.875 0.490 26652 26524 128 791 
cyclohexane 2.01 0.555 26602 26458 144 857 
decalin 2.165; 0.620 26512 26349 163 966 
CCl, 2.22 0.642 26378 26136 242 1179 
p-xylene 2.27 0.660 26410 26155 255 1160 
dioxane 2.28 0.665 26468 26294 174 1021 
benzene 2.28 0.665 26378 26136 242 1179 
mesitylene 2.36 0.695 26419 26143 276 1172 
diphenyl 2.53% 0.757 26273 25933 340 1382 


naphthalene 2.54 0.760 26239 25875 364 1440 
phenanthrene 2.72° 0.823 26000 25580 420 1735 


‘aot when otherwise is indicated, the figures given refer to « 
at ° 

2¢ at 10°. Calculated from the data given by Uhlig, Kirkwood, and 
Keyes, J. Chem. Phys. 1, 155 (1933) and the density measurements of 


Lowry and Erickson, J. Am. Chem. Soc. 49, 2729 (1937). The result is 
in agreement with the former measurements of Linde, Wied. Ann. 


cule and the solvent is then given by D= — Ry’ or 


9(2e+n*)? 


a’ a’ 

It is above assumed that all or a part of this 
energy is dissipated and that the energy of the 
emitted radiation will be smaller by this amount 
than the energy absorbed. If this assumption is 
correct, the (0-0) separation, as far as it is caused 
by the dissipation of induction energy, will be 
proportional to the factor A. 

Since the dispersion effect also has to be con- 
sidered as an interaction of periodic dipole 
moments, this proportionality will also hold for 
the dispersion energies L; and Le, except in case 
of resonance, which will be dealt with later. 
It is to be expected that the (0-0) separation as 
a whole will be proportional to the factor A. 
As D, Ly, L2 are proportional to A, the total red 
shift Eq. (2b) should likewise be proportional to 
this factor. 


2 


Owing to the smallness of the (0-0) separation 
as compared with the width of the bands it can 
be measured accurately only if absorption of 
fluorescence light in the solution is prevented. 
An arrangement similar to that described in the 


previous paper was therefore used, giving a 
fluorescent layer not more than 1 mm in depth, 
with an anthracene concentration of 10-* mole/ 
liter. Wave-lengths were measured for both fluor- 
escence and absorption against a superimposed 
iron spectrum. Photograms were taken with a 
large Moll microphotometer. The solvents used 
and the corresponding values of ¢ and A are given 
in the first columns of Table I. For the calcula- 
tion of A, a value of the refractive index n of 
anthracene equal to é¢ for this substance was 
used (€=2.7). The range of the factor A is 
limited by the comparative small variation of « 
in dipole free solvents. In no suitable solvent 
does the dielectric constant exceed 2.8. Sub- 
stances which are liquid at room temperature 
were investigated at 20°C, whereas diphenyl, 
naphthalene, and phenanthrene were observed 
at temperatures lying roughly 10° above their 
melting point. The absorption and fluorescent 
cell was maintained at the required temperature 
in an electric furnace. A special cell was con- 
structed for the observation of anthracene in 
liquid COs, a cross section of which is given in 
Fig. 1. The quartz plates Q were cemented with 
piceine to two covers C screwed on a cylindrical 
iron ring R. Absorption was observed in the direc- 
tion AB, while for the observation of fluorescence 
a beam of light was directed into the aperture F 
and the fluorescence emitted through B meas- 
ured. A piece of copper T in which a thermo- 
couple could be introduced for temperature 
control was attached to one of the covers. For 
the filling of the cell one cover was unscrewed 
and the necessary quantity of anthracene de- 
posited by evaporation from a solution in acetone. 
A weighed quantity of solid CO. was then 
placed on the plate and the cover screwed back. 
The anthracene is quickly dissolved in liquid CO». 


3 


Quantitative results are given in Table I and 
Figs. 2 and 3. The key to the numerical symbols 
in Fig. 2 and 3 is given in the first column of 
Table I. In columns »,», and vz; the frequencies 
of the longest wave absorption band and the 
shortest wave fluorescence band are given. The 
difference between these figures is the (0-0) 
separation. The last column gives the total red 
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Fic. 1. Cross section of cell for observation of 
anthracene in liquid COs. 


shift, based on the value \=3660A; v=27315 
both for and vy; in the vapor state.” 
As can be seen from Figs. 2 and 3, in which 
the (0-0) separation and the total red shift of 
the fluorescence spectrum are plotted against A, 
the expected proportionality holds for a number 
of dipole free liquids. A number of the solvents 
display a larger separation and red shift than 
would be allowed for by a linear relation. In all 
of these cases the main region of absorption lies 
in the neighborhood of the absorption bands of 
anthracene itself, i.e., between 2500A and 2000A. 
The solvents obeying the linear law are those 
which absorb only in the far ultraviolet, e.g., 
hexane near 1500A and CO, near 1600A and 
1000A. These circumstances suggest that the 
additional effects may be explained as a result of 
resonance. It cannot be expected that Onsager’s 
theory dealing with static dipoles should give a 
full account of the dispersion effect as far as it 
7 This value was adopted on the basis of absorption 
measurements of anthracene vapor by Gépfert, Zeits. 
wiss. Phot. 34, 156 (1935) and measurements by Pring- 
sheim, Ann. Acad. Sci. Techn. Varsovie 5, 29 (1938) in 
the fluorescence spectrum. The two maxima given by 
Pringsheim—3640A and 3690A—are probably imaginary 
and due to reabsorption, the real maximum apparently 
being situated in the gap between them. With A=3660A 


as maximum the distance to the next maximum becomes 
about 1400 cm™ as was to be expected. 


is dependent on frequency. This must be found 
by making a direct use of Eq. (1). 

Since the fluorescing molecule is found in the 
upper state gy’ before fluorescence occurs, the 
difference between the dispersion energy in the 
ground state ¢ and in the excited state ¢’ is due 
to the transition gy’. Owing to this transition 
the last term in Eq. (1) will be changed mainly 
by vey changing in v,-,. Hence the difference 
L.—L, for the dispersion energy is in the first 
place given by 


Since »,, is greater than v,,- in all observed cases, 
there results an increase of the dispersion energy 
in the excited state of anthracene and conse- 
quently an additional separation and red shift. 


Since Oar ‘Ver iS proportional to the molar absorp- 
tion coefficient for the frequencies »,,, one can 
estimate the value of the sum in Eq. (4) from the 
absorption coefficients of the solvents used. Both 
for the (0-0) separation and the red shift a 
proportionality is to be expected between 

2 


and the deviations of the actual measurements 
from the straight lines in Figs. 2 and 3. The 
agreement is fairly good in all cases where data 
for the absorption coefficients are obtainable up 
to at least 2000A. 
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It should be possible to explain the non-linear 
effects directly from Onsager’s theory, if one 
assumes that e varies with the frequency and 
calculates its value from that of the refractive 
index of the solvent at the maximum absorption 
wave-length of anthracene. Unfortunately such a 
calculation is possible only for benzene and even 
here only approximately, while data on the re- 
fractive indices in the ultraviolet are not avail- 
able for the other solvents absorbing in the 
nearer ultraviolet. Values of n for benzene*® have 
been measured up to 2700A, where n=1.65. By 
extrapolation into the main absorption region 
near 2500A, a value for m may be assumed 
lying between 1.7 and 1.8 which would give 
e=n? as roughly equal to 3. This makes the 
factor A in Eq. (3) a little greater than 0.9. 
This is the value actually obtained if the point 
for benzene in Fig. 2 is moved horizontally to 
cut the straight line. The same procedure carried 
out in Fig. 3 gives a correction somewhat too 
large. 

It may also be pointed out that the broadening 
of the bands depends on ¢ and that this broaden- 
ing becomes more pronounced in the cases of 
resonance. With CO, the bands are still narrow 
enough for the secondary maxima to be resolved, 
on the violet side in absorption and on the red 
side in emission, whereas for example with 
naphthalene the maxima are very shallow and 
blurred. The values of the half-breadth for the 
first bands in fluorescence are approximately 
400 for 500 for dioxane, and 
900 cm- for naphthalene. 

It is not to be expected that the interpretation 
of red shift and (0-0) separation which has been 


8 T. M. Lowry and C. B. Allsopp, Proc. Roy. Soc. A133, 
26 (1931). - 
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Fic. 3. Red shift plotted against the factor A in Eq. (3). 


advanced above will provide a satisfactory quan- 
titative fit for all cases observed. Both Onsager’s 
model which replaces the molecules of the solvent 
by a continuum and the simplifications involved 
in Eq. (4) for the dispersion effect are approxima- 
tions which necessarily introduce a certain dis- 
crepancy between experiment and theory. Thus 
in the case of CO: the red shift found is smaller 
and the (0-0) separation greater than was ex- 
pected from theory. In the case of benzene, 
carbon tetrachloride, p-xylene, and in particular 
mesitylene, again, the values found for the 
(0-0) separation were somewhat greater than as 
calculated from Eq. (4). 

It is planned to extend these investigations to 
the case of polar liquids and polar fluorescent 
substances, whereby it is hoped to shed further 
light on the mechanism of interaction between 
fluorescent solutes and their solvents. 

The writers are grateful to Dr. M. Schiffer 
for much helpful discussion and wish also to 
thank the Elizabeth Thompson Science Fund for 
its assistance in providing a wave-length com- 
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The Energy Values of the 3d*4 Electronic Configuration of Cobalt 
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(Received July 6, 1942) 


The energy matrix of electrostatic and magnetic spin-orbit interactions is set up for d*p. 
The secular equation is fitted approximately to the experimental values of the energy levels 
of this configuration in Co I. The calculated energy values are given. The interaction param- 
eters are compared with those of the ‘‘deep” configuration of CoI and Nil. The calculated 
values of the energy are compared with the empirical ones published by Russell, King, and 
Moore, the conclusion being drawn that there is an agreement between them good enough 


to confirm the assignments of these authors. 


HE uneven energy levels of atomic spectra 


have somewhat resisted attempts made to *Ga: 
associate them in electronic terms and to at- 
tribute electronic configurations to them. Among 
other things, this fact is due to the large number "Go: 
of uneven terms existing in the large majority of 
spectra. ‘F; 
Russell, King, and Moore represent in their 
analysis of the arc spectrum of cobalt! values of ?F,: 
the energy levels of the configuration 3d*4p con- 
veniently associated in terms. 
The object of the present work was to attempt ?F;: 
to confirm through calculation the assignments 
made by Russell, King, and Moore in this 
configuration. 
We never expected to resolve thoroughly the 
secular equation of the energy; but we were con- , D.: 
fident that some values of the energy, although . 
without great precision, would contribute help- 
fully to the realization of such a confirmation. ‘Ds: 
This hope was based fundamentally on the good 
results already obtained by Catalan and An- 
tunes? in the case of the deep configurations 2D): 
d™~*s2, d™—'s and d”. 
THE ENERGY MATRIX 2D: 
The elements of the energy matrix of electro- 
static interaction of d*p are the following: 
2D.: 
A. Diagonal Elements (*G Referred to as Zero) 
1 Henry Norris Russell, Robert B. King, and Charlotte *Pa: 
E. Moore, Phys. Rev. 58, 407 (1940). 
2 M. A. Catalan and M. T. Antunes, Zeits. f. Physik 102, 
432 (1936). b: 
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225 
3 
2 
15 
12F2+ 10F,+12F,’ + 
4F,’ 


54 9 
12F2+ 


13 2 216 
81 
10 
F.- 


7 
14F,— 70F,— 


7 
45Fy— 2F'2+ 


12 27 378 


1SF2,— 75Fi— 
27 


22F2+135Fit+ 2Gi+ 21G; 
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SFat+ 45 Fat 4F'2+ 3G; 


4S: 15F,— 75 


2S: 15F,— 75Fi+15F'2 


B. Non-Diagonal Elements 


15(15)# 
= 


6(3)! 
3\3 
2 Fle? 2G G 
(?F|e?/rij|? Fa) = roy “| 1— 3G3} 


2D, |e? 2D G,—12G 


a 2 1 3 


= -{ Gi-— 7G3} 
-{ 4F2 — Gi-— 
= 3(21)* - { —3Gi— G;} 


In the notation used above, that employed by 
Condon and Shortley,’ the value of F; is 


F,=F,(n2,n2); Fy(n2, n2); 
F’.= F,(n2, n’1); 
G,=G,(n2, n’1); G3=G;(n2, n'1). 


The elements of the energy matrix of the spin- 
orbit interaction grouped by successive values of 
the internal quantic number J are given below: 


7 
= (14)8 —F,+5Fi— —F'2 J=5} 
3 2H 4G 
= { i +9G; ele 
J=4} 
*H ‘G °Ga ‘F 
3 2 11\! 11)! 
——b ———a -(=)« 0 
5 (15)! 15 5 
2 1 (11)! | (11)! 
11\ (11)! 1 1 
| —(—) ———(3a+20)| —-(6a+8) —-(a—2b 
1/11\! 1 1 15)! 
5 4\15 2(15)! 10 4 
7 (11)! (15)4 1 
0 (a+b) —-(a—26) —-—a | —-(l1la—d) 
8 4 8 


3E. U. Condon and G. H. Shortley, The Theory of Atomic Spectra (1935), p. 177. 


| (9+?) — | | (q+0) 0 0 0 
I 672 I 
(q+0) — | » (<)- (q+0)— 0 0 0 
if 4(Z) 4\9 4(Z) 
6 (17) L)8 (¢)8 8 L\8 
4(Z) if I ¢ I 
1Z (PLE 4(1Z) Z L 
D—— D--— — q- 0 D—— 0 ($)— 
62 I I I ¢ ¢ 
Zz 4\9 ¢ (SE) if 
(q+). (q+?) » 0 (q+?) (<)- 
8 8 8 8 4(ST)8 £\8 
0 0 0 q— D q-— 
(Se e(ST) I 
x Z)8 
L\8 l €\9S L\Z ¢ 
0 0 (97 (—)- D4 » (q+?) » (97+0¢) (q—¢)- 
"dy "d+ "Az We "dz dy Dy 
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(q—»)-— 0 (q+) - 0 0 0 0 0 
if I (Le 
#(9)Z 9 e(L-S-€)Z #(0£)9 +(0¢) 
(qz7—?) (q+07)-— 0 (q7+2) 0 0 0 0 
I I l l 
p— -— q- — — 0 — 0 
6 #(0£)6 6 L 6 
0 0 (q—?g)-— |(@—07) | -}-| (qz—»)- 0 
(q+?) (97+?) p—— —|p———— | | (q+2) 0 
(q+) »——-—_ |(q—07)—--| (q+) | (97-0?) |o—-!| (q+0)—— 0 
0 0 p—— (q—27)—— (97-0) | (9+0zz)— | | 
0 0 0 — — D—— 0 q- D—— 
#\¢ #(9) I I if 
6 (0176 (S16 9¢ (L)6 
0 0 D (¢+0)——— | (q+0)———._ |(97+ 911) D |(¢+0)—— 
0 0 0 0 0 0 (q7—») D D—— | (q+0)—— | (q—2¢)- 
dy "dz "dz "dz "dy We "di ds 
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"di 


‘dr 
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"dz 


"dt 


"ds 


"ds 


"di 


"dy 


4(Z) #(9) 9 
0 o—— | o— | (Qz-—2)—— | (+?) 0 0 0 0 0 0 
+(Z) 
l I € #\9 if 
0 q- D— — — — 0 0 D—— D-— 0 
I I (PT) I 
p—— |v—-—| (9+07)—— | | 0 | 0 
I I 4(1Z) L if 
(q+0) p———| | (q—0)- — 0 (q+0)———_ | (+2) 0 
| 4(S) if if tI 1 
| 9(ST) ral (O1Z)Z (SP 
0 p——— | | (q—07)—-— 0 | (97+) 0 
4(S) if if l if 
0 q- Di — — p--—- — 1) 
0 0 0 0 0 (q—0g)- (q—27)———— 97) —| (@z—»)- 
#\9 I if +1 
#(Z)9 (olz¢ (O1Z)Z| (ST 6 I 
1 tI L 1 if 
0 p——— | (q+0)——_ — ome (q- (q+?) (q+?) 
I I I (1) I (012) 
0 0 0 0 0 0 (q+?) (q+2)——- (q—?I1)- 
8 I I #(01Z) 
S> ‘ds "di dy "di dz "dz “dy dy 


= 
| 
| | | | | | 


ELECTRONIC CONFIGURATION OF COBALT 367 


J=} 
‘D, ‘4D, 4p 2p, 2P, 2P, 25 
(230-458) b) (a+b) oy | 0 
a| = +5) | ~(sa—48) +s) | 0 
Ds (a ) (a—2b) 
(a+b) (a+25)| — (0-25 
12(3)) (6(42)) 12" 6(2)! ( | 242)! 
(2)! 2(7)! 2 2 10 
*P, —a —a —b 0 - 
3 3 (6)! (3)! 3(6)! 
5(14)! 1 7 1 (21)! 
*P. -——6 0 —b - 
12 6 2(42)! 2(3)3 2 3(2)4 
1 10 (21)! 
0 0 —-(a—2b) |———(2a—5b) |-——a | -——a 0 
3 3(2)4 3(6)# 3(2)! 


In these matrices, a and Db represent, respec- 
tively, the integrals of the spin-orbit interaction 
in the configuration d*p of an electron d and an 
electron p. 


THE VALUES OF THE ENERGY 


The elements of the complete energy matrix 
are the sum of the electrostatic elements together 
with the corresponding elements of the spin- 
orbit interaction (we disregard the interaction 
between configurations). 

The secular equation, obtained by subtracting 
the unknown energy W from the diagonal ele- 
ments of the complete matrix and equating to 
zero the corresponding determinant, has been 
treated as follows: all the elements depending 
exclusively on the energy of the spin-orbit inter- 
action have been replaced by zero. The secular 
equation, thus simplified, has been divided into 
a certain number of equations, the majority of 
them linear for W. The empirical values given 
by Russell, King, and Moore, were replaced in 
the linear equations referring to ‘Gs, as zero, and 
the values of all the parameters F2, Fy, ---, a, 6, 
with the exception of G; were determined by the 
method of least squares; the latter (G:) was ob- 


TaBLe I. Calculated values and empirical values of the 
energy of cobalt in the configuration 3d*4p. 


Observed term Calculated Observedterm Calculated 
‘Gs, 32431 32431 45971 44534 
‘Gy 32465 32919 ‘Dry 46330 44734 
{Gy 333173 33318 ‘D>, 46260 44877 
33674 33629 46502 44964 
‘Fy 32842 32449 44480 44648 
‘Fy, 33467 32963 44658 44774 
‘Fy 33946 33363 ‘P, 44858 44850 


‘Fy 34196 33649 
‘De, 32027 33501 


‘S 46563 47291 


‘De, 32654 34101 "Hy 50703 50846 
‘De, 33151 34529 “Hy 50903 50760 
‘Dy 33449 «34785 

33440 34176 43922 45494 
2G2,, 34134 35067 43911 45496 
53276 52502 Dey 45688 47453 
2G, 53374 52480 Dey, 46455 47840 
35451 34919 2Pe,, 43538 45585 
2Fe,, 36330 35730 43130 45330 
47225 45473 48334 48660 
47129 45378 2P, 48160 49157 
2Fey, 50579 «51217 78755 
50712 51301 2Pey 78663 
36092 35742 47978 47291 


36875 36496 
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TaBLe II. Interaction parameters of the configuration 
3d*4p of cobalt. 


M. T. ANTUNES 


TABLE III. Interaction parameters of the deep configura- 
tions of Co I and Ni I and of 3d%4p of Co I. 


107.46 
147.63 
199.89 
14.385 


tained by means of non-linear equations. The 
parameters F», ---, a, b, entering in all the equa- 
tions, were substituted by the numerical values 
obtained. Finally the equations were resolved in 
order of the energy W. Thus were obtained the 
values of the energy (referring to ‘Gs,;=32431 
cm) which are found in Table I together with 
the empirical values given by Russell, King, and 
Moore. 

The values of the parameters are shown in 
Table II. 

It would certainly be of interest to resolve more 
thoroughly the secular equation. In any case, the 
values of the energy which we have just pre- 
sented already confirm the assignments made by 
Russell, King, and Moore. Note that the differ- 
ence between the calculated values and the em- 
pirical ones attain, at the most, 2000 cm-! (?P,) 
which already shows a notable restriction of those 
regions in which may be found each term. Fur- 
ther, there are already calculated values for the 
energy which are closer to the respective em- 
pirical values than to others. The values of the 


3d74s2 Col 3d84s 3d* Col 3d%4s? Nil 3d%4p Co! 


Fi(n2,n2) 1522 1225 1653 1380 
F,(n2, n2) 112.3 110.9 124.0 107.5 
a —520 — 465 —493 —655.8 381.08 


@ The lack of accordance concerning the sign is explained by the fact 
of Marvin's having represented by a the integral of an electron d with 
the opposite sign, while we represent by a the integral itself. 


parameters are more similar to those of other 
configurations, as can be seen in Table III; the 
parameters of 3d74s? of Co I, 3d%4s of Co I, 3d* of 
Col, and 3d*4s? of Nil are those cited by 
Marvin.‘ 

It should be noted in passing that the values 
of some of the parameters may, because of the 
degree of approximation of our calculations, be 
varied within relatively large limits without pro- 
ducing excessively large variations of energy. 
Finally, it should be noted that a calculation of 
this kind necessarily shows up those large devia- 
tions between terms which, because unnoted, 
often deceive empirical classifiers, leading them 
into error; it is evident, then, why Russell, King, 
and Moore did not find the term *P.; the present 
calculation shows that this term is quite high. 

I wish to express my sincere gratitude to Pro- 
fessor Henry Norris Russell for his kindly interest 
in this work. I thank the Instituto Para a Alta 
Cultura of Portugal for the financial aid it has 
given me. 


*H. H. Marvin, Phys. Rev. 47, 521 (1935). 
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A point-to-hemisphere discharge was used for measuring free electron mobility coefficients 
in mixtures of hydrogen and nitrogen, and in pure nitrogen. Consistent results required the 
use of at least a trace of hydrogen or a preliminary clean-up of the tube with hydrogen. 


DIFFICULTY encountered in measuring 

the free electron mobility coefficient in ni- 

trogen at high current densities has been the 

tendency of the discharge wire to become too hot 
and melt.! 

In some measurements of the angular distribu- 
tion of current from a point,’ it was found that 
the use of a point-to-hemisphere discharge gives 
a fairly good hemispherical distribution of cur- 
rent, and that changes in the extent of collecting 
surface produced only minor changes in the po- 
tential required for a given current while redis- 
tributing the collected distribution considerably. 

For this reason, a tube of the form shown in 
Fig. 1 was used in which there is a probe which 
can be advanced along the axis of symmetry and 
the current to which can be measured separately 
from the rest of the collector. As the probe is 
advanced toward the discharging wire end, the 
solid angle subtended by the probe at the wire 
end increases, and the fraction of the current the 
probe must collect to avoid disturbing the angu- 
lar distribution of current increases accordingly. 

The tube was used in the circuit shown in Fig. 
2. Potential from a rectifier P was adjusted to 
draw a total discharge current measured on J 
(usually 0.3 ma) and the potential required was 
measured with the multiplier R and meter V. 
The resistance R, was adjusted until the probe 
current 7 was in the same proportion to the total 
current J as the solid angle subtended by the 
probe was to the total half-sphere of the collector. 
The potential of the probe due to the potential 
drop through R, of i was then measured by using 
a separate rectifier p, multiplier r, and meter », 
by varying the potential from p until an RCA 
991 neon tube N ceased discharging and a 0-110 


(1942) H. Bennett and L. H. Thomas, Phys. Rev. 62, 41 
?W. H. Bennett, Phys. Rev. 61, 53 (1942). 


volt electrostatic voltmeter EV read a minimum, 
usually less than the lowest scale division. The 
neon tube was necessary to avoid injury to the 
electrostatic voltmeter while adjusting the po- 
tential from p. 

The method of measuring probe potential 
seems to be quite accurate because changes of i 
by a factor of four either way produced changes 
in observed probe potential which changed the 
calculated value of mobility coefficient K by only 
a few percent. 

Although it was not expected that the poten- 
tial distribution would follow very closely the 
form for the space change limited free electron 
current between concentric spheres for twice the 
value of the current to the hemispherical col- 
lector, the experimental observations showed that 
the use of the probe on the axis in the above 
manner gave remarkably good consistency and 
good agreement with the results obtained with 
the cylindric tube. 

If the ionization sheath with the inner sphere 
negative has a radius 7o, the field at the surface 
of the sheath is Eo, the potential drop through the 
sheath is Vo, and the electrons have a velocity 
u= KE}, the current J» to the outer sphere with 
radius r is related to the potential V by 


r 
V= v+(—) (—) —1.145 
4K R 


4K 


where 


R=n|1- 


From this, the potential difference Vj; be- 
tween the probe at radius r; and the hemisphere 
with radius rz in terms of hemisphere current 
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I= is approximately for r2>r1>ro 


3I\! 
(=) 


or 


which does not contain"R. 

In Table I is shown a typical set of data. These 
are for pure hydrogen. The consistency of these 
data as well as all of that for the various mixtures 


Fic. 1. Experimental tube, point-to-hemisphere, with 
probe, for use at pressures of 1 atmosphere and more. 


= 


= 


Fic. 2. Circuit used for mobility measurements with probe. 


is noteworthy as showing that this method is 
valid even with the probe only a few millimeters 
from the discharging point. 

In Table II are given the values for K for 
various mixtures of hydrogen and nitrogen, and 
it is noted that the initial increase and subsequent 
decrease in K noted previously! is confirmed. 

The measurements of K for pure nitrogen were 
consistent only if the tube were first filled with 
hydrogen and conditioned, before filling with 
nitrogen. Filling with nitrogen after the tube had 
been open to the air resulted in a refusal of the 
nitrogen-filled tube to condition, and the data 
were erratic. 

A complete duplicate set of measurements for 
each of the gases and mixtures was made at 
3-atmospheres pressure. These results were con- 
sistent and gave values for the mobility coeff- 
cient (K;3) related to that at 1 atmosphere (K,) by 


K,=v3-K3. 


In some early experiments, a discharge was 
drawn from a steel point in pure nitrogen and it 
was noticed that a black deposit formed on the 
point. When the tube was opened and air was 


Taste I. Experimental data. 


n 1/8 K 
0.369 0.687 960 7.5X 10° e.s.u 
0.495 ‘0.614 840 7.5 
0.621 0.550 770 7.4 
0.878 0.447 590 7.9 
1.515 0.255 330 8.3 
(2.78 =r2) Average 7.6X 105 e.s.u 


Ia ® eo = DBA 


ev 
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allowed to get to the point, there was a sharp 
explosion and the black deposit turned rust- 
colored. This suggests that a nitride of iron was 
formed which is unstable in air. 

In some later experiments, it was found that 
drawing a current of about 0.6 to 2.0 milliampere 
in mixtures of hydrogen and nitrogen produced 
ammonia which gave a strong odor if the tube 
was opened, and this ammonia behaved as an 
electron-attaching impurity deconditioning the 
tube. After ammonia had been formed in such a 
mixture, reduction of current to 0.3 milliampere 
decomposes the ammonia and conditions the tube 
again. Opening the tube, no odor of ammonia 
could be detected. 


TABLE IT. Composition of gases and mobility of electrons. 


Composition K 

100% He 7.6X 105 e.s.u. 
30% H2+70% No 7.8 
10% H:+90% Ne 6.9 
3% H2+97% Ne 6.8 
1% H2+99% 6.2 
0.3% H2+99.7% Ne 5.7 
0.1% H2+99.9% Ne 5.0 
> Ne 3.9 


In mixtures containing 1 percent or less of 
hydrogen, the wire end becomes white hot at 
currents which produce ammonia and it may be 
that the formation of ammonia is thermal. 
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Measurements of the Hall constant and conductivity of cuprous oxide show that the 
exponential law of temperature dependence is not obeyed. It is shown that the departure from 
this law is caused by a loss of conduction holes with time, and an anomalous decrease in the 
mean free path in the vicinity of 100°C. An experiment is described which indicates that the 
rate of aging at 100°C is increased by the application of an electric field. 


INTRODUCTION 


HE purpose of this paper is to give data on 

the Hall effect and conductivity of cuprous 

oxide which contribute to the understanding of 

the mechanism of electrical conduction in this 
substance. 

It is generally believed that cuprous oxide 
containing a stoichiometric excess of oxygen is a 
semiconductor. Semiconductors are substances 
which possess a conductivity which is inter- 
mediate between that of insulators and metals. 
Furthermore, their conductivity approaches zero 
as the temperature approaches absolute zero. 
The published literature on semiconductors is 
extensive. Good accounts of their properties 
are given by F. Seitz! and N. F. Mott and R. W. 
Gurney. 


1F, Seitz, Modern Theory of Solids (McGraw-Hill Book 


Company, Inc., New York, 1940), pp. 62 ff. and 186 ff. 
?N. F. Mott and R. W. Gurney, Electronic Processes in 
Tonic Crystals (Oxford, 1940), pp. 153 ff. 


The application of a magnetic field to a sub- 
stance carrying an electronic current causes 
several electrical and thermal phenomena,* 
of which the most important is the isothermal 
Hall effect. The Hall effect is a rearrangement 


Fic. 1. Schematic diagram showing the electrical con- 
nections used in the measurements of Hall effect and con- 
ductivity in cuprous oxide. Ey, Hall voltage. Ec, conduc- 
tivity voltage. H, magnetic field. J, sample current. 


3 For a review see A. Sommerfeld and N. H. Frank, Rev. 
Mod. Phys. 3, 1 (1931). 
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TYPE K 
POTENTIOMETER 


ov 


4sv's" 


LS VARIABLE 
| RESISTANCE 
BOX 


SAMPLE HOLDER 
Fic. 2. Circuit diagram of the apparatus. 


of the equipotential surfaces in a solid carrying a 
current, caused by the application of an external 
magnetic field perpendicular to the current flow. 
Thus, if the conductor be a flat strip, such as 
that shown in Fig. 1, in which the current flow 
is along the length, the application of a magnetic 
field perpendicular to the strip will cause the 
equipotential surfaces to shift from positions 
perpendicular to the axis of the strip to those 
making a slight angle with the original positions. 
As a result of this shift, two points situated at 
opposite edges of the strip which are initially at 
the same potential will be at different potentials 
after the application of the magnetic field. This 
difference of potential is known as the Hall e.m.f. 
The magnitude of the shift is proportional to 
the magnetic field strength and the current 
density. The constant of proportionality is 
known as the Hall constant. A simple theory! 
of the Hall effect shows that the Hall constant R 
is given for a semiconductor by the equation 


where n is the number of electrons per unit 
volume and ¢ is the electronic charge. 

In cuprous oxide the sign of R is positive, 
meaning that the current carriers are positively 
charged. This anomalous result was first ob- 


tained by W. Vogt.‘ A satisfactory explanation 
is given by the band theory of solids.! In cuprous 
oxide small regions of electron deficiency are 
mobile, and play an important role in the con- 
duction. A region of electron deficiency behaves 
like an electron with a positive sign. 

In a study of electrical conduction in a solid, 
the Hall effect is an important quantity. From 
the equation above it can be seen that a knowl- 
edge of the isothermal Hall coefficient yields the 
number of conduction electrons per unit volume 
and their sign. The conductivity is given by the 
relation 

o=Constant X nl, 


where / is the mean free path. Knowing 2, o, 
and the constant, one can calculate the mean free 
path /, Finally, it can be shown that the mobility 
of the conduction electrons is given by the 
equation 


Notable researches on the Hall effect in cu- 
prous oxide have been carried out by E. Engel- 
hard,> W. Vogt,4 and H. Diinwald and C. 
Wagner.* Many other references may be found 
in bibliographies in these papers. 

4W. Vogt, Ann. d. Physik 7, 183 (1930). ; 

5 E. Engelhard, Ann. d. Physik 17, 501 (1933). 


6 H. Diinwald and C. Wagner, Zeits. f. physik. Chemie 
B22, 212 (1933). 


f 


_ 

| 
| 
| 
a 
U 
d= Reo. k 
y 
3m 1 
Pp 
lc 


nemie 


HALL EFFECT AND CONDUCTIVITY OF CUPROUS OXIDE 373 


APPARATUS 


To measure the Hall effect in cuprous oxide it 
is necessary to have a device which will measure 
voltage of the order of 10 millivolts without 
disturbing the electrical system. A very satis- 
factory instrument is the Leeds and Northrup 
Type K potentiometer. The conductivity was 
determined by measuring the current through 
the sample and the corresponding voltage drop 
between two points 1.25 inch apart in the line 
of current flow. The temperature of the sample 
was determined by a copper-constantan thermo- 
couple fastened to the cuprous oxide with Scotch 
tape. All of these quantities involve voltage 
measurements, and they were made with the 
same potentiometer and a switching system. 
The circuit diagram is given in Fig. 2. 

With the apparatus used, it was possible to 
measure the Hall effect and conductivity in the 
temperature range from —40°C to 100°C with 
continuous variation. The range of temperature 
between 25°C and 100°C was obtained by passing 
dried air through an electric heater, and thence 
through a box containing the sample. The lower 
temperatures were reached by bubbling dry 
air through liquid nitrogen. Various temperatures 
were obtained by varying the rate of air flow 
through the refrigerant. 

The magnetic field employed was 8850 gauss 
from an electromagnet. The field was calibrated 
and shown to be uniform over the area used, by 
means of a Hartmann and Braun bismuth spiral. 

To check the performance of the apparatus, 
use was made of the relation 


E,=RHI/d, 


where E, is the measured Hall effect voltage, J 
is the total current in the sample, and d is the 
thickness of the sample. First, J was varied 
keeping H and the temperature constant. This 
yielded a linear relation between E, and J. 
Second, with fixed J and temperature, H was 
varied. The result was a linear dependence of 
E, on H. Ohm's law was checked at room tem- 
perature for field strengths of the order of 0.5 
volt per cm. As a final check the conductivity 
of titanium dioxide was measured as a function 
of temperature from —40°C to 100°C. The 
logarithm of the conductivity plotted against the 


reciprocal of the absolute temperature yielded a 
straight line. 

In practice it is impossible to evaporate the 
Hall contacts on the sample on an equipotential 
line. To determine the true Hall effect voltage, 
the following procedure was carried out: (a) 
The Hall effect was measured with the magnetic 
field in a certain direction. This yielded V; in the 
equation 

Vi=Ent+ Veter, 


where E, is the true Hall voltage, Ve is the 
voltage drop due to improper electrode align- 
ment, and er is the contact e.m.f. (b) The meas- 
urement was repeated with the magnetic field 
in the opposite direction. This gave V2 in the 
expression 

Ve= —Ent+Verter. 


E, was calculated using the relation 
E,=(Vi-— V2)/2. 


It can be seen from the above equations that 
the contact e.m.f.’s do not affect the Hall effect 
determination. The voltage measurement for the 
calculation of the conductivity is affected by 
contact e.m.f., and these measurements were 
corrected by direct determination of the con- 
tact e.m.f. In no case did this correction amount 
to more than 2 percent. 

The accuracy of the data presented in this 
paper is 2 percent or better. 


PREPARATION OF THE CUPROUS 
OXIDE SAMPLES 


The cuprous oxide used in this research was 
prepared by complete oxidation of Chilean 


TABLE I.* Analysis of copper strip. 


99.968% 
0.0023 
0.0017 
0.0003 
0.0001 
0.00005 


* This analysis was made by the American Brass Company and was 
communicated to the author by C. C. Hein of the Westinghouse Re- 
search Laboratories, East Pittsburgh, Pennsylvania. 
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-010 DIA.GOLD SPOTS FOR 


OxIDE CONDUCTIVITY VOLTAGE 
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DURING OXIDATION 


DETAIL = GRAPHITE LAYER 
OXIDE (AQUADAG) 


SANO BLASTED SURFACE SPRAYED METAL 


Fic. 3. Structure of the electrical contacts to a typical 
cuprous oxide sample. 


strips 2}”’ by 3” by 7”. An analysis of a typical 
batch of Chilean copper is given in Table I. 
Complete oxidation is accomplished by sus- 
pending the clean copper blanks in a furnace at 
1000°C in air for 16 hours. After oxidation, the 
samples are removed from the oxidizing furnace 


Fics. 4A-F. Graphs showing the conductivity and Hall 
constant as a function of reciprocal absolute temperature 
for cuprous oxide quenched from different annealing tem- 
peratures. The discontinuity in the R and o curves of Fig. 
4F was probably caused by the closing of cracks in the sam- 
ple. The ends of the specimen were fixed, and expansion of 
the CuO at 100°C was sufficient to account for the effect. 
The experiment was repeatable, but the data at high tem- 
peratures were not reproducible. 
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and are placed immediately in an annealing 
furnace at a lower temperature. Annealing takes 
place in air for 2 hours, after which the samples 
are quenched immediately in water at 25°C. 
The temperature of the annealing furnace has an 
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' of reciprocal absolute temperature for cuprous oxide if 
quenched from 250°C and aged completely. 
important influence on the Hall effect and con- . 
se ductivity of the cuprous oxide. This temperature reduced copper. (b) The end contacts were a 
are is specified on all curves in this paper. formed by painting with Aquadag and spraying { 
Electrodes were prepared for the sample in the with hot solder, the rest of the sample being jl 
aling following manner: protected with draftsman’s Scotch tape. (c) The q 
takes (a) The oxide was lightly sandblasted to Hall conductivity electrodes were formed by i 
nples remove cupric oxide scales. It was then dipped evaporating gold spots on the sample in an 
25°C. in hot sulfuric acid to remove a layer of cuprous evacuated jar. The structure of the electrodes is 
as an oxide, then immersed in nitric acid to remove — shown in Fig. 3. } 
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Fic. 6. The mean free path as a function of reciprocal absolute temperature 
calculated from the data given by Fig. 5. The dotted line is the true mean free path, 
and the heavy line is calculated assuming a straight line for the conductivity curve 


of Fig. 5. 


RESULTS AND INTERPRETATIONS 


The conductivity and the Hall constant of 
cuprous oxide, treated in the manner previously 
described, does not obey the exponential law of 
temperature dependence. Figures 4 A, B, C, D, 
E, and F show the actual temperature dependence 
in the range from —40°C to 100°C for samples 
quenched from different temperatures. The 
exponential law is obeyed from —40°C to about 
35°C. The kinks in the curves in this region have 
been adequately discussed by Engelhard® and 
Mott and Gurney.” 

The bends in the R and o curves in the high 
temperature region seem to be caused by two 
factors; namely, a loss of current carriers with 
time, and an anomalous decrease of the mean 
free path with increasing temperature. 

Inspection of the R curves reveals that the 
Hall constant is larger at 100°C than the simple 
theory predicts. This indicates that in some way 
conduction ‘‘holes’’ have been removed from 
their role as current carriers. 

That loss of holes or ‘‘aging”’ is not sufficient 


‘to explain the departure from linearity is evident 


from Fig. 5. The data for this curve were taken 
on a sample of cuprous oxide which was com- 
pletely aged; that is, the loss of holes was reduced 
to a negligible number per unit time. The 
linearity of the R curve gives evidence that 
aging has ceased. However, the o curve shows a 
distinct bend. Figure 6, which is a plot of the 


mean free path as a function of temperature’ 
shows clearly that the actual mean free path 
(dotted line) falls considerably below the mean 
free path which would be expected if the o curve 
were linear. 

An interesting experiment was carried out 
which gives a deeper insight into the nature of 
aging. Cuprous oxide quenched from 250°C 
was aged in an oven at 135°C for 5 days. A 
sample treated in this manner ages only very 
slightly after 12 hours; hence, in 5 days equi- 
librium is certainly reached. The sample was 
removed from the oven and stored at room 
temperature for 6 months. The values of R and 
o for this sample after this time are given by the 
initial points in Fig. 7. Aging in an oven at 100°C 
for one hour produced a further drop in con- 
ductivity, which ceased in another hour. During 
the one-hour interval after the fifth hour, an 
electric field of 30 volts per centimeter was ap- 
plied to the sample while it was at 100°C. 
A large decrease in conductivity and a large 
increase in the Hall constant were the results. 
Further aging changed these values very little. 

The following picture of the mechanism of 
aging is proposed: All of the samples of cuprous 
oxide contain cracks which are too small to be 
visible in the samples quenched from low tem- 
peratures (samples quenched from 750°C show 
visible cracks). By diffusion, excess oxygen in 
the lattice escapes, leaving small regions of 
cuprous oxide having a low conductivity. The 
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Fic. 7. The Hall effect and conductivity as a function of time for cuprous oxide quenched 
from 250°C. The dotted curves show aging periods in which an electric field of 30 volts per 


centimeter was applied. 


net effect is to lower the conductivity of the 
whole sample and increase the Hall constant. 
Such depleted regions could be replenished by 
diffusion in an interval of time of the order of 6 
months. Further aging at 100°C rapidly depletes 
them again. The application of an electric field 
sweeps the current carriers out of the lattice more 
rapidly than ordinary diffusion processes. 


In conclusion I wish to thank Professor F. 
Seitz of the University of Pennsylvania for 
guidance and helpful criticism during the work, 
WPA project 25792 for assistance, Mr. C. C. 
Hein of the Westinghouse Research Laboratories 
for cuprous oxide samples and valuable sugges- 
tions, and the Westinghouse Electric and Manu- 
facturing Company for financial assistance. 
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Ionization of Gases by Collisions of Their Own Accelerated Molecules 


H. W. BERRY 
Department of Physics, Washington University, St. Louis, Missouri 
(Received May 25, 1942) 


The ionization of argon, nitrogen, helium, and hydrogen by fast neutrals in their own gases 
has been found to occur in the 1000- to 8000-electron-volt range. The efficiencies of ionization 
are represented in terms of an approximate cross section. The values of the cross sections at 
5000 electron volts are, respectively, for A, Ne, He, and He, 1.5, 0.9, 0.2, and 0.05 cm*/cm! at 
1-mm pressure. In this range, A shows a continuously decreasing efficiency with energy of the 
particle, while H2 and He both show increasing functions. Nz shows both a maximum and a 
minimum. The apparatus for production of the neutral beam and measurement of the ionization 
is described, as are possible defects in the measurements. 


INTRODUCTION 


HE ionization of gases by their own neutral 

atoms has been definitely established for 
the noble gases. Previous work! has shown the 
onset of ionization to occur when the impinging 
atom has an energy of only 30 to 70 ev. This 
energy, just sufficient to detach an electron from 
another atom, is found in each case to be 
roughly three times the ionization potential of 
the particular gas, but for other than noble 
gases ionization by neutrals is not known to 
occur, nor has the ionization function of the 
noble gases been measured above 1000 ev. Since 
the earlier work on other than noble gases was 
confined to less than 1000-ev energy of the im- 
pinging particle, there is no assurance that 
ionization will not occur in an energy range 
higher than this but still below the range charac- 
teristic of the alpha-particle. Experimental 
information in this range is necessary to an 
understanding of the ionizing process and its 
connection with that of such high energy particles 
as the alpha. 

A preliminary report’ of this work gives a 
qualitative description of the ionization of 
argon and nitrogen by their own neutral mole- 
cules in the 1000- to 8000-ev range. Since then 
these measurements have been refined and 
extended to more gases. 


). 
5H. W. Berry, R. N. Varney, and S. Newberry, Phys. 
Rev. 61, 63 (1942). 


EXPERIMENTAL PROCEDURE 


Since the discovery by Rosen and Kallmann‘ 
of the high efficiency of charge exchange or 
neutralization of positive ions as they pass 
through their own gas, the formation of high 
velocity molecular beams has been possible. 
This charge exchange occurs at cross sections 
many times the kinetic theory values for the 
gas molecules themselves; hence a_ classical 
collision is not necessary for the process to take 
place. There need be then no deviation of the 
positive ion or reduction in its velocity in order 
that the neutralization occur. 

To form such a neutral beam which would have 
energies up to 8000 ev it is necessary first to form 
a beam of positive ions with the desired energy. 
These ions are formed in the ion gun (see Fig. 1) 


VELOGITY 
SELECTOR 


NEUTRALIZING 


1ON GUN CHAMBER 


Fic. 1. Diagram of apparatus. 


6H. Kallmann and B. Rosen, Zeits. f. Physik 61, 61 
(1930). 
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by electron bombardment. The electrons emitted 
by a spiral filament are drawn to a cylindrical 
electrode by a potential of about 300 volts. 
After passing through the cylinder they are 
retarded by a reverse voltage sufficient to stop 
and reverse their direction. The positive ions 
formed from the collisions of the electrons and 
gas atoms are drawn by this reverse field to the 
deep slit. The filament is coiled coaxially with 
the gun so that positive ions formed between it 
and the cylinder will be drawn through the 
filament thus neutralizing the negative space 
charge and effectively increasing the emission. 

The ions are accelerated by a voltage between 
the deep slit and the slit in the end plate of the 
electrostatic “‘velocity’’ selector’ which homogen- 
izes the beam. From the velocity selector, the 
ion beam enters the neutralizing chamber 
through a second deep slit. This chamber is 
about 10 cm long and is kept field free by two 
parallel plates at the same potential as the deep 
slit. Positive ions remaining in the beam at the 
end of the neutralizing chamber are removed by 
a reverse field between the parallel plates and the 
plate at the end of the chamber. This plate is at 
a potential 90 volts more positive than the cyl- 
inder from which the ions start, and hence is 
capable of removing all the ions from the beam. 
Those which are neutralized go through a hole 
in this and a second plate which is 45 volts 
above the potential of the cylinder, thus stopping 
low voltage secondary electrons removed from 
the edges of the holes by the beam. To remove 
any secondaries formed in the gas in the neutral- 
izing chamber and hence having the energy of 
the retarding field for the positive ions, a mag- 
netic field is set up across the region between 
the plates. 

Finally the beam passes through a_ space 
charge type of positive ion detector* where the 
ionizing collisions can be measured. Usually 
the detector is composed of two identical cyl- 
inders with tungsten cathodes placed in opposite 
arms of a Wheatstone bridge and operates so 
that the fluctuations common to both cancel out. 


as 2. L. Hughes and V. Rojansky, Phys. Rev. 34, 284 

). 

a b= N. Varney, Phys. Rev. 47, 483 (1935) and 53, 732 
8). 
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Fic. 2. Electrical circuit. 


It was found in an earlier study® of such a 
detector with tungsten filament cathodes that 
the troublesome fluctuations observed were 
caused by the presence of residual gas chiefly 
oxygen and were entirely random in character. 
Further, the use of an oxide coated cathode 
instead of the tungsten wire proved satisfactory.'° 
Hence a single cylinder with an oxide coated 
indirectly heated cathode is used in the present 
experiments, and a fixed resistor substituted for 
the other cylinder. After passing through the 
space charge detector, the beam strikes a nickel 
plate thereby knocking off secondary electrons. 
These secondaries are collected on the space 
charge cylinder by a potential of 10 to 45 volts 
and the resulting current, read on a galva- 
nometer, provides a measure of the neutral atom 
beam intensity. 

The gases used in this experiment are the 
commercial tank product and it is necessary to 
purify them partially. The gas is first passed 
over copper at 450°C for the removal of O» and 
then over phosphorus pentoxide and potassium 
hydroxide for the removal of H:O and COs. 
In the case of He, in addition to the above, the 
gas passes through charcoal at liquid-air tem- 
perature. For He, the charcoal is at solid CO, 
temperature. The gas enters the system through 
capillary leaks. A schematic diagram of the 
apparatus and the electrical circuit is shown in 
Fig. 2. 


RESULTS 


The ratio of the two galvanometer readings is 
used as a measure of the ionization function. 
The galvanometer in the bridge circuit, since 


( om W. Berry and R. N. Varney, Phys. Rev. 57, 1063A 
1 

1H. Karr and R. N. Varney, Phys. Rev. 57, 1064A 
(1940). 


| 
E 62 
ann® 
pass 
high 
sible. 
tions 
sical 
take 
the 
der 
have 
orm i 
rgy. 
g. 1) q 
{ 
] 


380 H. W. BERRY 


3 
x 
Nico 
30 40 50 60 70 80 


Vv VOLTS 


Fic. 3. Ionization of argon neutrals in argon. 


it measures the unbalance of the bridge, is an 
indication of the change in the space charge 
limited current caused by the presence of posi- 
tive ions. The other galvanometer reads directly 
the secondary electron current produced by the 
neutral beam on striking the plate at the end of 
the beam path, or indirectly the size of the 
neutral beam producing the ionization. The 
ratio of the space charge galvanometer reading 
to the secondary electron current then measures 
the relative ionization. Both galvanometer 
readings are taken at the same time, hence such 
slow fluctuations of the beam strength that the 
galvanometers can follow the variation will 
have no effect on the results. 

The curve of the ionization of argon by argon 
neutrals as a function of the square root of the 
energy is shown in Fig. 3. The ordinates are the 
ratios of the detector galvanometer reading to 
the secondary electron current corrected to give 
an approximate value of the cross section per cm* 
of gas at 1 mm Hg pressure. This cross section 
is obtained by first observing the sensitivity of 
the space charge detector to argon ions pro- 
duced by electron ionization. Since the efficiency 
for ionization by electrons is already known, 
the number of ions formed by the known elec- 
tron current can be calculated. The ratio of the 
detector galvanometer reading to this calcu- 
lated ion current is the sensitivity. This value 
of the sensitivity is then used to calculate the 
ion current produced by the neutral atom beam. 
Knowing the size of the neutral beam, one may 
readily find the cross section in units of Nie in 
cm?/cm* at 1-mm pressure. 

The ionization function for A in A is seen to 
fall with increasing velocity. Although the 
various curves obtained at different times do not 
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Fic. 4. Ionization of helium neutrals in helium. 


fall off at quite the same rate, they all have the 
same general shape. Figure 4 shows the curve 
of the ionization efficiency of helium neutrals 
in helium. The cross sections denoted by the 
ordinates are calculated in the same way as for 
argon. The curve shows a continuous rise but 
with a cross section much less than that for A. 

Figure 5 shows a similar curve for the ioniza- 
tion of Hz by hydrogen neutrals. Since in the 
ionization of hydrogen by electrons there are 
obtained three kinds of ions, it should be ex- 
pected that the neutral beam have several kinds 
of neutrals. Hogness and Lunn" show that 
chiefly H2+ and H;* are produced, hence H: and 
H; would be the components of the neutral 
beam. As the gun pressure is changed, the rela- 
tive amounts of the two ions change. However, 
no significant change in the shape of the ioniza- 
tion function is observed for the different gun 
pressures. 

For nitrogen, as in the case of hydrogen, there 
is more than one kind of particle in the neutral 
beam. The yield of N2+ and N* for different 
pressures has also been investigated by Hogness 
and Lunn.” They show that as the gas pressure 
increases, the ratio of N+/N2* increases. Figure 6 
shows the efficiency of ionization of nitrogen by 
nitrogen neutrals. It is to be noted that the 
curve has both a maximum and a minimum in 
the range studied. Since the ionization is caused 
by two kinds of particles, it should be expected 
that the experimental curve be composite of the 
ionization of nitrogen by Nz and N. When 


" T. R. Hogness and E. G. Lunn, Phys. Rev. 26, 44 


(1925). 
2% T. R. Hogness and E. G. Lunn, Phys. Rev. 26, 786 


(1925). 
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different gun pressures are used, the following 
is found: As the pressure is lowered, the maxi- 
mum of the curve becomes less pronounced with 
a slight shift toward the higher energy, and there 
is a faster rise at the high voltage end of the 
curve. This would indicate that the maximum 
belongs to ionization by N but no definite 
conclusions can be drawn without knowing how 
the neutralization process affects the relative 
composition of the beam. 

It should be noted that the change in ioniza- 
tion over the whole voltage range for N¢ is slight 
and therefore more subject to the masking effects 
of the errors inherent in the experimental method. 
However, of the fifty-odd sets of data for nitro- 
gen, all those which may be considered reliable— 
that is, those which show no change in sensi- 
tivity during the run, show a linearity in re- 
sponse of detector to various amounts of ioniza- 
tion, and have well-grouped readings—have a 
maximum and minimum. The position of the 
maximum varies from 2000 to 3000 ev, which 
can be accounted for in the variation of the 
gun pressure, as can be the different rates of 
rise at the high energy end of the curve. 


DISCUSSION OF RESULTS 


Although the ionization of argon by argon 
neutrals was previously known, this is the first 
definite evidence of the ionization of nitrogen, 
helium, and hydrogen by neutrals of their own 
gas. A previous work on argon by Varney™ 
shows a continuous rise from onset to 400 ev, the 
limit of his experiment, Since the ionization 
observed in the present experiments is decreasing 
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Fic. 5. Ionization of hydrogen neutrals in hydrogen. 


1% R. N. Varney, Phys. Rev. 50, 159 (1936). 


continuously with voltage in the range of 1000 
to 8000 volts, a maximum must lie somewhere 
between 400 and 1000 volts. The newly obtained 
curve for ionization of nitrogen actually shows 
such a maximum. 

The question of the sensitivity of the space 
charge detector to various ions traveling at 
different speeds is an important one in that 
it is conceivable that the drop of the ionization 
function for A at higher energies may be a 
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Fic. 6. Ionization of nitrogen neutrals in nitrogen. 


spurious apparatus defect. The following evi- 
dence on the matter exists."“"* In vacuum the 
sensitivity of the space charge detector is as 
high as 50,000 for very slow ions (less than 1 
volt). But the sensitivity falls off to a value as 
low as 10 when the ions have several hundred 
volts energy. With gas in the detector the sensi- 
tivity to slow ions is less, having a value for the 
ions formed by electron bombardment in this 
experiment of 5000 to 10,000. This is at a pres- 
sure of about 5X10-* mm. The sensitivity to 
faster ions with gas in the tube is not so well 
known. The high efficiency of the detector for the 
argon ions produced by neutral atoms in the 
high energy range of the present experiments im- 
plies either that little kinetic energy is trans- 
ferred in ionizing impacts or that the detector 
is retaining its sensitivity to faster ions. The 
questionable drop in the ionization function for 
A begins at neutral atom energies of about 1000 
volts. It seems highly unlikely that a transfer of 
kinetic energy from the beam particle to the 
newly formed ion should set in only at these 
high energies. Until definite measurements of 
the sensitivity of the detector to fast ions in the 


™R. N. Varney, Phys. Rev. 47, 483 (1935). 
1H. Karr, unpublished M.S. Thesis, Washington 
University. 
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presence of gas can be made, such a possibility 
can only be indicated. It does seem, though, that 
the observed experimental results are not caused 
by this variation of the detector sensitivity to 
the velocity of the ion produced. 

Since the size of the neutral beam is measured 
by the number of secondaries that are produced 
when it hits a plate after passing through the 
ionizing chamber, the drop in efficiency may in 
reality be caused by a rise in the number of 
secondaries produced per neutral as the energy 
of the neutral increases. As the number of 
secondaries produced by a particle as it strikes 
a metal surface is very much a function of 
surface conditions, numerical results of other 
workers are likely to be neither reliable nor 
applicable. Qualitative measurements made dur- 
ing the present experiments seem to indicate 
that in the mean velocity range, about one 
secondary electron is produced for every two 
positive ions striking the surface. Rostagni!® has 
examined the production of secondaries by both 
positive ions and neutral atoms in the energy 
range up to 600 volts. His results indicate that 
as the energy of the particle increases, the num- 
ber of secondaries per ion and per neutral in- 
creases but with a decreasing rate. Judging from 
a rough extrapolation of the A curve showing 
secondaries per neutral as a function of voltage, 
the mean value for the 1000- to 5000-volt range 
would seem to be about 0.4 secondary electron 
per impinging particle. For helium, the value is 
about 0.5. This value of 0.5 was also used for 
hydrogen and nitrogen in the calculation of their 
cross sections. If there should be a variation 
with energy in the number of secondaries pro- 
duced by the neutrals in the range of this experi- 
ment, it would, for the case of He and He, 
merely increase their rate of rise, and similarly 
for Ne. It is not conceivable that this change 
could cause a maximum and minimum in the 
curves. For A, an increase in secondary emission 
with voltage would produce an appreciable de- 


16 A. Rostagni, Zeits. f. Physik 88, 55 (1934). 
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crease in the rate of drop, but it would not seem 
sufficient to remove it entirely. 

It should be mentioned again that the calcula- 
tion of the cross section in absolute units is only 
approximate. It involves a knowledge both of the 
number of secondaries liberated by the neutrals 
in order to determine the beam intensity, and 
also of the secondaries liberated by electrons 
from this plate at the end of the beam path. 
This last is necessary since it is not possible to 
maintain the plate positive with respect to the 
space charge cylinder; hence instead of the 
primary, the secondary electron current is 
measured. Moreover, in using the apparatus for 
measurement of ionization by electrons, the gas 
flow through the gun must be shut off and the 
voltages reversed. It is possible that these 
readjustments may affect the sensitivity of the 
space charge detector, changing it from its 
value for the ionization by the neutrals, the 
measurement of which is done just preceding 
the electron measurement. 

Previous work by Rostagni!’ on ionization of 
A by A neutrals gives a cross section of 3 cm?/cm? 
at 1-mm pressure at 600 volts. The results here 
are then in general accord with such a value, 
indicating that the indirect method of calculat- 
ing the cross section is essentially correct. It is 
interesting to note that the ionization by He 
neutrals occurs well below the range character- 
istic of the alpha-particle. And also for H: 
neutrals in He, ionization is found at an energy 
much less than that for the proton. For ioniza- 
tion by electrons, it is found that A has the 
largest cross section and decreasing in order 
through Ne, He, to He. The efficiencies found 
here for ionization by their own neutrals are 
similarly arranged. 

The writer wishes to express his gratitude for 
the guidance and help of Professor R. N. Varney, 
under whose direction the research was done, 
and also for the valuable suggestions of Professor 
A. L. Hughes. 
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Formulas for the initial slope of the discharge voltage and return voltage curves in absorptive 
capacitors are given and it is shown that the discharge voltage curves are limited by two 
threshold curves, one of them being directly proportional to the slope of the other one. Some 
applications for the determination of the dielectric relaxation function, the dielectric loss, 
the time constant, and the electrical charge of absorptive capacitors are given. 


BONING! has recently established some 
¢ relations connecting the initial slope of the 
return voltage curves with the dielectric loss 
and other characteristic properties of absorptive 
capacitors. A brief account of some of his 
deductions and results can be found in a paper 
by Whitehead and Eager;? these authors have 
performed an experimental test of one of Béning’s 
relations. Béning’s analysis is based upon Max- 
well’s theory of the two-layer dielectric. It has 
been recognized that this theory yields qualita- 
tive results only—the behavior of real dielectrics 
is much too complex to be accounted for by so 
simple a model. Therefore Béning’s formulas 
claim an approximate validity only. Yet by 
the application of the principle of superposition 
there can be obtained a rigorous treatment of 
the discharge voltage and return voltage curves 
observed in absorptive capacitors and of the 
way these curves may find useful applications 
in the technique of measurement. The validity 
of the principle of superposition for capacitors 
with solid dielectrics seems to be sufficiently 
established on experimental as well as theoretical 
grounds. 


I. THE SUDDEN OPEN CIRCUIT 


Consider a capacitor forming part of an elec- 
trical network. At a reference time ‘=0 there 
shall be a sudden open circuit in the branch 
containing the capacitor so that the terminals 
of the capacitor become completely isolated one 
from another and from the rest of the network. 
After the opening of the branch, the total cur- 
rent across the capacitor must be zero. If the 


1 P, Boning, Zeits. f. tech. Physik 19, 241 (1938). 


neti B. Whitehead and G. S. Eager, J. App. Phys. 13, 43 
2). 


principle of superposition is valid, this current 
can be expressed in terms of the geometric 
capacitance C, the ohmic resistance R, and the 
dielectric relaxation function g(t). Hence* 


dU U 
t>0: J(t)=C—+— 
dt R 


‘dU 
J (1) 


The solution of (1) gives the voltage at the 
terminals of the capacitor after the sudden open 
circuit. 

With the introduction of a linear operator 
FLU] this equation can be written in the 
shortened form 


FLU)+i(t) =0. (2) 


During the period of time previous to the 
sudden open circuit, i.e., for ¢<0, the capacitor 
shall have been energized by a given voltage 
U(t). In an anomalous dielectric, every variation 
of electric stress produces a current which con- 
tinues flowing even after the stress has ceased. 
The term i(¢) denotes the sum of the currents 
created during the period of energizing and still 
continuing after the insulating of the terminals; 
i(t) is given as a function of U(t) as follows: 


°dU 
i= f (3) 


It is this term that causes the rise of voltage 
at the terminals of a capacitor which has been 
short-circuited temporarily. Equation (1) demon- 
strates that this voltage must be different from 


3B. Gross, Phys. Rev. 57, 57 (1940). 


383 


seem 
only 
f the 
trals | 
and 
rons 
le to 
» the 
the 
| 
s for . 
> gas 
| the 
hese 
its 
the q 
of 
| 
here 
alue, 
ulat- | 
It is 
He 
cter- 
He 
ergy 
1iza- 
the 
der 
und 
are 
q 
for 
ney, 
one, 
»ssor 
| 
|| 


384 B. GROSS 


zero so long as 1(t) is different from zero. U(t) can 
remain zero only if the short circuit has lasted 
so long that at its end i(¢) can be neglected. 

The term j/o' in Eq. (1) represents the absorp- 
tion current which surges in the dielectric in 
consequence of the voltage variations which 
occur after the sudden open circuit. 

Immediately before and after the sudden 
open circuit the voltage at the terminals must be 
the same, for it can never jump spontaneously. 
So for t=0 one obtains the condition of con- 
tinuity 

lim U(0—e)=lim U(0+e). (4) 


It may be stated that this equation is not 
valid for a circuit where an inductor is put in 
series with the capacitor; in this case there 
would be a voltage jump at ¢=0. 

It is easily seen that (1) can also be written in 


the form 
(dQ/dt)+U/R=0, (5) 


where 


‘dU 
Q=CU+ f f ¢(a)dodr. (6) 


With the introduction of the quantity Q, 
the equation of the absorptive capacitor has 
the same form as that of the non-absorptive 
capacitor. 

In the case of a leakfree capacitor (R= ~), no 
loss of electricity could occur after the open 
circuit because there would be no ohmic con- 
duction. The quantity of electricity held by the 
capacitor at the moment of the open circuit 
would remain constant so long as the terminals 
remain insulated. If in (5) we put U/R=0, we 
have Q=const.=(Q(0). In the case of a leakfree 
capacitor, Q would obey a law of conservation. 
We therefore identify Q(t) with the total amount 
of electricity stored at any moment by the 
capacitor. CU is the quantity of electricity 
contained by the geometric capacitance and 
J«' the quantity of electricity absorbed by the 
dielectric. 

Consideration should be given to the dis- 
continuous variations of the applied voltage. 
At the moment of a voltage jump, dU/dt be- 
comes infinite. In consequence, the integrals (1), 


(3), and (6) apparently become indeterminate. 
But it is easy to see what values are to be at- 
tributed to the integrals in these cases. Every 
voltage jump AU(r;) is followed by an absorp- 
tion current AUg(t—7;). In the present case, 
for every discontinuous variation of U/(t), the 
value of an integral of the type fdU/drg(t—1r)dr, 
as calculated by ordinary methods, must be 
increased by a term AU;g(t—7;). It is in this 
sense we are using the symbol /. It is easily seen 
that the definition of the integral to which we are 
led by empirical considerations coincides with 
the definition of Stieltjes’ integral. Indeed, as it 
has been pointed out to us by F. M. de Oliveira 
Castro, the treatment can be made rigorous 
from a mathematical point of view if one writes 
the principle of superposition in the form of 
Stieltjes’ integral 


U. 


Il. THE INITIAL SLOPE OF THE DISCHARGE 
VOLTAGE AND RETURN VOLTAGE CURVES 
AND THE MEASUREMENT OF THE 
RELAXATION FUNCTION 


For the limit ¢=0, i.e., just at the moment 
when the sudden open circuit occurs and the 
terminals become insulated, the integral in (1) 
disappears. U has still the value it possessed at 
the end of the energizing period according to (3). 


Thus 


dU U(O) 7(0) 


Now consider the curve of return voltage. The 
capacitor has been charged during a very long 
time under a constant voltage E and then 
short-circuited during an interval T. The voltage 
jump of magnitude — £, occurring at the begin- 
ning of the short circuit, is followed by an ab- 
sorption current which at an instant ¢ is given 
by i(t)=—Eg(t+T). U(0) is zero. After the 
insulating of the terminals, the voltage rises 
again. The initial slope of this return voltage 
curve is therefore 
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One can just as easily obtain the initial slope 
of the discharge voltage curve. In this case, 
before the terminals have become insulated the 
capacitor has been charged under a constant 
voltage E for a period 7. The voltage jump 
at the beginning of the energizing period, of 
magnitude +£, produces an absorption current 
given by i(t) = Eg(t+7). Furthermore, U(0) = E. 
After the terminals have been insulated, the 
voltage begins to drop steadily. The initial slope 
of this discharge voltage curve is given by the 
equation 


(— = (9) 


The initial slope of the discharge curve ob- 
tained after a complete charge (T= © ; o(T) =0) 
is expressed by —E/RC. Its determination 
enables one therefore to obtain the characteristic 
time constant RC of the capacitor. The time it 
takes to charge completely an absorptive capaci- 
tor is generally of the order of many hours, even 
for low loss materials like mica and sulfur. 

The initial slope of the discharge voltage curve 
observed after the capacitor has been charged 
during an extremely short time (T->0) gives a 
measure of the initial value of the relaxation 
function. 

It is interesting to note that at the beginning 
of the discharge the voltage at the terminals of 
an absorptive capacitor decreases more rapidly 
than does the voltage of that non-absorptive 
capacitor which has the same ohmic resistance 
and geometric capacitance as the first one. 
Afterwards the discharge of the first one becomes 
slower. The initial slopes of the discharge voltage 
and return voltage curves are related one to 
another in a very simple way. Their sum is con- 
stant and given by —E/RC. 

Once the value of the geometric capacitance is 
known, the measurement of the initial slopes of 
the return voltage or discharge voltage curves 
permits the calculation of the relaxation func- 
tion by the formulas (8) and (9). In many cases 
for the direct determination of this function 
obtained by recording the discharge current of a 
short-circuited capacitor, this method may be 
substituted advantageously. When concerned 


with large values of t, g(t) becomes so small that 
only an electrometric method can give satis- 
factory results. 

The dielectric loss, characterized by the 
tangent of the loss angle 6, can be split up into 
two parts, one due to the ohmic conductivity 
and another one due to the absorption current. 
The first part is given by 1/wRC. According to 
Schweidler* the second one can be calculated if 
the relaxation function is known. It is 


1/oR+ ¢(u) sin wudu 
tan 6= (10) 
c+f ¢g(u) cos wudu 


The relations (8) and (9) give the connection 
between the dielectric relaxation function and the 
initial slopes of the return voltage or discharge 
voltage curves. We are therefore able to corre- 
late the dielectric loss with these functions, 
substituting in (10) for g(u) its values obtained, 
respectively, from (8) or (9). The formulas so 
obtained do not depend on any model circuit 
but suppose only the validity of the principle 
of superposition. 

For very high frequencies, the peak value of 
the component of the current in phase with the 
applied voltage is given by A ¢(0)+A/R, and the 
peak value of the component in quadrature with 
the applied voltage is given by wAC where A is 
the peak value of the applied voltage. Then we 
have Béning’s formula 


1 
¢(0)+— 
R tana 1 


tan 6= = +-——, 
wC wE wRC 
where tan @ denotes the initial slope of the re- 
turn voltage curve observed after a short circuit 
of extremely small duration. The critical value 


of w, for which this relation becomes valid, can be 
inferred from the condition 


=¢(0). 


* E. v. Schweidler, Ann. d. Physik 24, 711 (1907). 
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Ill. THE THRESHOLD VALUES FOR THE DIS- 
CHARGE VOLTAGE AND RETURN VOLTAGE 
CURVES AND THE MEASUREMENT OF 
THE TIME CONSTANT 


The application of the relations (8) and (9) 
requires the knowledge of the time constant RC 
and the geometric capacitance C. There are 
many other cases where one needs to know RC, 
which is perhaps the most important quantity 
for characterizing the behavior of a capacitor. 


Fic. 1. 


A conventional laboratory method for measur- 
ing RC consists in observing the discharge volt- 
age curve of the capacitor. If the voltage during 
the time ¢ drops from V; to V2, RC is calculated 
by the well-known formula ¢/In (V;/V2). In 
absorptive capacitors this formula cannot be 
applied because there the simple exponential law, 
upon which it is based, is no longer valid. One 
of the first tasks of all theories concerning the 
general behavior of absorptive capacitors has to 
be the deduction of a sufficiently precise method 
for determining RC. The method referred to 
above does not fulfill this requirement about 
precision because it implies the derivation of an 
experimentally given curve. In the following, 
there are developed some simple and general 
laws which provide a method of the kind needed, 
and help at the same time to come to a closer 
understanding of the behavior of absorptive 
capacitors. 

The discharge voltage curve of an absorptive 
capacitor depends to a great extent on the dura- 
tion of the previous charging period, decreasing 
rapidly if this period has been short and de- 
creasing slowly if it has been long. A similar 
relation exists between the duration of the short 


circuit and the rise of the return voltage, which 
increases with decreasing duration of the short 
circuit. In both cases there exist threshold 
values. 

The slowest of all discharge voltage curves U,(t) 
is observed after the capacitor has been charged 
under a continuous voltage E during a time T 
so long that any further charging would no 
more increase the quantity of electricity held 
by this capacitor. There 


f i()=0 


and 


FLU, (11a) 
U.(0) =E. (11b) 


The fastest of all discharge voltage curves U;,(t) 
is observed if the voltage E has been applied 
during a time 7 so short that there has not yet 
been any absorption of electricity in the dielec- 
tric. The quantity of electricity held by the 
capacitor is still equal to EC. There 


f i(t)=E(t) 


and 
FLU, J+£E¢(t) =0, (12a) 


U,(0) =E. (12b) 


The fastest of all return voltage curves U,(t) is 
observed if the capacitor, after having been 
charged under the continuous voltage E during 
an extremely long time, is short-circuited during 
a time T so short that it has lost only the quan- 
tity of electricity CE contained by the geometric 
capacitance, whereas the quantity of electricity 
stored in the dielectric remains unchanged. 
There 


T 
f elt)dt=0; i(t)= —Ee(t) 


and 
FLU. ]—E¢(t) =0, (13a) 


U.(0) =0. (13b) 


To obtain a direct relationship between these 
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voltage curves, differentiate (11a). Noting that 


d 

— t—r)d 

dt J dr 
0 dr? dt 


and taking into account (9), we obtain 


dU. 


dU, 
(Rc ) =E. (15b) 
dt 


Comparing Eqs. (12) and (15) one arrives at 
the relationship between U, and U, 


(16) 


Comparing (11), (12), and (13) one observes*® 


U.= U,- U, (17) 
and therefore 


dU, 
U.= Uat+RC (18) 
t 


All discharge voltage curves are limited by two 
threshold curves, one of them being directly propor- 
tional to the slope of the other one. The difference 
of these curves gives the threshold for all the return 
voltage curves. (See Fig. 1.) 

The relations (16) and (18) are generalizations 
of the differential equation of the non-absorptive 
capacitor 


U=—RC(dU/dt), (19) 


5B. Gross, reference 3, Eq. (12). 


which is obtained as a particular case for 
U, = U, and, consequently, U.=0. 
Integrating (16), we have 


Cc 


RC= 
U,(0) U,(t) 


(20) 


This relationship seems to be very suitable for 
an experimental determination of RC and of 
R and C alone if the measurements are repeated 
with a loss-free capacitor of known capacity 
connected in parallel to the test capacitor. 

Closely related to the measurement of the 
geometric capacitance is the determination of the 
quantity of electricity held by the capacitor. 
For this purpose, Eq. (15) suggests a simple 
method. By integrating (15) from 0 to © we 
obtain 


1 
Q(0) J Udt. (21) 


1/R times the integral over the whole dis- 
charge voltage (or return voltage) curve gives 
directly the quantity of electricity the capacitor 
held at the beginning of these curves. On the 
other hand, Q could be calculated if one per- 
forms an integration over the absorption cur- 
rent. Because it is simpler to determine 


f Udt 
0 


than it is to determine 


J vit, 


the first method seems to be more suitable. 
I am greatly indebted to Professor E. L. da 
Fonseca Costa who made this study possible. 
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The problem of contacts between metals and between a metal and a semiconductor is treated 
classically with the help of the results of wave mechanical theory of electron energy states in 
solids. The potential and electron density distributions in the two bodies near the contact are 
discussed. The bodies are assumed to be in immediate contact. The problem of a body in 
vacuum and the problem of two bodies separated by a gap are discussed qualitatively. 


INTRODUCTION 


E know that when two bodies are in con- 

tact their thermodynamic potentials of 
electrons should be the same and that the con- 
tact potential difference is equal to the difference 
of work functions at T=0°K. We do not yet 
have a detailed picture of how the potential 
and electron density distributions near the 
surface of each body are affected by the contact. 
Frenkel! using the Thomas-Fermi method has 
discussed the situation at the surface of a metal 
in vacuum. On account of the approximate 
nature of the method he gets zero work function. 
Mott? has considered the contact between a 
metal and a semiconductor. He is mainly con- 
cerned with the situation in the semiconductor 
when certain boundary conditions at the con- 
tact are assumed. In a previous paper*® we have 
discussed qualitatively contacts between metals 
and between a metal and a _ semiconductor. 
We shall now discuss the problem quantitatively 
treating it classically with the help of the results 
of wave mechanical theory of electron energy 


states in solids. 


GENERAL CONSIDERATION 


The potential energy of an electron in a solid 
is periodic with the periodicity of the lattice, 
approaching large negative values near the 
nuclei at the lattice points. Bloch’s treatment 
shows that energy levels of the electrons in such 
a potential field are divided into separate bands 
with gaps of forbidden energy caused by the 
periodic nature of the field. Changing the aver- 
age potential simply shifts the whole energy 


1 i: Frenkel, Zeits. f. Physik 51, 232 (1928). 
2 N. F. Mott, Proc. Camb. Phil. Soc. 34, 568 (1938). 


3H. Y. Fan, Phys. Rev. 61, 365 (1942). 


spectrum up or down. If there is a slow variation 
of the average potential such that it changes 
little over many unit cells, then to the first 
approximation the whole energy spectrum with 
the top and bottom levels of each band simply 
follows the variation and the density distribu- 
tion of energy states in each band remains con- 
stant. We shall thus assume: (1) The top and 
bottom levels of each energy band follow the 
variation of the average potential. (2) The varia- 
tion of the average potential affects only the 
density of electrons in the conduction band, 
since the number of energy states and hence the 
number of electrons per unit volume in the 
completely filled bands remains constant. These 
assumptions result in great simplification. They 
are, however, rough approximations especially 
when the average potential changes very rapidly 
as at the surface or the contact with another 
body, problems with which we are concerned. 

We turn now to the question of electron 
density in the conduction band. The wave func- 
tion of the electrons is of the form y= e‘*-'u(r), 
the energy being a function of k. The number of 
energy states per unit volume in the energy 
interval is given by (ix*)Q, where Q 
is the volume in k space between two surfaces 
of constant energy corresponding to W,; and W2. 
W is a complicated function of k and can be 
determined only by a complete wave mechanical 
solution for the given solid. Therefore it is not 
easy to determine 2. We shall adopt the cus- 
tomary approximation 


(1) 


where E is the energy of the bottom level and m 
is the effective mass of electrons for the band 
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W= E+-—.- 
822m 
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under consideration, which can be appreciably 
larger than the ordinary mass of an electron. 
This expression for energy is of the form for 
electrons moving in a space of constant potential 
(—E/e). It must, however, be borne in mind 
that we use this expression only for finding 
approximately the number of energy states in a 
given energy interval. The actual average poten- 
tial in the solid is quite different from (—E/e) 
and we should not identify the two with each 
other when we are concerned with potential. 

According to Eq. (1) the distribution of 
energy states is 


1 
(W-—E)**dw. 
4r°L h? 


The number of electrons in the energy states 
between W and E is 


n= 


With constant average potential the average 
density of electrons with energy in the interval 
W to Eis 


1 
n-—| f (2) 
We shall treat our problem classically by extend- 
ing the application of Eq. (2) to the case of 
variable average potential. 


The electrons in the conduction band of 
metals are highly degenerate. We have 


1 
n-—| 3(¢—£) 


3h 


(¢—£)*. (3) 


E varies with the average potential and n 
varies with E according to this equation. 

In a semiconductor there are few electrons in 
the conduction band. The lower energy bands 
are full. According to present theory the few 
electrons in the conduction band come mainly 
from impurity atoms. In the semiconductors 
with normal decrease of resistance with increase 
of temperature, the electron energy level of the 
impurity atoms is situated in the energy gap 


below the conduction band. If the impurity is 
electropositive it can supply electrons to the 
conduction band and we have an efficacious 
semiconductor. With electronegative impurity 
atoms electrons from the lower band are ab- 
sorbed by them leaving holes in the band and 
we have a deficient semiconductor. In general 
the impurity atoms may supply electrons to or 
absorb electrons from not only the normal 
energy bands but also localized energy levels 
due to imperfection of the crystal.‘ This may be 
of great significance in the theory of semi- 
conductors. We shall, however, consider only 
efficacious semiconductors and shall use the 
simple model: the impurity atoms supply elec- 
trons only to the conduction band and in the 
normal state all electrons in the conduction band 
are supplied by the impurity atoms. Since there 
are few electrons in the conduction band we 
have exp[(W—£)/kT ]>1. Equation (2) gives 


n= 8/2 ¢(f—E)/kT 
3 
=n) exp[—(E—E)/kT], (4) 


where mp is the normal value of n for E= Eo. Let 
the density of impurity atoms be N with local- 
ized energy levels at AE below the bottom of the 
conduction band. The density of electrons in the 
impurity levels is 


1 


In the normal state with E= Ey the total density 
of electrons in the conduction band and the 
impurity levels is N. When E is different from 
E, the density of excess electrons is 


N 
n+ 


1 
X [mo exp[(Eo—AE—$)/kT] 
+n exp[ —(E—Eo)/kT] 
— N exp[(Eo—AE—$)/kT] 
Xexp[(E—Eo)/kT]. (5) 
‘B. R. A. Nijiboer, Proc. Phys. Soc. 51, 575 (1939) ; 


A. H. Wilson, “Semiconductors and metals," Cambridge 
physical tract. 


= 
4 
ion 
ges 
ith 
| 
yn- 
nd 
ia- 
id, | 
he | 
ey | 
ey | 
lly 
lly 
ier NV 
on 
of 
1s. 
| 
al j 
ot | 
1) | 
m 


390 H. Y. FAN 


This should be equal to zero when E= Ey 
mo exp[ (Eo —AE—£)/kT ]+n0 
= N exp[ (Ey) -—AE—£)/kT] 
exp[ (Eo AE 


Since m)<N the first term in the bracket of the 
right-hand side of Eq. (5) is small compared with 
the other two and can be neglected. The number 
of excess electrons is then 


No 
exp[(E—AE—¢)/kT]+1 
X[exp[ — (EZ —Eo)/kT ]—exp[(E— Eo)/kT] 


exp[ —(E—£o)/kT] 
—exp[(E—Eo)/kT]} (6) 


provided exp(E—AE—{)/kT<1. It should be 
noted that when E— &£, is too large, then if it is 
negative the number of electrons in the conduc- 
tion band may become so large that Eq. (4) 
based on the approximation of Maxwellian 
distribution ceases to hold and if it is positive the 
condition exp(E —AE—{)/kT<1 may no longer 
be true. 


EXCHANGE AND CORRELATION ENERGIES 


The energy expression (1) can be improved 
by taking into account the effect of exchange 
and correlation energies. When the electrons are 
perfectly free and can be considered as a Fermi 
gas filling up each energy state with two elec- 
trons of opposite spins up to a maximum energy 
level, the exchange energy is given by® 


k?—km? Rn—k 

—| 2+ In 

2 Rm Rm+k 
where k is the wave number of the electron under 
consideration and k, is that corresponding to 
the maximum energy. For the electron of maxi- 
mum energy this expression gives e*k»/7. 
Strictly this expression applies to electrons in a 
field free space. We shall follow the Thomas- 
Fermi-Dirac method® and extend its application 


5 P. A. M. Dirac, Proc. Camb. Phil. Soc. 26, 376 (1930) ; 
A. H. Wilson, The Theory of Metals. 
( 945) C. Slater and H. M. Krutter, Phys. Rev. 47, 559 
1 


to electrons in a variable potential field. The 
correlation energy is more difficult to estimate. 
Wigner’ has shown that in the limiting case 
when the electrons have no kinetic energy the 
correlation energy bears the ratio of 0.292 to 
0.458 to the total exchange energy. For the ordi- 
nary density of electrons in the conduction 
band of a metal the correlation energy is much 
smaller.’ We assume it to be equal to 25 percent 
of the exchange energy. For the maximum energy 
electron the sum of correlation and exchange 
energies is 1.25¢k»/7. 

The above discussion may be applied to the 
electrons in the conduction band of a metal, 
which can be considered approximately as a 
Fermi gas. We do not have to consider the 
electrons in the filled bands. Their density being 
considered constant, they contribute a constant 
amount to the exchange and correlation energies 
of an electron. The highest occupied energy level 
in the conduction band is now 


1.25e? h? 
Rm+ 


8r 


Rm”, (7) 


E- 


where E should follow the variation of the aver- 
age potential. The volume in k space between 
two constant energy surfaces corresponding to 
W(k=km) and W(k=0) is 4xk,,*/3. The density 
of electrons in the conduction band is therefore 


n= (82/3)km*. 
With the help of Eq. (7) we get 
(2.5e2m [6.25e4m? 
oh 


1/2) 3 
n +2m(¢—B) | . (8) 
In a semiconductor there are few electrons 
in the conduction band, which can be considered 
as a non-degenerate gas. The exchange and corre- 
lation energies are small and we shall not intro- 
duce such correction. 


POTENTIAL EQUATION 


We shall smear the positive charge of the 
nuclei and the negative charge of the electrons 
in the filled bands into a fixed uniform positive 
charge, the density of which is equal to the 


7E. Wigner, Phys. Rev. 46, 1002 (1934). 
8 E. Wigner, Phys. Rev. 46, 1002 (1934), curve in Fig. 7. 
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normal negative charge density of the conduc- 
tion electrons in the case of metals or equal to 
the sum of normal negative charge densities of 
the conduction electrons and the electrons in the 
impurity levels in the case of semiconductors. 

Consider first a metal. With the help of Eq. (8) 
the Poisson equation can be written as (for one 
dimensional problem) 


av Sa | 2.5e2m 
de 
6.25e*m? 2.5e7m 
= +2m(¢— 


6.25e4m? 


where V is the average potential inside the metal. 
It has been pointed out that E follows the varia- 
tion of V; i.e., dE/dx = —edV/dx. Let 


6.25e4m? 


+2m(¢—E). (10) 


J= 


Equation (9) can be rewritten 


(A+ ve) | 
h ¥ h Yo 
= (11) 


Integrating and using the boundary condition 
that dV/dx=0 where E= Eo, we get 


x 
— (yo? 
+ +3 Koyo? (12) 


The right-hand side of this equation can be 
expressed in terms of the electron density n. 
From Eq. (8) and (10) we have 


(13) 


where 


h 3 1 3h' 
v=( = —n. (14) 
2.5e°m/ 


Substituting (13) into (12) we get 


dy\2 
= N N23 
(15) 


Put N=No+AN and expand the right-hand 
side of (15) in terms of AN 


dy\? 
(=) AN? 
(16) 
On the other hand we have from (13) 
Ay = — AN 
— (17) 


With AN/N<1 both (16) and (17) converge 
rapidly. By taking the first term of each series 
and substituting (17) into (16) we get 


Now —1/3__ 


Ay*=PAy*, (18) 


where 


3 1/2 
(; 


h? \ -) (19) 


Equation (18) gives after integration 
AyaA'e#" (20) 


AVa—e#" =Ac*", (21) 
2me 


or 


For semiconductors with the help of (6) the 
Poisson equation is 
—— = (22) 
dx? 


Here we have assumed that the bottom level 
of the conduction band £ follows the variation 
of V, i.e., dE/dx= —edV/dx. We can therefore 
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write 
(23) 
= e? a 
dx? 
where 
2=(E—Eo)/kT (24) 
and 
K=4re'no/kT. (25) 


Integrating and using the boundary condition 
that dV/dx=0 where E= Ep, we get 


dz/dx = +(8K)! sinh(z/2). (26) 
Integrating again we have 
tanh(2/4) = tanh(—eAV/4kT) = Be#@4*2, (27) 


CONTACT BETWEEN TWO BODIES 


Our problem is to determine the charge 
density and potential distributions near the 
contact surface. We have a second-order differ- 
ential Eq. (9) or (22) for each of the two bodies. 
The constants of the first integration are de- 
termined by the boundary condition inside 
each body: dV/dx=0 where V= Vo. Two more 
conditions are required for determining the 
constants of second integration, constant A of 
Eq. (21) or constant B of Eq. (27). These are 
supplied by the requirement of continuity of V 
and its derivative at the contact. The latter is 
the condition for the two bodies being electrically 
neutral as a whole. Let body 1 extend from — » 
to 0 and body 2 from 0 to ~. We have 


dV, 


- (=) - . 
- (=) 


z=0 
If gi+q2=0, then 
dV, dV. 
dx z=0 dx z=0 
The condition of continuity of V is 
Vit(AVi)2-0= V2t(AV2)2—0 
or 
(A Vi A V2)2—0 
= (29) 


since {;={¢2 is the condition of equilibrium. 
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Let us consider first two metals in contact. 
If we use the approximate Eq. (18) then the 
condition (28) gives 


AV,/A V2= —1,/l,. (30) 


The minus sign is due to the fact that we have to 
choose the negative sign when taking the square 
root of the right-hand side of Eq. (18) for metal 
2 and the positive sign for metal 1. From (29) 
and (30) the values of AV; and AV, and hence 
the integration constants A for the two metals 
can be determined. For the solution of the prob- 
lem we need to know [¢—(—eV) ] and the normal 
density mo and effective mass of the conduction 
electrons. According to the approximate Eq. (21) 
at a distance x=1// from the contact, AV falls 
to 1/e of its value. Taking mp=2.62 X10” (the 
value for sodium) and the ordinary electron mass 
for m we get 1/1=0.29X10-§ cm. Thus the po- 
tential and charge density differ from their 
normal values only in the immediate neighbor- 
hood of the contact, within one atomic distance. 
Bethe’ has shown that the Thomas Fermi method 
used by Frenkel gives 1//=0.74X10-° cm. 
Within such short distance from the contact the 
electron density can, however, differ considerably 
from the normal value. If [¢—(—eV)] of the 
two metals differ by a fraction of an electron 
volt, AV will be of the order of a few tenths of a 
volt. With the numerical data as assumed 
above, for AV=0.1 volt, An is about twenty 
percent of mo. Owing to the fact that the charge 
due to An is concentrated over a short distance 
near the surface, a large amount of charge of 
opposite signs on the two sides of the contact 
is required to form a sufficient double layer for 
compensating the difference between [¢ —(—eV)] 
of the two metals. 

In the case of contact between a metal and a 
semiconductor the condition of continuity of 
dV/dx gives 
(31) 


—(8K)}! sinh =lAVm. 
e 2kT 


If eAV,/2kT<1 we may put sinh(eAV,/2kT) 
FeAV,/2kT. Then 


AV,/OVm= —1/(2K)}. (32) 


°H. Bethe, Handbuch der Physik 24/2, p. 420. 
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Metal Vacuum 
ectron poten- 
tial energy 
B 


average-eV, 
V 


Fic. 1. Conditions at boundary between a vacuum 
and metal surface. 


Putting m»=10" and T=300°K we get (2K)! 
= 1.18 10°. With the value of / estimated above 
we have |AV,/AVm| =2.96X10%. Thus AV, is 
negligible compared with AV,; all the difference 
between [¢—(—eV) ] of the two bodies is evened 
out by potential variation in the semiconductor. 
At T=300°K, eAV,/2kT =19.3AV, where AV, 
is in volts. It is possible that this factor is much 
larger than one, in which case we may put 


—eA V, 


eAV, l e(—AV ») 
+In 
2kT (2K)! 2kT 


Hence 


(33) 


with the same values for / and K used above. 
According to this equation |AV,/AV,|a1 at 
eAV,/2kT =10.3 corresponding to AV,=0.515 
volt. Thus for the case of small difference be- 
tween [¢—(—eV) ] of the two bodies most of the 
potential drop takes place in the semiconductor, 
whereas for a large difference the potential drop 
in the two substances may become comparable. 
Let us assume that eAV,/4kT <1; we may put 
tanheAV,/4kTHeAV,/4kT in Eq. (27), whence 


e 
AV,= ———B exp[ — (2K) !x]. 
4kT 


In a distance of 1/(2K)!, AV, drops to 1/e of its 
value. Comparing the estimated values of / 
and 1/(2K)! we see that the potential variation 
and space charge extend over a much wider 
region in a semiconductor than in a metal. The 
poorer the conductivity of the semiconductor 


or the smaller mo the larger will be 1/(2K)! and 
|AV./AV |. 

In the above discussion we took account of the 
exchange and correlation energies for metals 
by Eq. (7) and made them proportional to km 
or to the third power of the electron density. 
This applies strictly only to electrons in field 
free space with constant density. When the 
density varies these energies will depend not 
only upon the density at the given point but also 
upon the density distribution. Hence our ap- 
proximation is the better justified the smaller the 
density variation or when the difference between 
[¢—(—eV) ] of the two bodies in contact is small. 


CONTACT WITH VACUUM 


For the vacuum space in contact with a body 
[¢—-(—eV)] is equal to the negative work 
function of the substance, — ¢. 


]+¢ 


is of the order of several electron volts. Within a 
very short distance at the surface of the body 
the electron density changes from its normal 
value inside the substance to a negligible value. 
A more accurate method of evaluating the ex- . 
change and correlation energies is therefore 
necessary. In fact it can be easily shown that 
unless this is done we should get zero work 
function.!° 

Consider the surface surrounding the body, 
outside of which there are no electrons. If the 
body is not charged then according to Gauss’s law 
SdV/dnds=0. It follows that dV/dn=0 at this 
surface. Moreover, outside this surface V?V =0. 
Therefore the potential at this surface is equal 
to the potential at infinity. Our expression for 
electron density, Eq. (8), gives [f -—(—eV)]=0" 
for n=0. Hence at this surface {= —eV = —eV,, 
and the work function is equal to zero. 

The actual condition is represented by Fig. 1. 
At surface A {=—eV but the electron density 
has not yet dropped to zero, contrary to Eq. (8). 
The potential V continues to drop (—eV con- 
tinues to rise). At surface B the maximum energy 
level is equal to the potential energy which is the 
sum of —eV and the exchange and correlation 


0 Cf. H. Bethe, Handbuch der Physik 24/2, p. 417. 
For vacuum E in Eq. (8) is equal to (—eV). 
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(b) 


-eV, | el, 
(c) 


Fic. 2. Conditions at gap between two bodies. 


energies. Beyond B there are no electrons and V 
remains constant. The continuity of V and its 
derivative serve as boundary condition at the 
surface. It is clear from this discussion that for 
the treatment of this problem we have to wait 
for the development of a better method of taking 
care of the exchange and correlation energies. 
Bardeen” has given a wave mechanical calcula- 
tion of the work function for sodium. 

Consider two bodies separated by a gap. If 
there is no interchange of electrons the condition 
will be shown by Fig. 2a. Except at absolute 
zero temperature there will be an_ electron 
atmosphere in the gap, which is in equilibrium 
with electrons in the two bodies. This requires 
¢:=f2. Under the condition shown by Fig. 2a 
more electrons flow from body 2 to 1. As a 
result the two bodies become oppositely charged 
and a potential difference is set up across the gap 
bringing (:=fs. If the gap is not too narrow 


2 J. Bardeen, Phys. Rev. 49, 653 (1936). 


H. ¥. FAN 


the amount of charge acquired by each body 
will be too small to affect appreciably the elec- 
tron density and potential distributions at its 
surface. We may then assume the work functions 
to remain constant and the potential difference 
across the gap is gi— ¢2 (Fig. 2b). 

With decreasing distance the charges acquired 
by the two bodies increase in order to maintain 
the same potential difference across the gap. 
For parallel plates the surface charge density is 


¢2)/4nd. 


If the charge is concentrated within a layer of 
10-*-cm thickness the corresponding electron 
density will be of the order 


¢gi-—¢21 


4nd e 


With a difference of one volt between the work 
functions and a gap of 10-7 cm the charge ac- 
quired due to the contact gives rise to an addi- 
tional electron density of the order of 10. 
This approaches the order of magnitude of the 
normal density of electrons in metals. Therefore 
with such small gaps we can no longer assume 
that the potential and electron density distribu- 
tions at the surfaces remain unchanged (Fig. 2c). 
The treatment of such problem has difficulties 
explained above. 

The treatment developed in the previous 
paragraphs applies to two bodies in immediate 
contact. We have to know [¢—(—eV)] for 
each of the two bodies or the difference. This 
quantity for a substance is given by a wave 
mechanical solution of the electron energy 
states. It should be related to the work function 
and we would be able to determine it from the 
experimental data of the latter if we could 
derive the relation between the two. Experiments 
on contact phenomena such as contact rectifica- 
tion should yield useful information. 
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A new technique is described for quickly locating and measuring weak resonance points of a 
plate of piezoelectric material. The method involves driving the resonator with a frequency 
modulated carrier and demodulating and amplifying a voltage related to the motion of the 
resonator. The output of the amplifier is placed on one set of plates of an oscilloscope and the 
other plates are charged with the modulating voltage. The resulting pattern on the oscilloscope 
screen permits easy location and measurement of resonant frequencies. This arrangement was 
employed with a dynamical method to determine all elastic constants of beta-quartz. Edge 
effects were eliminated by the use of high harmonics as described by Atanasoff and Hart. 
Several hundred frequencies related to six modes of four cuts were measured, this number of 
modes being sufficient to permit the simultaneous solution of the resulting secular equations 
with an over determination of the five elastic constants. Table III gives the values obtained 
for these constants and also contains all the known results of other experimenters. 


HE experimental techniques available for 

evaluating the elastic properties of aelo- 
tropic or isotropic solids can be divided into two 
general classes. The static method consists in 
observing the deformation of loaded bars or 
plates. The constants measured in this way are 
termed isothermal since in the process the small 
temperature changes caused by flexing the 
specimen have ample opportunity to become 
equalized. In the dynamic method, on the 
other hand, the elastic properties are inferred 
from the frequency of standing waves set up in 
a given portion of the material. From knowledge 
of the frequency, dimensions, density, and mode 
of vibration, the elastic constants can be calcu- 
lated. In this case small temperature differences 
will exist between regions of nodes and loops. 
The constants will differ slightly from those 
computed by the static method and are termed 
adiabatic. 

The advantages of utilizing high order har- 
monics in calculations of adiabatic elastic 
constants from vibrating plates were first pointed 
out by Atanasoff and Hart.! They have clearly 
demonstrated that as the number of nodal 
planes between two opposite surfaces of a 
vibrating plate are increased, the perturbing 
effect of the edges has a diminishing influence 
on the frequency f/n, where f is the observed 
frequency and m the order of the frequency. 


1], V. Atanasoff and P. J. Hart, Phys. Rev. 59, 85-96 
(1941). 


The limiting frequency f/n then depends only 
on the thickness of the plate, the mode of 
vibration, and the elastic properties of the 
medium. 

Beta-quartz, also known as “high’’ quartz, 
exists between 573°C and 870°C. Although the 
static method of determining elastic constants 
becomes increasingly difficult at these tempera- 
tures, an experiment in which this technique is 
used on beta-quartz is described by Perrier and 
de Mandrot.? They have published data on 
Young’s modulus in directions parallel and 
perpendicular to the principal axis and also at 
an angle of 50° to the principal axis. This 
information is not sufficient to permit a calcula- 
tion of the five independent constants for this 
substance. 

The announcement of Osterberg and Cookson® 
that elastic vibrations of considerable vigor could 
be maintained in high quartz by means of its 
piezoelectric effect has been confirmed by a 
recent research* conducted in this laboratory. 
In this latter investigation, the resonant fre- 
quencies of a particularly accessible mode of 
vibration were examined. A measurement of 
these frequencies made possible a calculation of 
the Cys adiabatic elastic constant for this 


2 A. Perrier and B. de Mandrot, Memoires de la Societé 
Vaudoise des Sciences Naturelles, No. 7 (Imprimeries 
Reunies S. A., Lousanne, 1923). 

3H. Osterberg and J. W. Cookson, J. Frank. Inst. 220, 
361-371 (1935). 

4]. V. Atanasoff and E. W. Kammer, Phys. Rev. 59, 
97-99 (1941). 
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Fic. 1. Schematic diagram of experimental apparatus used 
in measuring crystal resonance frequencies. 


material. The success of this measurement 
suggested the feasibility of determining the 
remaining four constants with this experimental 
technique. 

The theoretical treatment of the present 
subject is the same as that used by Atanasoff 
and Hart in their work on alpha-quartz. Proper 
modifications must be made, however, for the 
different symmetry conditions encountered in 
beta-quartz which is classified as Ds while 
alpha-quartz is D3. 

The method of detecting and comparing the 
frequency of crystal resonances in the present 
experiment differs considerably from that of 
Atanasoff and Hart. In their technique a crystal 
was driven by an oscillator of adjustable fre- 
quency. The resonance points of the crystal 
were observed by using a vacuum-tube voltmeter 
and then the frequency of the oscillator was 
compared with a standard. The measurement of 
harmonic resonant frequencies of the crystal 
could thus be carried to a high degree of precision. 
However, with beta-quartz the piezoelectric 
effect is much weaker and the losses inherent in 
the crystal holder circuit components when at 
high temperatures made the deflections of the 
vacuum-tube voltmeter unreadably small. <A 
method permitting considerable amplification of 
the piezoelectric reaction in the circuit was 
essential to the success of a dynamic method in 
this application. The necessary amplification, as 
well as certain other advantages, was obtained 
in the technique to be described next, but with 
a reduction in precision of measurement to 
about one part in 5000. 

The experimental method used in the present 
investigation can, perhaps, best be explained by 
reference to the schematic block diagram Fig. 1 
and the circuit diagram Fig. 2. The output 
voltage of a high frequency oscillator is applied 


V. ATANASOFF 


RFo 


> 
Fic. 2. Apparatus for measuring resonance frequencies. 


to the electrodes EE between which the specimen 
is placed. The oscillator is frequency modulated 
by a 60-c.p.s. tone, causing the carrier frequency 
to sweep periodically across a band of con- 
trollable deviation (of the modulated frequency 
from its mean). If a resonance frequency of the 
crystal lies within the range of maximum and 
minimum values of frequency attained by this 
modulated carrier, it will absorb some energy 
from the field. The exact treatment of the 
response of a quartz resonator to a frequency 
modulated input is rather complicated ; however, 
at a low modulation frequency the following 
intuitive interpretation is asymptoticaily correct. 

Since the frequency of the applied field is 
continually changing, it remains only for an 
instant at the resonance frequency of the crystal. 
After this impulse, the quartz, having very low 
damping, will continue to vibrate at its natural 
harmonic frequency for some little time. During 
this time, two radiofrequencies are present in the 
circuit; the first, a continually changing one of 
constant amplitude due to the oscillator; the 
second, a fixed one due to the vibrating crystal 
whose amplitude is decreasing because of 
damping. These two frequencies are applied to 
the detector, and the result of the mixing and 
demodulation is a single varying-pitch audiobeat 


note. This audiotone is amplified and applied 


to the vertical plates of a cathode-ray tube 
whose sweep frequency is the same as that of the 
tone which modulated the carrier of the oscil- 
lator, in this case 60 c.p.s. 

A photograph of a typical trace obtained on 
the screen is reproduced in Fig. 3. The abscissa 
of this figure represents frequency increasing 
from left to right, the resonant frequency of the 
crystal being at the center of the pattern. The 


| 
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fact that the frequency deviation of the modu- 
lated signal can be controlled at will is a great 
convenience. A compressed pattern (large fre- 
quency deviation) for resonant points permits a 
rapid survey. A more precise examination of a 
single resonance point is achieved by means of 
an expanded pattern. 

A parallel LC circuit is used between the 
crystal holder and detector. The circuit details 
are shown in Fig. 2. This unit is adjusted to 
resonate at the center of the frequency band 
covered by the modulated carrier. The same 
network provides a convenient method for 
introducing a standard comparison frequency 
into the channel, and it also suppresses undesired 
harmonics in the standard oscillator output. To 
compare the frequency of the crystal resonant 
position on the fluorescent screen with a known 
standard, the signal from a standard oscillator 
is introduced by magnetic coupling and a pattern 
obtained somewhat similar to the one shown in 
Fig. 3. This trace is superposed on the screen 
with the crystal response. When the standard 
oscillator is tuned until the two minima of the 
patterns coincide, the standard oscillator is 
tuned to a resonant point of the crystal. 

Since the specimen between the electrodes EE 
must be at a temperature above 573°C if 
observations are to be made on beta-quartz, a 
special holder was devised. Nickel electrodes 
were used with iron rods as leads to the exterior 
of the furnace. A sheet of mica placed over the 
face of each electrode prevents contact with the 
specimen and reduces likelihood of breaking the 
plate due to non-uniform conduction of heat 
from the electrode. Unglazed porcelain plates 
and cylinders make up the balance of structural 
material. The holder assembly fits snugly in a 
resistance furnace heated by direct current to 
avoid induction signals being picked up by the 
detector. A Chromel-Alumel thermocouple was 
used to measure the temperature of the holder 
adjacent to the specimen. 

With the apparatus described one may easily 
measure the frequencies of resonance of a given 
specimen which fall in any selected frequency 
range. The problem of classifying these raw data 
is somewhat similar to the technique which must 
be employed in atomic spectra. If the crystal is 
examined in the proper frequency range one 


can select from the observed resonances a 
succession of frequencies which are nearly 
proportional to a number of successive odd 
integers. Such a series of frequencies can only 
accurately occur under the following conditions: 
1. That the crystal resonator must have parallel 
faces so placed that plane waves can travel 
between them. 2. That the frequency range must 
be high so as to make the odd integers large. 
It is natural to interpret such a series of fre- 
quencies as constituting the odd harmonics of a 
single mode of vibration but it is easy to see 
that these frequencies are not simply related to 
a fundamental frequency of a finite block but 
rather that they asymptotically approach the 
harmonics of a crystal plate with faces which 
are infinite continuations of the above-mentioned 
parallel pair. From these data one may easily 
estimate! the fundamental frequency of the 
infinite plate which yields the elastic constants 
according to a rather simple theory. When the 
parallel faces of a crystal resonator are separated 
by about one-half centimeter several members of 
each series (corresponding to these faces) occur 
in the octave below 10 mc. If the crystal plate 
has lateral dimensions of say two centimeters, 
it will be found that these frequencies already 
closely exhibit a constant difference and hence 
accurately determine the fundamental frequency 
of an infinite plate. 

Some of the plates used in this experiment 
were selected from those prepared by Atanasoff 
and Hart for their work on alpha-quartz. 
Several additional plates were made up with 
orientation not found among their collection. 


Fic. 3. Typical oscillograph of crystal response. 
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The material from which these new cuts were 
obtained was carefully examined for optical 
twinning by the use of polarized light. Moreover, 
the surfaces of the final specimens were etched 
with hydrofluoric acid in a further effort to 
detect twinning and non-uniformity of structure. 
The samples surviving these tests are described 
in Table I. 

Thus far the measurements relating to linear 
dimensions and angles of orientations have been 


TABLE I. The dimensions of the quartz plates 
in centimeters. 


Lateral dimen- 
sions in 


Orientation Thickness x's 
Y 30° 90° 0.4444 2.197 2.178 
R 30°, 51° 473’ 0.4500 2.176 2.222 
73° 30’, 59° 29’ 0.8081 1.469 2.214 
0° 45° 0.4442 2.227 2.200 


made on the specimen as alpha-quartz at room 
temperature. Since the material is to be heated 
to 600°C, expansion will change the dimensions 
of the plate as well as the orientation. Quartz 
expands at a different rate parallel to the optic 
axis as compared to that at any direction 
perpendicular to this axis, the latter expansion 
being the larger. These coefficients have been 
measured by LeChatelier' and also by Jay.‘ 
The values published by Jay will be used in the 
calculations to follow. 

To compute the thickness of the plate as 
beta-quartz from the thickness d at room 
temperature, use is made of the equation: 


The superscripts a and 8 on the direction 
cosines \ signify that they are the cosines of the 
angles between the normal of the plate and the 
optic or z axis as alpha- and beta-quartz, re- 
spectively. The value \° is not measureable and 
must be computed from the expansion factors. 
If » is the angle between the normal to the 
plate and the optic axis at room temperature, 
then the correction in radians to be added to 
this angle due to the effect of increased tempera- 


5 H. LeChatelier, Comptes rendus 108, 1045-1049 (1889). 
a 3) H. Jay, Proc. Roy. Soc. London A142, 237-247 
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ture will be very nearly: 
sin 2p. 


The quantities y, and y. are the changes in 
length per unit length over the temperature 
range considered for directions perpendicular or 
parallel to the optic axis, respectively. As a 
consequence of the change in angle yu, the 
direction cosines throughout the determinantal 
equations must be adjusted accordingly. Table 
II contains the corrected values for thickness of 
the cuts based on Jay’s results y,=0.00998, 
0.01733. 

The manner in which secular equations can 
be obtained for each plate of Table I has been 
explained in some detail by Atanasoff and Hart! 
in their paper on alpha-quartz. These secular 
equations contain the elastic constants as 
unknowns and for beta-quartz have the general 
form: 


| — K’, 0, C' i131 ° | 
| 90, — K?, 0 (=0 (1) 
C'rsn, 0, Chan 


The simultaneous solution of a set of equations 
of the type (1) begins with the evaluation of the 
characteristic value K? for the various plates. 


TaBLe II. Computed values of A? for the plates, with 
p=2.517 g/cc; temperature 600°C. 


Plate d(v) K? 
Orientation v thickness cm/sec. dynes/cm? 
of plate KC/sec. “d" cm «108 «100 
30°, 90° 417.49 4521 188.74 35.85 
0°, 44° 48’ 457.63 4502 206.02 42.70 
30°, 51° 35’ 
Mode I 436.60 4565 199.30 39.99 
Mode ITI 723.08 4565 330.08 109.6 
73° 30’, 
59° 18’ 
Mode I 240.1 .8205 197.0 39.06 
Mode ITI 407.6 .8205 334.4 112.5 


To do this it is necessary to know the density 
of quartz as well as the thickness of the plate at 
the temperature in question and finally, the 
frequency v of the mode of vibration being 
observed. The value 2.517 g/cm* used for the 
density of quartz at 600°C was taken from the 
work of Day, Sosman, and Hostetter.’ The 
thicknesses of the plates as beta-quartz after 


7 Day, Sosman, and Hostetter, Am. J. Sci. 37, 16 (1914). 
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making correction for expansion are listed in 
Table II. This table also shows the magnitude of 
K? as computed from the relation: K?=4pd?(v)*. 

In the secular Eq. (1) which applies to the 
Y plate, a factor appears, namely, (Cy,— K*) =0, 
which readily yields the constant C4, equal to 
35.85 X10" dynes/cm*. The variation from the 
earlier published value,‘ of 35.75 X 10" dynes/cm? 
is due to the use of Jay’s result of 0.01733 for 
the y, expansion coefficient for beta-quartz 
rather than the older value of 0.0161 estimated 
from the data given by LeChatelier. 

For the 0°, 44° 48’ plate, the constant C’122 
which occurs in the linear term of Eq. (1) is set 
equal to 42.7010" dynes/cm?. Thus: 


(0.70453)?C66+ (0.70967)?Ci4=42.70 X10". (2) 


When the above value for Cy is substituted, 
the elastic constant Cs. becomes: 


(3) (Cir — Ci2) = 49.7 X10" dynes/cm?. (3) 


In order to make certain that the value for K? 
used to find Ce¢ is the proper one, an examination 
of the constant Cy, is necessary and also the 
direction of the exciting electric field must be 
noted. For the 0°, 44° 48’ plate, the C4, constant 
which appears in C’j22 will be C3223. The piezo- 
electric constant é,4 for beta-quartz is the only 
one with strain subscripts ‘‘23.’’ Hence, the 
mode belonging to the linear term for this cut 
can only be excited by an alternating electric 
field with a component in the direction of the 
electric axis. Furthermore, the constants C4, 
which appear in the quadratic portion of Eq. (1) 
when applied to the 0°, 44° 48’ plate are of the 
form C3113. Since és is the only piezoelectric 
constant with the strain subscripts ‘13’ the 
modes belonging to the quadratic term require 
an alternating electric field with a component 
in the direction of the mechanical axis for 
excitation. 

Both modes predicted for the quadratic 
portion of the R secular equation were observed. 
Substitution for K* and the direction cosines 
together with the known value for Cy, yields 
the equations: 


(0.23700) — (53.825) Cir — (37.008) C33 
— (16.993) Cis— (23700) C213+8099=0; (4) 

(0.23700) C11Cs3— (11.039) Cy: — (10.097) C33 
— (16.993) Cis — (0.23700) C213 +165.6=0. (5) 


When the plate is taken from the crystal in a 
manner which prevents the normal from being 
perpendicular to any of the crystallographic 
axes, as is the case for the 73° 30’, 59° 18’ cut, 
it is not obvious which modes belong to the 
quadratic factor. In the determinantal equation 
both C3113 and C3223 are present in each factor 
so that all three modes could appear with the 
electric field applied along a lateral dimension. 
However, only two modes were observed in 


TABLE III. The adiabatic elastic constants and moduli 
for beta-quartz at a temperature 600°C. 


Perrier and 
Elastic Lawson de Mandrot 
constants (adiabatic) (isothermal) 
10" Elastic moduli X<10-" cm?/dyne 
dynes/cm? X10-" cm?/dyne cm?/dyne 
Cu =118.4 Su= 0.9257 Su =0,9345 
19.0 Siz = —0.0802 
Cis= 32.0 Sis = —0.252 Sis = —0.226 
C33 =107.0 1.085 S33 = 1.050 
Cu= 35.85* Su= 2.789 
S(45°) = 1.073 S(45°) = 1.067 
'S(50°) = 1,057 S(50°) =1.075 


* A value of 19.36 X10" dynes/cm? has been given for this constant by 
Osterberg and Cookson (reference 3). We regard this value as untenable 
and a similar opinion is held by Lawson (reference 8). 


this instance. If we know what Cos should be 
from Eq. (3), a substitution for each K? in the 
linear portion of Eq. (1) involving C’122: verifies 
that both observed modes belong to the quad- 
ratic factor. This yields two more equations 
relating Ci:, Cis, and C33, namely : 


(0.19267) C11C33 — (63.621) Cy, — (26.902) C35 
— (13.814) Cis— (0.19267) C?;3+8634=0; (6) 
(0.19267) — (9.2840) — (7.7460) Cos 
— (13.814) Cys — (0.19267)C?,3+125.5=0. (7) 


Equations (4), (5), (6), and (7) are more than 
enough to determine the three unknown con- 
stants they contain. The solution was achieved 
by a trial and error process of successive approxi- 
mations. The values for the elastic constants 
shown in Table III satisfy the four equations 
within 1 percent. 

It is of interest to compare these results with 
the isothermal measurements made by Perrier 
and de Mandrot. Their values for Young’s 
modulus perpendicular and parallel 
to the optic axis are expressed in kilograms per 
square millimeter, but in Table III these units 
are converted to dynes per square centimeter. 
After making a linear interpolation to the 
temperature 600°C, the reciprocals of E, and 
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Ey, namely, Si; and S33, become 0.934 and 
1.05X10-" cm?/dyne, respectively. The com- 
puted modulus for an angle of 50° with respect 
to the optic axis is 1.05X10-" cm?/dyne. This 
quantity was also measured directly by Perrier 
and de Mandrot. Again interpolating on their 
table to 600°C, the estimated value for this 
modulus becomes 1.07X10-" cm?/dyne. The 
accuracy of this latter interpolated value in 
particular is doubtful since the experimental 
points are widely separated in the region of 
600°C and the curve is known to take a rather 
sharp bend in this same region. Considering 
the experimental difficulties inherent in the 
static method when used at these temperatures, 
it is felt that the agreement between isothermal 
and adiabatic moduli is satisfactory, the dis- 
crepancies being of the order of a few percent. 
Lawson has made a dynamic measurement of 
the modulus in a direction of 45° to the optic 
axis. This was accomplished by setting up 
longitudinal vibrations in a quartz bar and 
computing Young’s modulus from the frequency 
of vibration, density, and length of the bar. 
Corrections were calculated for the perturbing 
effect of the finite dimensions of the practical 
bar on the theoretical frequency of a thin bar 
having the same length. In this manner Lawson 
arrives at the value 1.067 X10-" cm?/dyne for 
the 45° modulus while the data in Table III 
yield 1.073X10-" cm?/dyne. This agreement 
between these two adiabatic measurements taken 
in such different ways is gratifying. Using Perrier 
and de Mandrot’s measurements of S33 and Si 


8 A. W. Lawson, Phys. Rev. 59, 608 (1941). 
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together with the recently published‘ value for 
Sis, Lawson proceeds further to estimate the 
modulus obtaining —0.225 x 10-" cm?/dyne. 
A greater discrepancy appears here with the 
value listed in Table III. Since Lawson’s 
measurement of the 45° modulus agrees so well 
with the present set of constants the dependence 
of his result for S;3; on the work of Perrier and 
de Mandrot seems to be the principal source of 
error. On the other hand it should be noted that 
Lawson has made a slight error in his reduction 
of the results of Perrier and de Mandrot to 
c.g.s. units, 

Each determination of elastic constants by a 
dynamical method requires a knowledge of the 
mode of vibration. This knowledge can be 
obtained by an experimental examination of the 
vibrations’ or it can be inferred from theoretical 
considerations.’ The method of Atanasoff and 
Hart throws a greater burden on the frequency 
spectrum by making the identification of modes 
depend on it. Since the frequency spectrum can 
be determined with high accuracy and since 
theory enters only in a simple limiting way this 
method seems to possess some real advantages. 
We feel that the principal source of error in the 
present investigation is the difficulty of con- 
trolling and measuring the temperature of the 
crystal with sufficient accuracy over the con- 
siderable period of time necessary to make 
readings. The slight discrepancies noted by 
Atanasoff and Hart and explained by Lawson*® 
are less important. 


9A. W. Lawson, Phys. Rev. 59, 838 (1941); 62, 71 
(1942). 
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Letter to the Editor 


ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. The closing date for this department is the 
third of the month. Because of the late closing date for the sec- 
tion no proof can be shown to authors. The Board of Editors 
does not hold itself responsible for the opinions expressed by 
the correspondents. Communications should not in general ex- 
ceed 600 words in length. 


Correction: Investigation of the X-Radiation 
from Te’! (125 Days) by Critical Absorp- 
tion and Fluorescence 


(Phys. Rev. 62, 83 (1942)) 
R. D. O'NEAL AND GERTRUD SCHARFF-GOLDHABER 
Department of Physics, University of Illinois, Urbana Illinois 
September 6, 1942 


HROUGH an inadvertency, Fig. 1 of the above Letter 

to the Editor included four large dots which were not 

a part of the original drawing. The corrected figure is 
shown below. 
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